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1 Introduction 

The classical three-body problem studies the motion of a system with three 
point masses under the action of the Newtonian gravitational potential. Let 
&i = OPi, i = 1,2,3, be the position vectors of the points Pi with masses 
m.j > 0, with respect to a chosen inertial coordinate system for the Euclidean 
3-space M 3 . Then a motion of the three point masses will be described as a 
curve in the (unrestricted) configuration space, namely the Euclidean space R 9 
consisting of all triples (ai,a2,a3). 

Newton's equation of motion is the following second order system of three 
vector differential equations 



to,, a, 



dU iriimj TOiTO/c 



9a, rf 



{a.j -a t ) + — — (a fc -a 4 ) (1) 

' ij r ik 

where {i,j, k} = {1, 2, 3} , and 

toito 2 to 2 to 3 mim 3 

U = 1 1 , = ai-aj , (2) 

t-12 r 23 ri 3 

is the Newtonian potential function. In brief, the central problem is to under- 
stand both the geometry and the analysis of the solutions of the above system 
ijljl. where each solution curve (or trajectory) is uniquely determined by the 
initial positions and velocities of the point masses. 



1.1 The classical conservation laws 

The equations Q amount to solving a dynamical system in phase space of 
dimension 18. It is easily seen, however, that the classical 3-body problem (in 
fact, the n-body problem for all n) is invariant under Galilean transformations, 
namely the 10-dimensional Galilean group of 4-dimensional space-time R 3 x BL 
Accordingly, there are 10 first integrals or conservation laws, and they actually 
reduce the integration problem to one of order 18 — 10 = 8. These are the 
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conservation of linear momentum, angular momentum and total energy, and 
they are easily deduced as follows. First, by adding the three vector equations 
we have 

miai + m 2 a 2 + m 3 a 3 = 
and hence the center of mass has the uniform motion 

acM = — (miax + m 2 a 2 + m 3 a 3 ) = p + tv , m = J2 m i 
m 

where po and Vo are two constant vectors (and hence six scalar conservation 
laws) determined by the initial data. To make effective use of these we recall 
the Galilean principle of relativity of Newtonian mechanics, according to which 
one may choose an equivalent inertial frame of reference with origin at a^M and 
axes in the same directions as before (or rotated). This justifies using a center of 
mass inertial reference frame, which effectively reduces the actual configuration 
space to the following 6-dimensional Euclidean space 

M = K 6 C K 9 : m ^ = ( 3 ) 

Next, the vector 

$7 = miai x a.i + ™ 2 a 2 x a 2 + ™ 3 a 3 x a 3 (4) 

is the total angular momentum of the system. It varies covariantly with rota- 
tions of 3-space, and it is also constant along a trajectory since by 



f2 = J2 rma-i x Sj = J2 m * a i x ( J2 ~ |~( a i _ a i)) = 

i i jjii r ij 



This is the law of conservation of angular momentum, whose geometric signifi- 
cance lies much deeper than that of conservation of linear momentum. 
Finally, the total energy is defined to be 

h = T-U (5) 

where 

T = - {mx |ai| 2 + m 2 |a 2 | 2 + m 3 |a 3 | 2 ) (6) 

is the kinetic energy and — U is the potential energy. Straightforward differen- 
tiation using Newton's equation Q gives 

h = f-il = J2 rriihi ■ a 4 - Y, &i = J2&i- ( m * a * ~ ) = 

OAi dsn 

Hence, h is constant along a solution curve of Q and this is the law of conser- 
vation of energy. 
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1.2 Least action principles 



Newton's equation provides a characterization of the motion from the "differen- 
tial " viewpoint, but it is also useful to characterize the motion as a boundary 
value problem, namely we ask 

for a given pair {p, q} in the configuration space and time interval (t\, t 2 ), 
what are those trajectories T(t),ti <t <t 2 , with T(ti) — p, r(t 2 ) = q? 

The idea of seeking solutions of the above problem as the extremals of a 
variational principle applied to virtual motions dates back to the 17th century, 
inspired by the success of Format's principle of least time in geometric optics. 
Thus, a type of least action principle for classical mechanics was proposed al- 
ready by Leibniz, Euler and Maupertuis. It was, however, Lagrange who finally 
provided a precise mathematical formulation : 

Lagrange's least action principle. The solutions of the above boundary 
value problem are characterized by the variational principle of extremizing 
the action 

Jx[r] = jrdt (7) 

among all virtual motions between a given pair of points and with the 
same constant total energy h. 

It should be noted that time is allowed to vary in the above integral, that is, 
the limit of integration is not fixed. This awkwardness led Jacobi to reformulate 
the least action principle to the problem of determining the geodesies on a 
suitably defined Riemannian manifold (in modern terminology). In his famous 
lectures on mechanics UJ, Jacobi essentially introduced the notion of kinematic 
metric on the configuration space Mq, namely the kinetic energy expression JfjJ 
defines a mass dependent Euclidean metric 

ds 2 = 2Tdt 2 = m t {dx 2 + dy 2 + dz 2 ) (8) 

where a^ = [xi,yi, Zi) is the position vector of the point Pi with mass m^. 

On the other hand, T = U + h is also a function on Mo for a fixed energy 
level h, and this explains the following step: 

Jacobi's reformulation of Lagrange's least action principle. The action 
integral 

V2 Ji [r] = y/2 J Tdt = J Vfds = J VWThds = J ds h (9) 

is the arc-length of the virtual motion T in the metric space (Afh,rfsl), 
namely in the subspace 

M h = { P eM ;U{ P ) + h>Q} CM (10) 

with the squared arc-length element 

ds\ = (U + h)ds 2 (11) 
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Thus, according to Jacobi's " geometrization trick" , trajectories of Newton's 
equation are precisely the geodesies in the above Riemannian space (M/,, cfe^). 
This also demonstrates why Jacobi in his study of mechanics, in fact, anticipated 
the general notion of a Riemannian metric. For more information on these issues 
we refer to Liitzen[S]. 

On the other hand, in 1840 Hamilton formulated another least action prin- 
ciple, also inspired by the results of geometric optics. 

Hamilton's principle of least action. The solutions of the above boundary 
value problem are characterized by the variational principle of extremizing 
the action integral 



among all virtual motions T(t) between a given pair of points, for a fixed 
time interval [t±, t%] ■ 

In general, the validity of the " integral" viewpoint represented by any chosen 
variational principle is verified by calculating its infinitesimal limit, which must 
coincide with (or be equivalent to) Newton's equation. In the case of ||j5J| and 
HI 2(1 respectively, this amounts to the calculation of the geodesic equations of 
the metric dsh and the associated Euler-Lagrange equations of the above La- 
grangian function L, respectively. In both cases it is easily checked that these 
are equivalent to Newton's equation. 

1.3 An alternative geometric approach 

Traditionally, the three-body problem is usually studied in the framework of 
Hamiltonian mechanics, canonical transformations and symplectic geometry, 
based on the least action principle of Hamilton and the Hamilton- Jacobi theory. 
Moreover, the specific dynamics due to the Newtonian forces is usually assumed 
from the very beginning. Our present approach is, however, different from this, 
roughly for two major reasons : 

• Firstly, we focus attention on the purely kinematic properties of virtual 
three-body motions in a Riemannian geometric setting and in the frame- 
work of equivariant differential geometry, and 

• secondly, the Newtonian dynamics is introduced as the final step, involving 
geometric reduction and conformal modification of the kinematic Rieman- 
nian structure, based on the least action principle of Lagrange and Jacobi 
(cf. Jacobi[B], Lecture 6). 

Guided by the above program, the first author initiated studies in 1993 and 
was joined by the second author in 1994. The basic material, covering their work 
up to the winter of 1995, was presented in the two preprints [3], 0], and further 




ti 



L = T + U 



(12) 
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studies of the n-body problem continued in the following years. However, the 
two basic preprints were never published and, unfortunately, they have had a 
rather limited circulation in the mathematical community. 

On the other hand, in the recent years new and beautiful results on the three- 
body problem, and the more general n-body problem as well, have appeared in 
the literature, some of which are deeply related to the above geometric approach. 
This clearly suggests that new and unsolved problems along these lines are now 
becoming more feasible. To further stimulate this trend we propose hereby 
a review of the works 0] from 1994-95, and Chapter 1 -7 of the present 
monograph is, indeed, merely a faithful presentation of their actual contents. 

The exposition has been updated by some changes in notation and termi- 
nology, together with a restructuring of ideas and proofs in order to unify and 
enhance the readiblity of the presentation. For the convenience of the reader, 
central results are now formulated as main theorems, such as Theorem A, B, 
...,F, G presented in Section 2.2. 

In the final Chapter 8 we have included some additional and selected un- 
published material from 1995, mostly concerning the moduli curves of triple 
collision motions in the special case of energy h = (and, for technical reasons, 
the case of uniform mass distribution). The more general case is stated as an 
open problem in the last section. In this chapter we have done some explicit 
calculations which also serve as an illustration of how to apply the setting and 
the results up to Chapter 7. 

The three-body problem has a vast literature with many excellent papers. In 
particular, the theoretical analysis of the problem and its influence during the 
last century is overwhelming. Our interest in the problem started with a study 
of Siegel's monumental analysis of triple collisions involving clever applications 
of canonical transformations in the Hamiltonian setting. However, we were also 
astounded by the lack of basic geometric reasoning more directly linked to the 
kinematic geometry of three-body configurations, and apparently, this seemed 
to be typical in the more recent literature (such as Marchal's book), of which 
we had only superficial knowledge. Perhaps it was, after all, worthwhile having 
a closer look at the underlying geometric structure of mass triangles and their 
motions in 3-space ? 

With this ambition we started, hopefully with no prejudice due to the exist- 
ing literature, and this lead us to the purely kinematic study described in the 
first five chapters of this monograph, together with some preliminary investi- 
gations of the ensuing dynamics due to gravitational forces, most of which is 
presented in Chapter 7 and 8. 

Thus, in our "blindfolded" study during 1994-95 we deliberately avoided and 
did not consult any paper on the three-body problem (except Siegel's work). This 
explains why the list of references were almost empty, and we must apologize 
for that. The present list purposely reflects this former state of affairs, but now 
there are at least some relevant titles which were available prior to 1995 and 
which may be useful for the reader. Some of these are also standard references 
of historical interest. In retrospect, various topics and results discussed in this 
monograph are certainly more or less treated by the many authors who have 



7 



contributed to the rich literature in classical or celestial mechanics. Therefore, 
we also apologize to those authors who may feel that we have failed to make the 
appropriate reference to their work appearing before 1995. On the other hand, 
in the present monograph references newer than 1995 have not been considered. 

2 The basic setting and a presentation of the 
Main Theorems 

2.1 Basic notions and terminology 

Let R 3 denote Euclidean 3-space with the standard basis {i, j, k}. A three-body 
system consists of three labelled point masses (Pi,TOj) in 3-space, and its geo- 
metric model is the triangle with vertices Pi and the mass TOj > attached to 
Pi, which we shall refer to as an m-triangle ( or simply a triangle) or configura- 
tion. An m-triangle will also be identified with its triple (ai,a2,a3) of position 

vectors = OPi. These triples are usually denoted by boldface letters such as 
X,Y, but occasionally we also use the notation S, Si etc. It is tacitly assumed 
that a fixed mass distribution (mi, mi, m 3 ) is given, and it is normalized so that 
mi + m 2 + 777-3 = 1. 

Thus, the abundance of individual motions of three point masses jointly 
combine to a rich variety of virtual motions of m-triangles, through which the 
triangle changes its kinematic invariants such as size, shape and orientation 
(position) and velocity. This is our starting point for a systematic investigation 
of the kinematics of m-triangles, as a basis for a geometric approach to the 
three-body problem. 

In addition to previously defined quantities, set 

a.j = the central angle (at origin) opposite to the vertex Pj 
ujj = the (scalar) angular velocity of slj for planary motion 
Ij = m 3 \sLjl 2 , I = h+I 2 + h= p 2 (13) 
Ci = —mill + m.2/2 + 7773/3 etc. (cyclic permutation of indices) 

Tj = \m 3 \kj | 2 , T = Ti + T 2 + T 3 , ci. © 

flj = mj&jXkj, fl = fli + fl 2 + 3 , cf. I@J 

A = the area of A(Pi,P 2 , P 3 ), Ai = the area of A(0, P 2 , P 3 ) etc. 

where I (resp. Ij) is the total (resp. individual) polar moment of inertia, and 
similarly T, H (resp. Tj,£lj) denote kinetic energy and angular momentum as 
in Section 1.1. See Figure 1. 

Certain functions of the mass distribution appear frequently, so we introduce 
the notation 

to* = —(1 — rrij): the dual mass distribution, with ^""^ 777* = 1 (14) 
rhi = m 2 m 3 , m 2 = m 3 mi , m 3 = mim 2 
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and the basic elementary symmetric functions of the symbols m; are 

nii = 1, m — "^2 iTT-i, fn — ra\m2vaz (15) 

2.1.1 Vector algebra and kinematics in the Euclidean space M 

We will assume a center of mass reference frame as in Section 1.1, and hence an 
m-triangle X is a vector of the 6-dimensional Euclidean configuration space Mo, 
cf. (PI. The zero vector X = represents the one-point triangle (or the "triple 
collision" configuration), and we say X is collinear or is an eclipse configuration 
(resp. is non-degenerate) if the subspace n(X) C M 3 spanned by the position 
vectors a^ has dimension 1 (resp. 2). A virtual motion X(i) (or S(t)) is a time 
parametrized curve in Mo, assumed to be (piecewise) differentiable so that its 
kinetic energy © is defined. The size of an m-triangle X is naturally measured 
by the Euclidean length 

p=VJ=|X| (16) 

in Mo with the kinematic metric JHJ, equivalently given by the following inner 
product of Jacobi type 

X • Y = mia! • bi + m 2 a 2 • b 2 + m 3 a 3 • b 3 , (17) 

where X = (ai,a2,a 3 ), Y = (bi,b2,b 3 ). For convenience, we define 

wxX = (itixa 1 ,wxa2,uxa3), w ei 3 (18) 

X x Y = mi(aj x b*) 

and observe the general triple product identity 

w x X- Y = w • X x Y, weM 3 (19) 

The infinitesimal generators of the 5'0(3)-action on Mo are the rotational 
(or Killing) vector fields 

X -> w x X, u) eM 3 ~ so(3) 

of fixed angular velocity u). These vectors are tangential to the S'0(3)-orbits. 
Thus, at each X the tangent space TxMq ~ Mo has an orthogonal decom- 
position into vertical and horizontal vectors, where the vertical ones are the 
above Killing vectors lj x X and the horizontal vectors Y are characterized by 
X x Y = 0, due to l(T§J). 

For any virtual motion X(t) in Mo, the velocity vector at each time t has 
the above type of splitting, namely 

X = —X = X w + ± h = fw x X) + x' 1 (20) 
dt 

where u> = oj(i) is commonly referred to as the (instantaneous) angular velocity 
of the motion. Correspondingly, kinetic energy splits as the sum 



1 9 1 
T=-|u;xX| 2 + - 

2 1 1 2 



X' 1 
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of purely rotational and horizontal kinetic energy, respectively. The motion is 
called horizontal if the velocity is always horizontal. 

Using i|19|) we also deduce the following relationship between the angular 
momentum and angular velocity of a virtual motion, namely 

fJ = XxX = XxX"=Xx(u;xX) (22) 

Indeed, to each m-triangle X is associated the inertia operator 

E x : R 3 -> R 3 , w -> X x (w x X) (23) 

relating the two vectors a; and fi. This operator is invertible when X is non- 
degenerate, whereas £1 determines u) modulo a summand along the line n(X) 
when X is an eclipse configuration. In any case, the rotational velocity com- 
ponent w x X in 12U|) is uniquely determined by and X. Consequently, the 
motion is horizontal if and only if fi vanishes; in particular, such a motion must 
be planary (see Remark [331). 

The above inertia operator corresponds uniquely to the associated inertia 
tensor 

Bx(u, v) = (u x X) • (v x X) = ^m,(u x a,) ■ (v x aj ), u,veR 3 (24) 

which is a bilinear symmetric form on Euclidean 3-space. They are related by 
the identity 

B x (u,v) = I x (u) ■ v 

For example, they provide an orthonormal eigenframe for each m-triangle, and 
hence a moving eigenframe for a motion of m-triangles, see Theorem D and 
Section 3.4. 

2.1.2 Oriented m-triangles and their configuration space M 

Since triangles in 3-space can be oriented we propose the following " refinement" 
of the notion of an m-triangle. Define an oriented m-triangle to be a pair (X, n), 
where n GS 2 C R 3 is a unit vector perpendicular to II (X). In particular, a non- 
degenerate triangle can be oriented in two ways, namely we say the orientation 
is positive (resp. negative) if (a l5 a 2 ,n) is a right-handed (resp. left-handed) 
frame. 

Clearly, for an m-triangle motion the orientation may be chosen so that 
the normals n(t) vary continuously along the motion, and then the orientation 
changes to the opposite one as the motion passes (transversely) through an 
eclipse configuration. In the study of planary motions we will (tacitly) assume 
the plane to be the xy-plane, with unit normals n = ±k. 

Set M+ (resp. M_) to be the set of positively (resp. negatively) oriented 
m-triangles, together with the set E of oriented eclipse configurations, where 
the latter includes the 2-sphere of orientations of the one-point triangle X = 0. 
Then the configuration space of oriented m-triangles is the union 

M = M+ U M- C M X S 2 , E = M+nM- (25) 
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and a virtual motion of oriented m-triangles is a parametrized curve on M 

r(*) = (X(i),n(t)) (26) 

As a submanifold of Mo xK 3 ~ R 9 , M inherits a Riemannian structure and 
a natural isometric action of 50(3), also referred to as the congruence group. 
Let us have a closer look at this manifold and the two projection maps 

M m n S\ 

where 7Ti is a 2-fold covering over the non-degenerate m-triangles. On the other 
side, each (X, n) is 50(3)-equivalent to some (Y, k), where Y lies in the xy- 
plane K 2 and hence belongs to 

W 4 CM :Y] mi&i = , a; e M 2 (27) 

Then it is a useful observation (not mentioned in the 1994-95 preprints) that 
7r 2 is the 50 (3)-equivariant projection of a homogeneous 4-plane bundle 

R 4 -» M ~ 50(3) x so(2) R 4 -> 50(3)/50(2) = 5 2 (28) 

where the 4-space in (|27|l is the fiber over the unit normal k of the xy-plane, 
and 50(2) is the rotation group fixing the z-axis. 

In M is expressed as the 50(2)-orbit space of 50(3) x R 4 , where 
h G 50(2) acts by {g,Y) -> (gh^^hY), and the 50(2)-orbit of ( 5 ,Y) is 
identified with the oriented m-triangle (gY,(7k). This is an example of a well 
known "twisted product" construction. 

2.1.3 Congruence moduli space and shape space 

The orbit space 

M = M/S0(3) = M+ U M_ , E = M+ n Af_ (29) 

is the (congruence) moduli space, where the union corresponds to the 50(3)- 
invariant splitting 1251) of M, and the shape space is the subspace 

M* = Ml U Ml , E* = Ml n M* (30) 

corresponding to m-triangles of fixed size p = 1. Since an m-triangle (and its 
congruence classes) is scaled by the size function p > (|16fl . M = 0(M*) has 
a natural structure of a cone over M* . So, we may regard M as the union 
of two identical cones or "half-spaces" M± glued together along their common 
boundary E. This will be made more precise below. 

Let us first investigate the topology of the above spaces from a trigonometric 
viewpoint, using the quadratic form I|b8f) representing the squared area of m- 
triangles. The triple (Ji, 12,13) is, indeed, a complete congruence invariant for 
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(unoricnted) m-triangles. There is only one "half-space" in the unoriented case 
and we may express it as the following cone in {/, }-coordinate 3-space 

M± ~ {(h, h,h) | Ij > 0,Q(h,I 2 ,h) > 0} (31) 

where the corresponding shape space M± is cut out by the plane I1+I2+I3 = 1- 
The eclipse variety E, defined by the condition Q = 0, is the cone over E*. 

Now it is not difficult to see that Mjjl is topologically a closed 2-disk with 
E* as boundary circle, and consequently the full shape space 1)31) [) is a 2-sphere 
M* ~ 5 2 with a distinguished equator circle E* separating the two hemispheres 
Mjl. Note that the triple of central angles 

(01^2, "3), I] a i = 27r 

is a complete system of invariants for the shape of unoriented m-triangles and 
hence these angles also yield coordinates for each of the hemispheres. 

It follows that the full moduli space M is the cone over a 2-sphere and hence 
is homeomorphic to 3-space, 

M = C(M*) ~ C(5 2 ) = R 3 , (32) 

in such a way that E ~ R 2 is the coordinate plane z = 0, M + is the upper half 
space z > and M_ is the lower half-space z < 0. 

Finally, from the viewpoint of equivariant geometry, we observe that the 
pair M D M* is the S l O(3)-orbit space of the vector bundle and its sphere 
bundle, namely 

M* = (50(3) x so(2) S 3 )/SO(3) = S 3 /SO(2) ~ S 2 (33) 
M = (50(3) x so(2) M 4 )/50(3) = R 4 /50(2) = 0(5 3 /50(2)) m M 3 

For comparison reasons, if we only consider unoricnted m-triangles, then the 
corresponding calculation of the moduli space as an orbit space will yield the 
closed half-space 

M± = M /SO(3) = M /O(3) = M 4 /0(2) w R 3 ± (34) 



2.2 Statement of the Main Theorems 

In this summary we focus attention on six main topics, each of which is centered 
around one or two main theorems, labelled by A, B, CI, C2... 



2.2.1 Kinematic geometry of m-triangles and universal sphericality 

For a virtual 3-body motion with vanishing angular momentum, that is, a hor- 
izontal motion T(t) in M, the kinetic energy T depends solely on the moduli 
curve T(t) in M, namely in terms of the local coordinates (Ji, it is the 
following "differential" expression 

T =\j + ^Jq\ mih+ili+ii? - CiIiii+Ji+2 J (35) 
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where Q = 16m\m2m 2 /S. 2 . This is, indeed, a positive definite quadratic form on 
the tangent bundle of the moduli space M and thus naturally defines a kinematic 
Riemannian metric 

ds 2 = 2Tdt 2 (36) 

For a general virtual motion the same expression l|35|) is, in fact, obtained 
from T by removing the rotational kinetic energy. Therefore, by (|21|) . the hori- 
zontal kinetic energy 



~~ 2 



1 ' 



d =, , 
dt T{t) 



= T-T U (37) 



may well be referred to as the kinetic energy in the moduli space. 

Both the definition and the formula for the above metric (|36|l are depen- 
dent on the given mass distribution in a rather intricate manner. Therefore, it 
is a pleasant surprise that such a kinematically defined Riemannian structure 
(M,ds 2 ) turns out to be not only independent of the mass distribution, but it 
is, in fact, isometric to the Riemannian cone of the Euclidean sphere of radius 
1/2, namely 

Theorem A Let I = p 2 be the moment of inertia and M* be the subspace 
of M with I — I, and set da 2 = ds 2 |m* to be the restriction of the kinematic 
metric. Then 

ds 2 = dp 2 + p 2 da 2 (38) 
(M*,da 2 ) ~ S 3 (l)/U{l) = CP 1 ~ S* 2 (l/2) (39) 

where S 3 (l) is the 3-sphere of radius 1 and S 3 (l) — > CP 1 is the classical Hopf 
fibration. 

The surprising emergence of spherical symmetry in the kinematic Rieman- 
nian space [M, ds 2 ) for arbitrary mass distribution naturally brings in the clas- 
sical spherical geometry as a useful tool in the study of the three-body problem. 
We propose to call this fundamental fact the universal sphericality of the kine- 
matic geometry of m-triangles. 

The orientation reversing map (X, n) — > (X, — n) of oriented m-triangles 
induces an isometric involution of (M*,da 2 ) with E* as its fixed point set, 
namely the distinguished equator which divides M* into two hemispheres M±. 
On this circle lie the three points py representing the shape of the three types of 
binary collisions, cf. (|97|l . Indeed, their relative positions on the circle determine 
the mass distribution uniquely, sec Section 4.4 and l|136|) . 



2.2.2 Unique lifting property 

Theorem B To a given curve T(t) in the moduli space M\{0}, together with 
a given constant vector fl and initial configuration F(to), there exists a unique 
curve T(t) in M with T(t) as its moduli curve and with f2 as its conserved 
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angular momentum. Moreover, the curve T(t) can be computed in terms of the 
C 1 -data of f(t). 

Consider the orbit map ir : M — > M, and observe that SO(3) acts freely 
outside the sphere 7r _1 (0) ~ S 2 and defines a principal bundle 

7T : M \ 7r _1 (0) -> M \ {0} (40) 

In particular, above the half-spaces M± ~ M.± there are locally trivializing 
diffeomorphisms 

SO (3) x M± -> M± (41) 

Geometrically speaking, a motion of m-triangles can be represented by a 
time parametrized curve T(t) in M (or M±) which (locally) consists of two 
components, namely a moduli curve T(t) that records the change of size and 
shape of the oriented m-triangles, and a, position curve j(t) in SO(3) that records 
the change of position. The latter curve is, of course, constrained by the fixed 
angular momentum. 

In the special case of fl = 0, namely the horizontal lifting of moduli curves, 
the proof of Theorem B follows directly from the standard theory of principal 
G-bundles with a connection (cf. e.g. 0), applied to the above principal 50(3)- 
bundlc. Therefore, we shall rather focus on the general case with f2 ^ and 
present two different proofs. The first proof involves the inertia operator (|23|) . 
and the second proof is an application of Theorem D stated below. We refer to 
Section 5.2.2 

2.2.3 Angular velocities and kinematic Gauss-Bonnet formula 

Theorem CI For a planary motion of oriented m-triangles with normal vector 
k and angular momentum f2 = f2k, let u>i = (f> i be the individual (scalar) 
angular velocity of the position vector a^ = OPi . Then 

Wi = w? + y, i = 1,2,3 (42) 

where ujf is a "differential" expression purely at the moduli space level, namely 

u)± = —{1^0(2 — I2&3) etc. (cyclic permutation of indices) (43) 

= ~ Q ^[(C 3 h ~ C 2 / 2 )^ - (Ci + 2m 2 I 3 )i 2 (44) 

+ (C 1 +2m 3 I 2 )i 3 ] etc. 

Remark 1 The proof of Theorem CI holds for non-planary motions as well, 
that is, the plane n(t) of the m-triangle is time dependent. Then u>i stands 
for the (scalar) angular velocity of the velocity component of &i in II(i), and in 
formula OM) fl must be replaced by the normal component Q -n(t) (cf. Theorem 
D). We refer to Section 3.2.1. 
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We introduce the following three kinematic 1-forms on the moduli space 
M - {0} : 

01 = —{I^doL2 — I^da^) etc. (cyclic •permutation of indices) 

= ■= ^[{Czh - £2/2)^ - (Ci + 2m 2 h)dl 2 (45) 

8mim 2 m3A7 ij 

+ (Ci + 2m 3 I 2 )dI 3 ] etc. 

In fact, they are invariant under scaling and may therefore be regarded as 1- 
forms on the shape space M* via the canonical retraction M \ {0} — > M* . They 
share the basic property 

d@i=2dA, i = 1,2,3 (46) 

where dA is the area form of the 2-sphere M* ~ 5 2 (l/2). Evidently, the 1-forms 
have singularities on the eclipse circle E*. 

By suitably combining the kinematic 1-forms on appropriate regions on 
M* and applying Green's theorem, the kinematic Gauss-Bonnet version as de- 
scribed by the next theorem follows immediately from Theorem CI and (|46l) . 

Theorem C2 Let the shape curve of a piecewise differ entiable motion of 
oriented m-triangles with f2 = constitute the oriented boundary of a region D 
in M* ~ 5 2 (l/2). Then the total change of position of the triangle is a rotation 
of angle equal to twice the oriented area of D, namely 

A^ = / W ?dt= / e i = Jf2dA (47) 

to dD D 

The above type of integral (|47|) is an example of the geometric phase in 
the literature. Its value depends only on the shape curve and is independent 
of its parametrization. In the case of a planary motion with nonzero angular 
momentum, however, the total change of position in the above case (I47|l has an 
additional term called the dynamical phase, namely as a consequence of l|42f) 

A<t> i = ff2dA + t fjdt (48) 

D t 1 

2.2.4 Moving eigenframe and Euler equations for m-triangles 

Let (X, n) be a nondegenerate oriented m-triangle. Then we can choose eigen- 
vectors of the inertia tensor l|24|l which constitute a positive orthonormal frame 

(ui,u 2 ,n) €50(3), 

where Ui,u 2 lie in the plane II(X) and Ui x 112 = n. By definition, 

£ x (ui,ui) = Ai, Bx(u 2 ,u 2 ) = A 2 , Bx(ui,u 2 )=0 (49) 
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where the two eigenvalues Ai may be expressed (cf. Section 3.4) as 



{Ai,A 2 } = i(J± y/P - 16m 1 m 2 m 3 A 2 ) = -(lisinp), (50) 

using spherical polar coordinates (p, 9) on the 2-sphere M* centered at the 
north pole Af, where ip is the colatitude with ip = at the pole. The eigenvalue 
in the normal direction ±n is the largest eigenvalue 

A3 = Ai + A2 = I. 

To a continuous motion of oriented m-triangles we may choose such an eigen- 
frame 

ff(i)={ui(t),u 3 (t),n(t)} (51) 

varying continuously with the motion. In particular, t — > $(t) is also a parametrized 
curve in SO (3). 

Theorem D Let $(t) in i|51|) be a moving eigenframe attached to a differen- 
tiable motion T(t) of m-triangles, with Q as the conserved angular momentum. 
Then the triple of inner products 



gx = SI • ui, g-2 = n ■ u 2 

satisfy the following system of ODE, namely 



93 



n n 



(52) 



91 = 92 

92 = 91 

93 



9i92 



1 

A3" 
1 

~i 
,±_ 

'\2 



namely 






1 , 
A2 


it 

2 


■> cos ip 


^3 


u 

2 


■> cos ip 















(53) 



where the numbers Xi(t) are the eigenvalues of the inertia tensor of T(t) and 
depend solely on the moduli curve T(t) = (p(i),ip(t),0(t)). 



Corollary 2 It follows from 1^5 3\) that 31 (to) — 32 (£0) — at just one time to 
implies that gi(t) — g2(t) = for all time. Thus, such a motion is planary if 
and only if the angular momentum vector is perpendicular to the m-triangle at 
just one time to. 

Remark 3 The system ^5,S\) is the exact generalization of the classical Euler 
equations for a rigid body, see e.g. Arnold' 1], p. 14-3, where Mi and the fixed 
numbers Ii correspond to our gi and Aj, respectively. In $5<'jp the additional 
terms are due to the change of shape, and the system is singular where the 
motion passes through an eclipse configuration, say, with Ai = and hence also 
gi = 0. In particular, the eclipse takes place along a line perpendicular to Q. 
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The triple (31,32,53) is the coordinate vector, with respect to the moving 
eigenframe, of the constant vector J~2 = fik. It determines the position of the 
m-triangle, in particular its normal vector n, up to a rotation around the k-axis 
by a specific precession angle x(t)- This angle is calculated by quadrature from 
the formula 

= n-kxn = + cf . Section 5.2.1 (54) 

|kxn| 2 .9?+.9 2 2V Ai \ 2 ' V 1 

It follows that the m-triangle motion T(t) is largely described by two curves 
on the 2-sphere, namely the shape curve T*(t) = (<p(t), 0{t)) and the precession 
curve, that is, the curve traced out by the normal vector n(t). Thus, for the 
study of non-planar motions it is a basic problem to investigate the relationship 
between these two curves. 



2.2.5 The reduced Newton's equations 

Now, let us focus attention on the dynamics of three-body motions, namely the 
Newtonian equation of motion 0J. Such motions are, of course, a very special 
subclass of all the virtual three-body motions considered before. 

The fl-reduced Newton's equation is a second order system of ODE for the 
moduli curves of three-body motions with a fixed angular momentum vector fi, 
namely the three equations (with index i mod 3) 

t A rr f m< + 2m i+1 ■m i + 2m i+2 \ , 
7 * = 4T * - 13 + -3 7 * ( 55 ) 

\ 'M+l '»,*+2 / 



1 



(rrii+ili+i - m i+2 Ii+2) 



' i,i+2 



which are easily derived from the system (JJJ. However, the individual kinetic 
energy terms T% depend on fi, of course, but are otherwise expressed solely at 
the level of the moduli space M. 

The two cases of planary and non-planary motions differ substantially in 
complexity, so we will consider them separately. In the following two theo- 
rems we assume the initial position T(t ) is not a collinear configuration (since 
otherwise the given initial data will be incomplete). 

Theorem El A planary three-body motion T(t) is completely determined by 
its moduli curve T(t), initial position T(to) and angular momentum vector. The 
curve T(t) is a solution of the fl-reduced Newton's equation $55p . with kinetic 
energy terms 

+ (56) 
where oji is the i-th individual angular velocity Km . 

Conversely, each solution curve of this Q-reduced Newton's equation can be 
realized as the moduli curve of a three-body motion in the xy-plane, with a given 
initial position and normal vector 17k as the conserved angular momentum. 
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Remark 4 The above tt-reduced Newton's equation may, of course, be expressed 
purely in terms of the coordinates {Ij} or the mutual distances {rij}, see {64V - 
In fact, such an fl-reduced Newton's equation in terms of coordinates {rij} was 
derived by Lagrange^. We refer to Section 4-3.1 for another version in terms 
of spherical coordinates of M ~ R 3 as a cone over the 2-sphere. 

For the statement of the general (e.g. non-planary) version of the above 
theorem, let {ui,u 2 ,n} denote a continuous eigenframe associated with the 
motion T(t), and let 

(51,52,53), where g\ + g\ + g\ = n 2 , (57) 

be the coordinate vector of l~i relative to this frame, as in Theorem D. The 
individual kinetic energies depend on the components gk, more precisely, they 
split into a tangential and normal component 

T l = T[ + TV (58) 

where the tangential term T[ depends on the normal component 173 and T? 
depends on g\ and (72 ■ We refer to Section 5.1 and 5.2.3 for a precise description 
of these quantities. 

Theorem E2 A general three-body motion T(t) is completely determined by 
its moduli curve T(t), initial position T(to) and angular momentum vector ft. 
The curve T(t) is characterized by the ft-reduced Newton's equations \55]) with 
kinetic energy terms Ti $581) depending on the moving frame coordinates \5 r /\j of 
ft and are thus coupled with the Euler equations, namely the first order ODE 

15$). 



2.2.6 Reduction of the least action principles 

Here we will only consider planary three-body motions and state the associ- 
ated fi-reduced least action principles, whose extremals are precisely the moduli 
curves of those planary three-body motions with a fixed angular momentum 
fi = ilk. In this case the total kinetic energy 

- D, 2 

T = T + T LJ =T+ — (59) 

and the Lagrange function L = T + U are, in fact, defined at the level of the 
moduli space M, and therefore the two action integrals 

J i n =jTdt, J 2 ,n = J Ldt (60) 
r *i 

apply to moduli curves T(t). 

Theorem F The solution curves of the planary tt-reduced Newton's equation 
can be characterized as the extremal curves of Ji t Q (resp. J2,n) applied to curves 
r in M with fixed end points together with fixed energy h (resp. fixed time 
interval [ti, £2] )■ 
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2.2.7 Shape curves of triple collision trajectories 

In Chapter 8 we initiate a general study of the geometry of moduli curves in 
the vicinity of a triple collision. According to a classical result of Sundman and 
Siegel, towards the collision these curves T(t) approach a ray solution, which 
in the generic case has the shape ±po of an (oriented) equilateral triangle. In 
the simplest case of equal masses, ±p are the poles of the 2-sphere M* ~ S 2 , 
and hence the correponding shape curves approach one of the poles. Thus, it 
is natural to focus attention on the family of shape curves representing a triple 
collision with the limit shape of po . Here we state a theorem which is a simplified 
version of Theorem Gi stated in Section 8.1. 

Theorem G Assume uniform mass distribution and zero total energy, and 
consider the family & of arc-length parametrized shape curves T*(s),s > 0, 
representing 3-body motions with a triple collision at s = 0, say T*(0) is the 
north pole of S 2 . The family has the following properties : 

(i) There is a unique curve Tq q for each initial longitude direction 6q at the 
pole, and Tq q and ^g a+7T / 3 are congruent modulo a rotation of the sphere. 

(ii) The six meridians representing isosceles triangles belong to the family &. 
They divide the sphere into six congruent sectors of angular width 7r/3, and each 
curve Tq o stays within a sector, at least until the first eclipse (at the equator). 

(Hi) The curves Tg o (s) are analytic in s, with no singularity before the first 
eclipse, and T* 8o+7r {s) = T* 9o (-s). 

(iv) The sign of the curvature of the above shape curves is the same inside 
a sector, and the sign is the opposite in neigboring sectors. 

3 Basic geometric and kinematic invariants of 
m-triangles 

3.1 Ceva-type trigonometry 

In classical Greek geometry individual triangles - not their motions and kine- 
matic relations - are the geometric objects of basic importance. A triangle 
A(Pi, Pi, P3) is specified by its three vertices Pi, and its congruence properties 
involve the fundamental geometric concepts " side" , " angle" and " area" , whose 
relationships are expressed by trigonometric identities and congruence theorems. 
In our study, however, we are rather concerned with m-triangles, that is, a pos- 
itive mass m,i is attached to Pi . Thus it is natural and useful to reformulate the 
usual trigonometry into a kind of Cev a- trigonometry, depending on the given 
mass distribution. 

Let us first establish the following three identities (cf. (|13|) ) : 



Ceva-area law : Aj = 



rrijA 



(61) 



Ceva-sine law : 



sin on 2 A 



i = 1,2,3 



(62) 



mi ||ai|| ||ai|| ||a 2 || ||a 3 | 
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Ceva-cosine law : 2^/mim 3 - \J Iilj cos a& = — C% (63) 
By calculating cross products such as 

= ai x m j a j — Tn2&i x a 2 + ?7i 3 ai x a3 
= (2m 2 A 3 - 2m 3 A 2 )n 
the first law 1)61(1 follows directly, and then the sine law 1(62(1 follows : 

sina3 2 A 



2m 3 A = 2A 3 = |ai| |a 2 | sina 3 



m 3 a 3 | lla^l ||a 2 || ||a 3 | 



Furthermore, consider the "small triangle" with side vectors {m^ai}, say, with 
one vertex at the center of mass O and an adjacent side along OPj for some j. 
The triangle has outer angles and by applying the usual cosine law to it we 
deduce the cosine law (fS3\i . 

Next, by combining the usual cosine law and its Ceva version, the relation- 
ship between the mutual distances Si and the moments of inertia is 

(1^)7^ { i jjjfc} = { i j2 ,3} J (64) 
J mjirik 

from which we also deduce 

/ = V mi m rl = -^C x + -^C 2 + -^C 3 (65) 
* — ' m 2 m 3 m 3 mi mim 2 

where the first identity is known as Lagrange's formula for the total moment of 
inertia with respect to the center of mass, and m* are the dual masses 1)140. 
Finally, consider the "Heron" quadratic form 

H(a, b, c) = 2(ab + bc + ca) - (a 2 + b 2 + c 2 ) (66) 

and recall the classical Heron's formula for the area A 

ff( S 2 , S 2 , S 2 ) = 16A 2 (67) 

Set 

Q{Ix,h,h) = H(m 1 I 1 ,m 2 l2,m 3 I 3 ) = 2 mjmjljlj -^rrijlj (68) 

i<3 3 

= ^ c i G 3 = 4mim 2 Ii/2 - C 2 etc. 

i<3 

and consider again the "small triangle" with side vectors m^. On the one 
hand, its area A is related to A by 

4A 2 = m 2 m 2 |ai x a 2 | 2 = 4rn 2 m 2 A 3 = 4m 2 m 2 m 2 A 2 

and on the other hand, 16A 2 = Q(Ii, J 2 , / 3 ) by (|67|l and 168(1 . Consequently, 
we obtain the 

Ceva- Heron formula : Q(Ii, 7 2 , 13) — 16r7i 2 m 2 m 2 A 2 (69) 
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3.1.1 A simple torque formula 

As a simple application of the Ceva-area law (|fc>l|) we prove the following re- 
sult concerning the individual torques due to gravitational forces acting at the 
vertices Pj of a nondegenerate m-triangle A(P 1; P 2 , P3). 

Lemma 5 Let ti be the torque of the Newtonian gravitational forces at Pi,i = 
1,2, 3, with respect to the center of mass O. Then 

ti = £li — 2mi77i2TO3A(-g 3 )n , i mod 3 

T i,i+1 r M+2 

where n is the unit normal vector so that (ai , a2 , n) is a right-handed frame. 

Proof. Let F12 and F13 be the gravitational forces due to the mass points P2 
and P3 acting on Pi, namely 

m 1 m 2 , . miro 3 , . 

F12 = — 3 — (a 2 - ai), F13 = — 3 — (a 3 - ai) 

r 12 r 13 

Then, by definition of torque and the area law l|tjl|) 

, . mimo , > mim 3 . 
t x = a : x (F 12 + F 13 ) = — 3 — (ai x a 2 ) H 3 — (ai x a 3 ) 

r Y2 r 13 
,171x1712 A ^ ,171x1713 . ., 1 1 

= ( — 3 — 2m3A)n — ( — 3 — 2m2A)n = 2mim2m3A(- 3 - g— )n 

r 12 r 13 r 12 r 13 

and similarly at the other two vertices. ■ 

Corollary 6 Corollary ti = if and only if rx2 — fx3, and all ti — if and 
only if the triangle is regular (i.e. equilateral). 

3.2 Analysis of angular velocities and kinetic energies 

In the orthogonal splitting l|2U|) of the velocity of a virtual motion X(t) = 
(ai(t), a2(t), a3(i)), the horizontal component further splits into two summands 

X' 1 = X" + X ff = + X°" (70) 

representing the change of size and shape, respectively, and correspondingly the 
total kinetic energy splits as 

T = T U + T h = T U + (T p + T a ) = - \u x X| 2 + {^p 2 + T a ) (71) 

In this chapter we will show that T h actually equals the expression in l|35|l , and 
in particular it depends only on the velocity of the image curve in M. This will 
justify our definition of T h as the kinetic energy T of the moduli curve, hence 
also our definition of the kinematic Riemannian metric on M 

ds 2 = 2Tdt 2 = dp 2 + 2T a dt 2 (72) 
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Our first proof of Theorem A is by showing that the metric l|72[) actually trans- 
forms to the metric (|138|) . 

The differential expression (|35f) . as a function on the tangent bundle of M, 
is calculated by eliminating from T its dependence on the angular momentum, 
namely the rotational energy. In fact, it suffices to consider a class of virtual 
motions whose term is easy to calculate and hence eliminate. For this single 
purpose we could as well assume T u = from the outset and simply express T 
at the moduli space level. However, it is also illuminating to analyze the class 
of planary motions with a broader perspective. 

3.2.1 Kinematics of planary motions and proof of Theorem CI and 
C2 

We assume the motion takes place in the xy-plane and write 

fi = ilk, lj = wk 

In this case (1711 reads 

in 2 i 

T = T U ' +T P + T a = -— + -p 2 +T a (73) 

On the other hand, from the orthogonal decomposition of each ay into its rota- 
tional and radial component 

k = oJj(k x a,-) + — &j , where p 2 — Ij, (74) 
Pj 

the total kinetic energy also adds up to 

Therefore, by combining (|73|l and (|75(l the "intricate" energy term T CT , respon- 
sible for the change of shape, is given by 

r = -V^ + (iy^--p 2 )--- (76) 

Now, start from the above expression to express T a purely in terms of the 
individual moments of inertia Ij . 

Lemma 7 Letujj be the (scalar) angular velocity ofa.j. Then 

u>l = — {I^a 2 — + — etc. (cyclic permutation of indices) (77) 

ai = - H-Wi+^ + ft^) etc. (78) 

8mim2m3A 12 13 
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Proof. Set 9j to be the angle of &j with respect to a chosen reference direction 
in the plane. Then ujj = 9j and 

on = 8 i+2 - 9i+i, on = uj i+2 - LOi+i , i mod 3 

The total (scalar) angular momentum sums up to 

= X x (w x X) • k = Iuj 

Consequently, 

IiUi = - I 2 uj 2 - I3W3 = Q — I 2 a 3 + I 3 a 2 - (I 2 uj\ + hu)\) 

and this proves formula J77J. 

Next, by differentiating the Ceva-cosine formula (|63|) for a\ with respect to 
t and use the expression for sinai from the Ceva-sine formula l|t)2l) for a%, we 
obtain the formula i|78|) . ■ 

Substitution of the expressions (|78J) into l(77jl also leads to the following 
formula involving only Ij's, namely 

wi = Z TrfKCah - C 2 I 2 )^ - (d + 2m 2 I 3 )i 2 (79) 

+ (Ci + 2m 3 / 2 ) j 3 ] + - etc. 

and this completes the proof of Theorem CI. 

Furthermore, using either 1)77]! . 178(1 and the relation di + ct2 + A3 = 0, or 
using (|75|) directly, we calculate 



^(E(^r + A-,)| +(2 -if)W (+2 ) 

\i mod 3 / 

and finally by insertion into l|76[l we deduce the formula 

T ° = ^Jq( E rnil i+1 l i+2 if - Ciliii+Ji+z) (81) 

\i mod 3 / 

Consequently, the metric l|72l) on M may be written 

ds 2 = dp 2 + p 2 da 2 (82) 

where 

d<J 2 = ( ^2 mil l+ il i+2 dl 2 - C i I l dI l+ idI i+2 J 

^ \i mod 3 / 
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is the induced metric on the shape space M* = (I = 1). Indeed, the metric 
expression da 2 is a tensor on M* since it is invariant under scaling in M. 

On the other hand, on M* the relation /1+J2+I3 = 1 implies dli+dl 2 +dls = 
0, and therefore the above metric on M* can be restated as 



rf 2 = I{ HI + (1 - m 2 )I 2 ] dlf + [-If + (1 - mi)/i] d/ 2 5 
' -[2/ 1 7 2 -(l-m2)Ii-(l-mi)/2+m 3 ]d/id7 2 



(83) 



where Q* denotes the restriction of Q to M* and we have used the mass nor- 
malization mi + 772 2 + 777,3 = 1- 

Lemma 8 TTie area form of(M*,da 2 ) is 

dA = — 7=dh A dl 2 

Proof. As usual, the area form expresses as 

dA = VDdh A dl 2 

where 

D = _L { 7i/ 2 (l-mi-/i)(l-m 2 -7 2 ) 1 _L_ 

Q*2 I -I[2/ 1 /2-(l-m2)/i-(l-m 1 )/2+m 3 ] 2 J 4Q* 

is the determinant of the metric (|83J) . ■ 

Finally, we turn to the kinematic 1-forms (|45|l on M, whose definition is 
suggested by the expressions (|79|) for the individual angular velocities. Regarded 
as 1-forms on M * they are related to the area form by 



dOi = A da 2 — dl 2 A da 



= 1 f dh A (-2m 2 d/2 + ^dl 3 + ^dl x ) 1 
2v^\ - rf/ 2 A (-2777 3 d/ 3 + ^d/i + ^dl 2 ) J 

= 1(^ + ^3 + 2(7772 + 7.^ 

2VQ 7 h 
= -^=dl 1 A d/ 2 = 2d A 

This proves formula (|46|l and, as observed in Section 2.2.3, this also completes 
the proof of Theorem C2. 



3.2.2 A purely kinematic proof of Theorem A 

From the metric expression ()82|l it follows that M is a Riemannian cone over 
the shape space (M*,da 2 ), expressed in (|30J) as the union of two isometric disks 
along their common boundary circle E* . Both disks are parametrized by the 
region Q* > in the (2i, /2)-plane, where 

Q*{h,h) = Q{h,h,l-h-h), <Ii<l (84) 
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is the quadratic form l|68|) with 23 = 1 — I\ — I2. In the following we will 
describe our original calculations in 0] leading to the discovery of the universal 
sphericality . 

At first glance, the mass distribution {777;} is intricately involved in the 
formula l|83[l of da 2 , so we will focus attention on the mass dependent quadratic 
form Q* . The major step of the proof is, in fact, the algebraic approach of 
seeking better coordinates by transforming the metric tensor da 2 into a simpler 
one. The geometric proof using the Hopf bundle (see Section 3.2.3 below) is, in 
fact, our second proof. 

Intuitively, one expects that optimal simplicity and maximal symmetry is 
achieved by a suitable affine transformation of the (ii, / 2 )-plane which trans- 
forms the region Q* > into the unit disk and makes the metric more "trans- 
parent". This simple idea was, indeed, the key leading to such a remarkable 
coordinate transformation. 

As indicated in Figure 2, Q* = defines an ellipse which is tangent to the 
triple of lines given by I\ = 0, 1% = and Z3 = 1 — I\ — I2 = 0. It is easy to see 
that its center of symmetry is the point {m\, m^), so we first set 

h=h-m\, I 2 =I 2 -m* 2 (85) 

and obtain 

Q* = mim 2 m 3 - (1 - 777 2 ) 2 / 2 - (1 - mi) 2 I 2 - 2(7773 - mim 2 )hl2 
This suggests a rotation through the angle 

_ 1 1 2(mim 2 - m 3 ) 
2 (mi - to 2 )(1 + m 3 ) 

and new coordinates x, y defined by 

l\= x cos ip — y sin ip , I 2 = x sin ip a + y cos ip 

V 2 f (86) 



m2)(l + 777-3 ) cos2t/) 



m 2 )(l + 777,3) cos2V>, 



and we notice the identity 

(—(1 — 77i 2 ) 2 + (1 — 7711) 2 ) sin 2ip + 2(7773 — mim 2 ) cos 2ip = 



Then 
where 



7.2 



Q* = 777 l77l 2 m3 — fl 1 X 



Mi = gK 1 - m i) 2 + I 1 - m 2) 2 ) + \( m i 



(7773 — 777im 2 ) sin 2ip 

M 2 = ^((1 - ™i) 2 + (1 - "72) 2 ) - 5 (mi - 
(m 3 — m 1 777 2 ) sin 2^ 
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Thus, by setting 



x, y=\ v 



Mi V M 2 

the expression (|86|l transforms to 

Q* = mi m 2 m 3 (l - x 2 - y 2 ) (87) 
Therefore, the following combined transformation 



m 1 m 2 m 3 . /mim 2 m 3 „ 

Mi V M2 



. /TOiTO 2 m 3 m 1 m 2 m 3 

1 2 = smipQ* x + cos Vo\/ 2/ + m 2 

V Mi V M2 

will transform the formula l|83|) of dcr 2 into 

2 l(l-y 2 )cfo 2 + (l- : r 2 )dy 2 +2xydxdy 
4 1 - x 2 - y 2 

From here, we simply set 

x — sin ip cos 9 , y = sin(^sm# (90) 
which will transform l|89l) into the metric (|1(J7|) . This proves that 

(M*,d<r 2 ) ~ 5 2 (l/2). 

3.2.3 The Hopf flbration and a geometric proof of Theorem A 

The moduli space M is, by definition, an 50(3)-orbit space with the induced 
differential structure, and according to (|33[) it is also the orbit space 

M = M/SO(3) ~ R 4 /50(2) « R 3 = C(S 2 ) (91) 

of the orthogonal transformation group (50(2), R 4 ). As a quotient of a Rieman- 
nian space by a compact group of isometries M has the induced orbital distance 
metric which measures the distance between orbits in K 4 (or M). 

Let 5 3 = 5 3 (1) C K 4 be the unit sphere and recall the well known classical 
Hopf fibration, which in the above metric setting reads 

50(2) -> 5 3 -> S 3 /SO(2) = CP 1 ~ 5 2 (l/2), (92) 

where the projection is a Riemannian submersion and the quotient space is the 
round 2-sphere of radius 1/2. Combined with l|91|l we have an isometry 

M ~ M 4 /50(2) = 0(5 3 /50(2)) = 0(5 2 (l/2)) 
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of Riemannian cones over the 2-spherc M* ~ S 2 (l/2). The cone M is homco- 
morphic to K 3 , but they are only diffeomorphic away from the cone vertex (or 
base point O) which corresponds to the origin € R 3 . 

Finally, to complete the proof of Theorem A it remains to observe that the 
above orbital distance metric actually coincides with the kinematically defined 
one. The two metrics are, for example, determined by the kinetic energy they 
associate to "motions" in M . These are the image curves of virtual m-triangle 
motions X(i), which can always be chosen with vanishing angular momentum, 
namely they are horizontal (cf. Section 2.1.1). These motions are planar, say 
X(t) is a curve in K 4 C M . Horizontal curves are those perpendicular to the 
S*0(2)-orbits, and at a point X ^ the horizontal tangent vectors constitute 
the subspace W(X) ~ R 3 consisting of all Y such that X x Y = 0. 

Now, the orbital distance metric on M is defined by demanding the projec- 
tion 7r : M 4 — * M 4 /50(2) = M to be a Riemannian submersion, that is, that 
the tangent map dir takes H(X) isometrically to the tangent space of M at 
7r(X). Equivalently, the kinetic energy associated to a moduli curve is the same 
as the kinetic energy of a horizontal lifting. On the other hand, the kinematic 
metric (|36fl also associates to a moduli curve the kinetic energy T of a lifting 
with vanishing angular momentum. Consequently, the two metrics on M are 
identical. 

3.3 Linear motions of m-triangles 

According to Newton's inertia law, in a center of mass reference frame and in 
the absence of forces, the trajectory of the three-body system in the Euclidean 
configuration space Mo will be a geodesic, namely a linear motion 

6(t) = (l-t)S 1 +tS 2 (93) 

where 6% = (a l5 a 2 , a 3 ), S 2 — (bi,b 2 ,b 3 ) are appropriate m-triangles. Such 
motions are also characterized by having constant velocity (5 2 — Si), and the 
motion (|93[1 has constant angular momentum 

Q = Si x S 2 (94) 

Moreover, twice the action integral J7J) of the motion from Si to S 2 is the squared 
distance 

2jTdt= \Si-S 2 \ 2 = J2m l \& i -b l \ 2 (95) 

On the other hand, it is clear that the two m-triangles Si lie in a common 
plane if the vector (|94|l is zero, namely the linear motion l|93ll has vanishing 
angular momentum. Then the motions (|93|l provide a useful tool in analyzing 
the kinematic geometry of m-triangles since their moduli curves are exactly the 
geodesies in the moduli space (M,ds 2 ). Their shape curves will be arcs along 
great circles (geodesies) on the round sphere M* = S 2 . 

Let us collect some simple facts, assuming Si and S 2 are m-triangles in the 
xy-plane (with normal vector k) and = 0. 
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Example 9 If Si and 82 have the same orientation, then the distance between 
their congruence classes Si in M is 

distC5 1: h) =\5 1 -S 2 \ (96) 

Example 10 Consider two congruence classes 6±, 82 in M whose shapes 8* are 
different and not antipodal points on S 2 . It is not difficult to see that represen- 
tative m-triangles 8i can be chosen in exactly two ways, modulo a rotation of 
the xy-plane. The two choices are (81,82) and (81,-62) for suitable Si and 82- 

The shape curves of the corresponding linear motions ifffg)) . for < t < 1, 
are the two geodesic arcs T± between 8* and 8* 2 whose union is a great circle. 
Each of the shape curves extends (as t — ► ±00) to the whole circle, minus the 
limit point (Si T ^2)* as \t\ - * 00, which lies on the opposite arc T%. 

We also remark that the relative position of Si and 82 can be calculated from 
the line integrals of the kinematic 1- forms 0j along the geodesic arc T* ± , see 
Theorem CI and j45j) . 

For easy reference, the shape of the three types of binary collisions are the 
following three points on E* 

bij : h = - , Ij = , h = (1 - m k ) (97) 

1 - m k l-m k 

where by = bji represents an m-triangle (a 1 ,a 2 ,a 3 ) with a^ = &j and 1=1, 
and {i,j, k} = {1, 2, 3}. There are two more points on the sphere M* which are 
of kinematic importance, namely the north pole and south pole {AA,S}. Each 
pole is, of course, the geometric center of the corresponding pole AfJ, and we 
refer to Corollary 1 1 31 for an intrinsic geometric characterization of their shape. 

Lemma 11 The north pole (resp. south pole) is the shape of the positively 
(resp. negatively) oriented m-triangle whose normalized individual moments of 
inertia equal the dual masses, namely 

N (or S) : Ij = m* = \(l~ mj), 3 - 1, 2, 3 (98) 

Proof. Let So = (ai,a 2 ,a 3 ) be an m-triangle with Ij as in (|98|l . It suffices to 
show that the three points l|97() on the equator E* have the same distance in M* 
to the point 8q £ MJ . Let <5i = (b 1 ,b 2 ,b 3 ) be the (unit size) m-triangle with 
bi = V2ai and b 2 = b 3 = - j-g^- bi. It is easily checked that J~i = in 
that is, the linear motion between Sq and <5i has vanishing angular momentum. 
The image of 81 in M is the point 8\ = 623 on E* , and according to l|96l) 



l*o " *l| = VE TO ^ k- b i| 2 = V 2 - = 2 sin - (99) 

equals the distance between Sq and 8\ in M. Therefore, their (spherical) distance 
in M* is equal to ir/4. On the other hand, it is clear from the above calculation 
(or by symmetry) that similar choices of Si with <5J = bi 2 or b 3 i lead to the 
same distance ir/4. ■ 



28 



3.4 Eigenvalues and eigenframe of the inertia tensor 

The bilinear form £>x denned by (|24|) is identical to the well known inertia 
tensor in classical mechanics, for the special case of an m-triangle X viewed 
as a rigid body. This is useful in the kinematic study of non-planary motions 
of m-triangles. The geometric interpretation of the quadratic form is that it 
calculates the moment of inertia I u of the body with respect to the central axis 
through the vector a;, hence also the rotational kinetic energy due to the angular 
velocity u>, namely 

2T W = B x (w,w) = \u x X| 2 = H 2 I u 

By an eigenframe of X we mean an orthonormal basis in 3-space consisting 
of eigenvectors (i.e., along the principal axes) of £?x, whose eigenvalues are 
the associated moments of inertia. We will use the following notation for the 
eigenvalues and associated eigenframe, 

Ai < A 2 < A 3 «— > (ui,u 2 ,n), (100) 

where Ui and u 2 span the plane n(X) if the triangle is nondegenerate, whereas 
in the collinear case Ai = and {u2,n} can be any orthonormal basis of the 
normal plane n(X)^. 

Lemma 12 The eigenvalues of Bx. are related by 

Ai + A2 = A3 = /, A1A2 = 4mim2TO3A 2 , 

and consequently 

1 



Ai = -{I±y/P - 16TOim 2 m 3 A 2 ), i=l,2 (101) 

Proof. We consider the case that X is nondegenerate, and then I is the moment 
of inertia with respect to the normal direction. Let {ui, u 2 } be an orthonormal 
frame of n(X) consisting of eigenvectors of £>x , namely 

B x (ui, u x ) = Ai, B x (u 2 , u 2 ) = A 2 , B x (ui, u 2 ) = 

Then 

Ai + A 2 = | Ul xX| 2 + |u 2 xX| 2 = ^ mj {2 |a/ - (u r a,) 2 - (u 2 - aj ) 2 ) 
= J2 m i l a j| 2 = 1 

Set 

ai = fflnUi + ai2U 2 , a 2 = a 2 iUi + a 22 u 2 (102) 
Then on the one hand 

Bx(a 1 ,a 1 )Bx(a2,a2)-Bx(a 1 ,a 2 ) 2 = (103) 

2 

A1A2 



|ai x X| 2 • I a a x X| 2 - [( &1 x X)-(a 2 X X)] 2 = 



an ai2 

121 «22 
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where by the Ceva-area law (I61|) 

2 

an 012 

«21 0,22 

and on the other hand, 



= |ai x a 2 | = 4moA 2 



Sx(ai,ai) = m 2 |ai x a 2 | 2 + ra 3 |ai x a 3 | 2 = 4m 2 m 3 (m 2 + m 3 )A 2 
Bx(a 2 ,a 2 ) = 4mim 3 (mi +m 3 )A 2 (104) 
Sx(ai,a 2 ) = -4mim 2 m 3 A 2 

When the expressions l|104|) are substituted into l|103|) we obtain 

AiA 2 = 4mim 2 m 3 A 2 

and then formula (|101fl follows. ■ 

Corollary 13 The poles $U8j) are the shapes uniquely characterized by any of 
the two equivalent conditions : 

(i) Ai = A 2 (i.e. "umbilical" shape) 

(ii) the m-triangle attains the maximal area 

A max = 1 (105) 
4 v /mim 2 m 3 

among all m-triangles with the same moment of inertia I . 

Proof. Clearly, Ai = A 2 if and only if the area A is given by the formula of l|105|) . 
On the other hand, let us maximize the area function Q = CiC 2 + C 2 C 3 + C 3 Ci 
(cf. Q65[l. (|68|) ^ . using Lagrange's multiplier method subject to the constraint 
1=1. It follows that 

Ci = —miin* + nijm* + mum* k 

or equivalently Ij = m*, by l|13|l . ■ 

The following result will also be useful. Briefly, it says that a linear motion 
whose shape curve is a meridian arc from a pole to the equator, has a constant 
eigenframe. 

Lemma 14 Let So = (ai,a 2 ,a 3 ) be an m-triangle with the shape of a pole 
$981) . and let {ui,u 2 } be an orthonormal frame of the plane II(oo). Moreover, 
let Si = (b l7 b 2 ,b 3 ) be a degenerate m-triangle satisfying Sq x Si =0 and 
u 2 ■ bj = for all j . Then 

S t (ui,u 2 ) = 

holds along the linear motion S t = (1 — t)So + 15\, where B t is the inertia tensor 
of St, cf. \2J$ - Hence, {ui,u 2 } is an eigenframe for St for each t. 
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Proof. From the above Corollary it follows that 

#o(ui,u 2 ) = Bi(ui,u 2 ) = 
Moreover, by the assumptions 

(u 2 x S ) ■ (ui x Si) = X) m j( u 2 x a-j) ■ (ui x bj) 

j 

= E m 3 [(ui • u 2 )(aj • bj) - (u 2 • bj)(ui • a,-] = 

j 

Therefore 

B t (ui, u 2 ) = (ui x [(1 - t)5 + tSt]) ■ (u 2 x [(1 - t)5 + t«i]) 
= (l-t) 2 B (ui,u 2 )+< 2 B 1 (u 1 ,u 2 ) 
+ (1 - t)t[(u! x S ) ■ (u 2 x Si) ± (u 2 x (5 ) • (ui x Si)] 
= (1 -i)i(ui x u 2 ) • (S xSi) = 



4 The spherical representation of shape space 

M* 

By the spherical representation we refer to an identification of M* with a round 
2-sphere, with a distinguished (northern) hemisphere whose natural orien- 
tation induces the positive orientation of the equator E* and hence the (east- 
ward) direction of increasing longitude. We also assume the (cyclic) ordering 
&23, &3i, i>i2 of the three binary collision points (lying on E*) is in the positive 
direction. Finally, the correspondence should represent the kinematic geometry 
and hence is an isometry 

M*^S 2 {l/2) (106) 

which identifies each shape S* with a specific point on the sphere. We will 
develop methods enabling us to express the spherical coordinates in terms of 
intrinsic invariants of 5* ', and conversely. 

Let (r,9) denote polar coordinates on S 2 (l/2), where r is the polar distance 
which measures the spherical distance from 5* to the north pole Af G and 9 
is the longitude angle. For convenience, we also introduce spherical coordinates 
(ip,9), where the angle <p — 2r is the colatitude with if — at the north pole. 
In these coordinates the Riemannian metric of the sphere M* expresses as 

da 2 = dr 2 + ^ sm 2 (2r)d9 2 = ^(d<p 2 + sin 2 tp d9 2 ) (107) 

0<r<|, 0<92<7T, 0<6»<2tt 

Remark 15 The choice of the zero meridian 9 = is a matter of convenience, 
and until further notice our convention is that it passes through b 23 . Indeed, only 
longitude differences (9 — 9') is an intrinsic property of shapes of m-triangles, 
see Section 4-2 and 4-3- 
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4.1 Geometric interpretation of the polar distance r 

Let S be a nonzero positively oriented m-triangle, that is, r < ir/4. We will 
investigate the relationship between the polar distance r and the geometric 
invariants of 5. Let 5* be the intersection point between the equator E* and 
the meridian passing through 6*. Choose unit size representatives 

So = (ai,a 2 ,a 3 ), S ± = (b x b 2 ,b 3 ) 

of the pole Af and d* with Sq x Si = 0, and observe that the above meridian is 
the shape curve of the linear motion <$ t = (1 — t)So +tSi, < t < 1. Henceforth, 
let t be the unique value such that 8* t — 5* . 

Let C be the cone surface (cf. also Definition I42|) in M spanned by the rays 
through points on the above meridian between Af and 5*. C is isometric to a 
Euclidean sector of angular width n/4, see Figure 3. The cord distance between 
Af and 5 \ is the number in (32} . On the other hand, the cord distance between 
Af and 5t can be computed in two different ways, namely 

— 7T 7T 7T 7T 

dist(Af, 6t) — 2t sin — = sin — — cos — tan(— — r) 
8 8 8 8 

Hence, 

* = ~(l-cot~tan(~-r)) (108) 

and moreover, 

cos 5 

dist(0,S t ) = tt^- (109) 

cos(| - r) 

Next, let us compute the eigenvalues of B t = Bs t , namely 

Ai = B^ux,^), X' 2 = B t (u 2 ,u 2 ), 

where by Lemma H^l we have chosen an orthonormal frame {ui, u 2 } of the plane 
n((5 ) with u 2 • hj — for all j. It follows that 

A' 1 = Bt(u 1 ,u 1 ) = \ui x 5 t \ 2 = |ui x (l-t)<5 | 2 

= (l-t^SoCux.m) = ^(i-0 2 

Therefore, since 5* and 5 t differ by the scaling factor (|109|) . the eigenvalues of 
are 

cos f , 1 . . , cos f 
A 1= 8 V x = - 2 8 110 

cos(| — r) 2 cos(| — r) 

and A 2 = 1 — Ai. 

Finally, we can use (|108|l . (|110|) and the identity 

AiA 2 = 4m im 2 m3A 2 

to solve for r as a function of the area A of S* . We state the final result as 
follows : 



32 



Lemma 16 The polar distance r for an arbitrary given positively oriented m- 
triangle S is given by the formula 



cos(2r) = VTOiTO 2 m 3 y (HI) 

where rrij = 1 and A (resp. I) is the area (resp. moment of inertia) of 8. 
Equivalcntly, by (|1U1|) there is the formula 

sin2r= j|Ai-A 2 | (112) 



4.2 Geometric interpretation of the longitude angle 6 

Let Sq be an oriented m-triangle whose shape <5q is the pole J\f or S (i.e . r = 
or 7r/2). Recall from Lemma ITU it is possible to deform So, through a linear 
motion with zero angular momentum, to the shape of any given degenerate m- 
triangle. Then the shape curve will be the meridian from the pole to a point SI 
on the equator circle E*. Moreover, the line spanned by the final configuration 
Si is uniquely determined by So, and there is a constant eigenframe throughout 
the deformation. 

Now, let us consider two points 6*, <5 2 011 E* and seek an interpretation of 
their spherical distance in M* = S 2 (l/2). 

Theorem 17 Let Sq — (ai,a 2 ,a 3 ) be an m-triangle, of maximal area for a 
fixed moment of inertia, and let Si , <5 2 be degenerate (but nonzero ) m-triangles 
satisfying the vanishing angular momentum condition 

S a x Si = So x S 2 = 

for the linear motions from Sq to Si, i — 1,2. Then the angle ip between the lines 
spanned by Si and <5 2 is equal to the distance between the associated points S\ 
and <5 2 in the shape space M* , namely 

ip = \\e 2 -6i\ 

Proof. We may assume all m-triangles are confined to the xy-plane, the shape 
Sq is the north pole and 5* has longitude angle 9j,i — 1,2. Consider the piece- 
wise linear motion whose associated shape curve is the spherical triangle D in 
M* with vertices <5J,(5 2 ,<5o. Starting from Si, the motion passes successively 
through the m-triangles SiS' 2 ,Sq, S[. Here 5' 2 is congruent to 82 and is situated 
in the same line as Si, whereas S 1 is congruent to Si but is actually situated in 
a line making the angle ip with the original line. 

Next, let us apply the Gauss-Bonnet formula Ij47(l to the region D, thus 
obtaining an equality between twice the area of D and the angle tj). Finally, we 
simply combine this with the fact that, as a geodesic triangle on the sphere of 
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radius 1/2, the area of D equals one quarter of its angle at 8q, namely the angle 

\02-0l\. ■ 

In general, it turns out that the longitude angle 9 of an m-triangle is de- 
termined by the relative position and size of the eigenframe and normalized 
area A/ 1 respectively. To make this relationship precise, let 8 = (ai,a 2 ,a 3 ) 
be a non-degenerate, positively oriented m-triangle in the xy-plane with normal 
vector 

ai x a 2 _ 
|ai x a 2 | 

and assume the shape 8* is not the pole Af. 

Theorem 18 Let {u 1; u 2 } be an eigenframe of 8, where \i\ is the eigenvector 
of the inertia tensor B$ associated with the smallest eigenvalue Ai, and let ip i 
be the (oriented) angle from &i to u^. Then the following identity 

9 l + sinw l + sin<^ 

tan — = tan?/>i = cottp 2 (113) 

2 cos if cos (p 

relates the angle ip i to the spherical coordinates (tp, 9) of the shape 8* on the 
2-sphere M* , where tp is the colatitude and 9 is the longitude (eastward, with 
9 = at b 23 ). 

Remark 19 The above formula holds for any choice of eigenframe since ip 1 
changes by ±7r if Ui is replaced by — Ui. 

Proof. The first step of the proof is to derive a formula which expresses ipi 
solely in terms of intrinsic invariants of the m-triangle, together with a simple 
recipe for calculating this angle. Then, by applying the Gauss-Bonnet formula 
(cf. Theorem C2) we shall deduce formula l|113|) . 

Since i/) 2 = ipx ± ir/2 we need only prove the first identity in (|113f) . First of 
all, in order to have the angle ip 1 uniquely defined we must specify the choice 
of eigenframe. Namely, let {u l5 u 2 } be a positive frame and hence U! xu 2 = k, 
and moreover, we assume ii! chosen so that 

-2<^i<2 ( 1U ) 

Let {ei,e 2 } be the orthonormal frame derived from {ai,a 2 } by the Gram- 
Schmidt algorithm, with ei = ai/|ai|. Using the expressions (|1(J4|) it is not 
difficult to deduce 

A 2 

A = Bsfei^i) = 4raim 2 m 3 (m2 + m 3 ) — 

B = B s {e 1 ,e 2 ) = ((m 2 + m 3 )C 3 - 2mim 2 h)— (115) 
C = Bs(e 2 ,e 2 ) = I — A 
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By writing 

Ui = cos^ej + sin^ea (H6) 
U2 = sin ^x e i + cos "0i e 2 

and inserting these expressions into Bs(u\, u 2 ) — 0, we deduce the formula 

tan2^ 1 = -^ = ^ 7 (117) 

and similarly 

Ai - A 2 = Bsiux,^) - B s {u 2 ,u 2 ) = {A- C) cos 2t/> x + 2£sm2<0 1 (118) 

How is V'i determined from (|117|1 and l|118|) ? Assume first B = 0. With the 
normalization J^toi = = 1, this happens when 

7i(m 3 - mim 2 ) + h(l - mi) 2 = (1 - mi)m 3 

For example, with uniform mass distribution this holds for the isosceles triangle 
with | as | = | &3 1 . Note that {ei,e 2 } is also a positive eigenframe, and the 
eigenvalues A^ equals A and C. Since Ai < A 2 , by assumption, we deduce from 
(|118|) that A > C implies V>i = §, and A < C implies ip 1 = 0. 

Next, assume B ^ 0. By combining l|117|) and (|118fl we eliminate 008 2?/^ 
and obtain the expression 

Ar - A 2 = ( A - C f + AB2 sin2 ^ (119) 
iti 

Consequently, sin 2^ has the opposite sign of B, namely 

< Vi < ~ if B < 0, - | < Vi < if S > (120) 

Finally, we turn to the proof of formula (|113|) . On the sphere M* , let <5J be 
the intersection point of E* and the meridian passing through 5* , and consider 
the spherical triangle on M* with vertices b23, 5*, S* , whose area is denoted by 
A. The right angle at the vertex <5J has adjacent edges of length 

101 7T , 7T 

s = — < — and r, 

2 ~ 2 4 

and by applying the spherical sine law and area formula to the magnified triangle 
on S' 2 (l) with area A = 4A, we have 

A sin2ssin(? -2r) 

tan — = 



2 1 + cos 2s + cos(f - 2r) + cos 2s cos(f - 2r) 

sin 2s sin(f- - 2r) 7r 

tanstan( r) 



1 + cos 2s 1 + cos(f - 2r) 4 
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On the other hand, applying the Gauss-Bonnet formula (Theorem C2) to 
the triangle in Figure 4, it is not difficult to see that the total rotation of the 
vector ai is through the angle hence 

A 

±V X = 2A = - 

and consequently 

6 TT 

tantp 1 = ±tan-tan(— — r) (121) 

We claim that the sign to be used in l|121fl is -1. This can be seen by 
considering the situation where 6* lies between 623 and b^i, by observing that 
4>i decreases as 9 increases (i.e. when 6* approaches b3i). This completes the 
proof of formula l|113|) . ■ 

4.3 Intrinsic form of the spherical representation 

We will focus attention on the "inverse" of the correspondence (|106fl . namely 

The first step in this direction was, in fact, our kinematic proof of Theorem A 
(cf. Section 3.2.2). Namely, by substituting (|§U)l into and considering the 
special case of 7713 = m\m%, it is easy to verify that the expressions l|88|l simplify 
to 

I\ = m*(l + sin (/3 cos 0), I2 — m^{\ + sin sin 0) (122) 

These formulas are, indeed, a special case of a general intrinsic description of 
the spherical representation, purely in terms of geometric concept. 

The longitude distance between 5* = (ip 1 ,9i) and 5% — (^21^2) is, by defi- 
nition, the angle \0i — 62] mod 2ir. Then it is easy to check that l|122|) can be 
stated as 

Ii — m*(l + sin ip cos 6 i), i = 1,2,3 (123) 

where Q\ (resp. 62 or #3) is the longitude distance between 6* — (>p,0) and the 
binary collision point 623 (resp. b^i or bi 2 ). 

On the other hand, consider the three distance functions on M* = S 2 (l/2) 

<j l = <Ti(6*) = dist(S*, b i+M+2 ) (i mod3) (124) 

which measure the (spherical) distances to the points by. By the spherical 
cosine law applied to the triangle with vertices Af, by, 5 , it follows that 

cos2c7i = s'm(pcos8i, £=1,2,3 

and consequently (|123f) has the invariant form 

h = m*(l + cos2cr 2 ), i = 1,2,3 (125) 
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This may be stated as 



cos2ct 4 = — - 1 = — L cf. ra (126) 

m* m* 



or equivalently 



Thus, in order to establish 1)125(1 or 1)123(1 as a general formula it suffices to 
verify formula ((127(1 in general. Again, the basic idea we use is to construct a 
suitable linear motion in Mq, as we did in Chapter 3, namely we consider the 
linear motion whose shape curve is the (shortest) geodesic from 6* to by. We 
shall calculate the length of this curve in the following way. 

Let 5 = (ai,a2,a3) be a given m-triangle, normalized with 1=1, and 
consider the linear motion of vanishing angular momentum 



„, , , . „ , 1 — mi nil mi , 

S(t) = (1 - t)S + t{J—±a. 1 ,- ai , - ax) 

V h y/[l-mi)Ii y/(l-mi)h 

between S = 5(0) and the normalized m-triangle Si = (5(1) with the shape 
Sj. = E>23- Its shape curve is the desired geodesic. 

The length \S — 5i\ of the segment in Mq from 5 to <5i is also the length cti 
of the chord in M between 6* and §*, see (|96f) and Figure 4. By applying the 
Ceva-cosine law (|63H to the calculation of inner products • a.j we arrive at the 
expression 

a 2 , = \5 - 5i\ 2 = 2(1 - ; /! 



1 — mi 



and consequently 



o • CT 1 / J l 

ai = 2 arcsm — 



2 V 1 - mi 

4.4 The reduced Newton's equation in spherical coordi- 
nates 

Let us utilize the structure of M as a cone over the 2-sphere to express the 
reduced Newton's equation of Theorem El in terms of the spherical coordinate 
system (p,<p,ff), where p = \fl measures the distance from the base point O 
(cone vertex). The relationship between coordinate functions and 
{p, ip, 8} is expressed by ifMt and (|123[l . thus enabling us to transform the equa- 
tion to a system purely in terms of p, ip, 9. This change of variable is, however, 
rather messy, but an equivalent system can be worked out in several ways. For 
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example, we obtain the following system of three 2. order equations 



(i) 



n 1 

p+- (U + 2h) 



P P 



{Hi) 



(m) 




(128) 



valid for planary three-body motions with a fixed energy level h. The angular 
momentum f2 constant is not implicit in these equations since it is an integration 
constant defined by the initial value problem. In fact, we have also the equations 



namely the energy equation and the Lagrange- Jacobi equation (cf. 112171) ). The 
latter is precisely equation (i) in ^ 128(1 . and the energy integral 



makes any of the two equations (ii) or (iii) superfluous. 

4.5 Ceva-type relations in the spherical representation 

The classical Ceva theorem tells us that the lines from the vertices to the center 
of mass of an m-triangle 8 divide the triangle into subtriangles whose areas are 
in the proportion 



On the other hand, a point 8* on a hemisphere divides it into three spherical 
triangles with areas Ai, with the common vertex 8 and the binary collision 
points bi2, &23, &3i as the other vertices, cf. Figure 5. In this way, various 
(normalized) geometric invariants of 8 such as sides, areas, angles (cf. Figure 1) 
have their spherical "dual" counterparts, although the dual quantity may be of 
a different type. There are, for example, the three central angles a* (resp. 7Tj) 
of (5 (resp. at the shape point 8*) with J2 a i = E"'" = 27r. According to the 
following lemma, the areas Ai are dual to the angles a, , and later we also show 
the areas Aj are dual to certain sides of the spherical triangles. 

Lemma 20 Let 8 = (ai,a 2 ,a 3 ) be an m-triangle with central angle ai opposite 
to the vector a^. Then the area of the spherical triangle in M* with vertices 
bi2, 8* , &3i equals 



Proof. Let £2,-^3 be the lines spanned by the vectors a 2 and a3 respectively. 
There is an obvious piecewise linear motion with f2 = which starts at 8 and 



T-U 



= h, 



I = 2T + 2h, 



u + h-^ = \p 2 + ^ 2 + {sn^f e 2 ) 



Ai : A 2 : A 3 = mi : m 2 : m 3 



Ax = -{it -ax) 
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collapses the triangle to a degenerate configuration of shape b 3 i along L 2 , and 
continues along L 2 until the shape of bi 2 is reached. This motion keeps the 
direction of a 2 unaltered. 

On the other hand, the linear motion with Q = which collapses S to a 
configuration of shape 612 along L 3 , will rotate a 2 to a vector along L 3 which 
lies opposite to a 3 . Hence, the total change of position when 5 is deformed 
according to the above piecewise linear motion whose shape curve encloses the 
spherical triangle, is equal to the angle n— a%. Finally, by the Gauss-Bonnet 
formula, this is twice the area Ai of the triangle. ■ 

Next, we turn to the mutual distances, that is, the sides of the m-triangle S 

si = r 23 = |a 2 - a 3 | etc. 

and ask for their spherical counterpart, namely the spherical distances <r, from 
S* to the binary collision points. The quantities Sj have, indeed, a nice geometric 
interpretation in the vector algebra representation described below, see <|141[) . 
But first, by combining (|64|) and i|126|) . the identity 

9 1 — m, — L 1 — to,- L 1 — TO; . o . . 

s 2 i = ^ = —^(1 - z ) = sin 2 a z (129) 

to, TOj 1 — rrii nrii 

holds, where we have assumed normalization 1 = 1 and (as usual) X) m « = 1> 
cf. also l|14fl for notation. By summation over i the condition 1=1 reads 

1 = rhisf = V^(l — mi) sin 2 ai or to* cos 2<T; t = (130) 

where the first identity is just the normalized version of Lagrange's formula (|65|) . 
Now, let So be an m-triangle whose shape is a pole Sq = A/" or S, and let 

{ai,a 2 ,a 3 }, {(3 1 , (3 2 , f3 3 } (131) 

denote the central angles ai of Sq and the central angles = 7r^ at Sq, respec- 
tively, see Figure 5. In particular, /3 X is the longitude distance between b 3 i and 
bi 2 , or equivalently, Pi/2 is their distance in AI* . The relationship between the 
two triples in 1)131(1 follows from the above lemma, namely 

Pi = AA % = 2tt - 2a, (132) 

On the other hand, the triple of angles and the (normalized) mass distribu- 
tion uniquely determine each other. In one direction, we apply the Ceva-cosine 
law to So and obtain 



cosai = cos/?! = — -, etc. (133) 

V 1 TO 2 y 1 TO 3 (1-TO 2 )(1-TO 3 ) 

In particular, the angles are in the range 

^ < a { < 7r , < p { <n (134) 
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We also note that cos/3^ can be calculated by applying a distance function 
<7j,j 7^ i, sec (|124|) . For example, cosf3 1 — cos(2o- 3 (b3i)), and by i|126l) and the 
fact that b3i has I3 = jz^, we deduce again the above expression for cosf3 1 . 

In the other direction, we would like to know which angles in the range (|134fl 
are actually realizable for some mass distribution. The condition on the angles 
is 

2 sin /? i < ^sin/3 3 , for each i (135) 

and the corresponding (normalized) mass distribution is defined by 

2 ship, l + Ecos/3, 
mi = 1 — ==— : — 7T = — ; ^ 7, i tt , etc. 136 

£sin/^ -1 + E c °s^ -2COS/?! ' v ; 

We omit the simple proof of this, remarking that the realizability condition 
(|135() also has a nice geometric interpretation. Namely, consider the three binary 
collision points bi2,b23,t ) 3i as the vertices of a triangle in the Euclidean disk 
(of radius 1/2) with the equator circle E* as boundary. In general, let us call a 
triangle central if its circumcenter O lies in its interior. Then the realizability 
condition simply says that the above triangle must be central. For another 
property of this triangle we also refer to Lemma |2~TI 

4.6 The vector algebra representation of the kinematic 
geometry 

Since the kinematic study of m-triangles and their motions essentially involves 
spherical geometry, wc are naturally led to the vector algebra in the Euclidean 
3-space M 3 , where inner products and determinants are the basic invariants. 
Therefore, it is sometimes convenient to represent M* by the sphere S 2 (l) of 
unit vectors in R 3 and hence its cone M becomes the whole Euclidean space 3- 
space. Thus we introduce the vector algebra representation of the moduli space 
M by constructing the following transformation between Riemannian cones 

M = C{S 2 (1/2))^C(S 2 (1))=M. 3 (137) 
V :(p,r,6)^(p 2 ,2r,6) = (l,<p,0) 

which magnifies M* to a sphere of radius 1 and squares the distance to the 
origin. It is a diffeomorphism away from the base point (or origin) O, namely 
the class of the triple collision p = 0. 

In (|137|l . (I,(p,6) are the usual spherical coordinates in 3-space associated 
with the Euclidean coordinates (x, y, z), that is, 

x = /sin if cos 9, y = I simp sin 6, z = I cos if 

where < <p < tt, < 8 < 2ir. The kinematic metric (13811 on M, expressed 
as a metric on K 3 , now becomes the following conformal modification of the 
Euclidean metric, namely 

,9 ,9 9,9 dx 2 + dy 2 + dz 2 

ds 2 = dp 2 + p 2 da 2 = —7== = (138) 
4-\/:zr + y A + z z 
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A variable point on M* will be represented by a unit vector 



5* -> p = (x,y,z) e S 2 (l) cM 3 

and we fix the following notation and location of binary collision points (cf. 
lfT33|) and Figure 6) 

bas->bi = (1,0,0) 

b 3 i^b 2 = (cos/?3,sin/?3,0) (139) 
i>i2 -> t>3 = (cos/3 2 ,-sin/3 2 ,0) 



Moreover, l|12t)[) now reads 



1 



p • b, ; = cos 2a i = -^Ii = { x cos /3 3 + y sin /3 3 i = 2 (140) 

x cos (3 2 ~ V sm P2 * 

Recall that 2<7j is the spherical distance between p and b^; hence, by (|129|l 
there is the simple formula 



P b, 



(141) 



which expresses the three mutual distances Si for the normalized m-triangle 5 as 
the Euclidean distances (modulo a fixed factor) from p GS 2 to the three fixed 
points hi lying on the circle x 2 + y 2 = 1, z = 0. 

Lemma 21 There is a unique mass distribution (rn'i, m 2 , m' 3 ) such that (bi, b 2 , b 3 ) 
becomes an m-triangle with center of mass at origin, that is, 

m'lhi + m 2 b 2 + rn' 3 i>3 = 0, 

namely the dual masses m'; L = m* , cf. 

Proof. Since the triangle is central the origin lies in its interior, so there are 
barycentric coordinates qi > 0, unique up to a common multiple, such that 
qSt = 0. On the other hand, by combining (|lH()fl and l|140|) . J2 m*bi • p = 
holds for all p and consequently m| bj = 0. ■ 



4.7 An integral formula for the distance function on M* 

The kinematic Riemannian metric da 2 , expressed in terms of coordinates 1\ , / 2 
as in H83|) . may be viewed as the infinitesimal version of an integral formula 
for the distance function a(p,p') on each hemisphere of M* — S 2 (l/2). 
Our calculation of such an integral formula will be based upon a special type of 
coordinates; namely, we choose the points {b23, &3i} as a bipolar system whose 
associated distance functions {cti,ct 2 } constitute a coordinate system on each 
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hemisphere M±. In our final formula, however, we shall express the distance in 
terms of I\ , 1% , and more simply in terms of their translates F — F — m* . 

To this end, it is convenient to apply the above vector algebra representation, 
where we use the unit sphere S 2 (l) representation of M* rather than the sphere 
5 2 (l/2), and the collision points b 2 3, 631 and variable points p,p' are replaced by 
the unit vectors bi, 62, p, p', respectively. Thus, a(p,p') is half of the spherical 
distance between p and p' on the unit sphere. 

We shall reduce the calculation of the distance function a to a simple vector 
algebra involving determinants and Lagrange's formula : 



dct(bi,b 2 ,p)det(bi,b 2 ,p') = 



1 bi b 2 bi p' 
b 2 • bi 1 b 2 • p' 
p • bi p b 2 p • p' 



(142) 



By writing the left side as 



det(b 1 ,b 2 ,p) 2 det(b 1 ,b 2 ,p') 5 



1/2 



and applying Lagrange's formula to each square, the right side of (|142|1 equals 
the product 



1 bx • b 2 bi • p 
b 2 • bi 1 b 2 • p 
p • bi p ■ b 2 1 



1/2 



1 bi • b 2 bi • p' 
b 2 • bi 1 b 2 • p' 
p' • bx p' • b 2 1 



1/2 



(143) 



and this gives us an identity where bi • b 2 = cos/3 3 is a constant, the inner 
product 

p p' = cos2er(p,p') 

appears linearly, and the other non-constant entries b» • p = cos2<7i, bj • p' = 
cos 2<7^ in the determinants are of the type l|126|) (or (11401) ). 

Let D,D' be the determinants in l| 143(1 . Solving the above determinant 
identity with respect to p • p' gives 



cos2ct(p,p') = — 2 — WDD' + F) 

sm (3 3 



(144) 



where F is a bilinear form of both vectors (/i,^) and J 2 ), and moreover, 
F = and D' = sin 2 /3 3 when p' is the pole P of the hemisphere. The latter 
observation implies 

cos2er(p,P) = cos2r = 4A^/TOim 2 TO3 
Consequently, by the Ceva-Heron formula l|()9|l 



sin/3 3 



cf. dUl 



D 



sin 2 03 
mim 2 m3 



Q* 



sin 2 P3 
mim 2 TO3 



(mim 2 m 3 - Q* Q {Ii,I 2 )) 
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where Q* = Q\m* is the restriction of the quadratic form l|68f) . and 

Q* Q = (1- m 2 fll + (1 - mxfll + 2(m 3 - mi m 2 )hh (145) 

On the other hand, taking p = p' in (1 14411 gives 

F _ mim 2 m3 - Q* _ Q„ 
sin 2 /? 3 miTO 2 m 3 mim 2 m3 

and therefore, F as a function of (h, I 2 , 1{, I 2 ), is a constant times the polar- 
ization of Qq in (|14B|> . This establishes the general spherical distance formula 

<7{p,p ) = - arccosf ) (146) 

2 m\m 2 m 3 

where 

-PoZ(QS) = (1 - m 2 ) 2 hi[ + (1 - mifhl' 2 + (013 - m^ih^ + J( J 2 ) 

Remark 22 spherical trigonometry it is easy to see that the polarization 
term in jl^6p can be expressed in polar coordinates as 



sin tp sin tp' cos(8 — 9') 

5 Motions of m-triangles with conserved angu- 
lar momentum 

In the previous chapters we have investigated the kinematic quantities, their 
general relationships, and the resulting kinematic identities are valid for any 
motion of m-triangles, with no explicit assumption on the invariance of the 
angular momentum vector SI. In dynamics, however, the motion is governed 
by a potential function and one is primarily interested in the trajectories of the 
equations of motion, which is a second order ODE. In these cases the invariance 
of O is generally seen as a consequence of (rotational) symmetry properties of 
the potential function, as in the Newtonian case discussed in the introductory 
chapter. From this viewpoint, the linear motions (cf. Section 3.3) are the 
trajectories in the trivial case of a constant potential function, but still we have 
found them to be useful in our survey of the kinematic geometry of the moduli 
space M. 

In this chapter we turn to the study of virtual m-triangle motions, in full 
generality except with the explicit assumption that the angular momentum is 
conserved. Our aim is also to complete the proofs of the Main Theorems B, D, 
El, E2 and F stated in Section 2.2. 

Remark 23 It is a classical result, dating (at least) back to Weierstrass, that 
a 3-body motion with O = must be planar. There is a simple and purely 
kinematic proof of this fact. Namely, if (X(i),n(t)) is a motion of oriented 



43 



m-triangles and X = (ai,a 2 ,a 3 ) is nondegenerate, then n = if and only if 
qi = (72 = 0, where qi = n • a^. Therefore, the proof follows from the identity 

±(ai x a 2 ) x ft = (ftai + g 2 a 2 , 

where for {i, j} = {1,2}, qi is a linear combination of qi and qj with coefficients 
2Ami(l — m,j) and2Amimj, respectively. 

5.1 Moving eigenframe and intrinsic decomposition of ve- 
locities 

Consider an m-triangle motion X(i) = (ai,a 2 ,a 3 ) with a (continuous) eigen- 
frame $(t) = (ui,u 2 ,n) of its inertia tensor i?x (I24|) . For convenience, let us 
assume X(t) is nondegenerate (say, for some time interval) and hence X spans a 
plane II(X) = lin {ui, u 2 }. Any vector v in 3-space has an orthogonal splitting, 
v = v r + v r ', where v r is the tangential component lying in the plane II (X) and 
v 7 ' is the normal component. 

We will combine the splitting l|20|) of the velocity of X(t) with its tangential 
and normal decomposition, namely we decompose the individual velocities a^ 
into their tangential and normal parts and write 

X = X7 + X" = (a[, hi, kl) + (a?, ag, ag) (147) 

Recall the roles of the (instantaneous) angular velocity vector u> — u>(t) and the 
angular momentum which by the inertia operator (|23|l essentially determine 
each other. They are responsible for the purely rotational (or rigid) motion 
of the m-triangle, and in accordance with l|147|) the rotational velocity has the 
splitting 

X u =wxX = u' 1 xX + u T xX = (f) T + (f)' ] (148) 

Lemma 24 The normal velocity component of the m-triangle motion is purely 
rotational, that is, 

X" = (f)' = u T xX 
where oj t is the tangential component of uj, and moreover, 

n = u) x n = u T xn 

Proof. Consider the 2-parameter family of degenerate m-triangles 

(Aui +ptu 2 ) x X = (cin,c 2 n,c 3 n) 

where the constants Ci = Cj(A, fi) are linear combinations of A and \i with the 
relation m i c i = 0. The last identity also expresses the range of the linear 
transformation (A,/i) — * (01,02,03). 

On the other hand, for certain constants fcj 

XP = (kin,k 2 n,k 3 n), ^TOifc., = 0, 
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and consequently the system of equations c; = ki,i — 1,2,3, has a unique 
solution (A, /it), that is, there is a unique vector u/ such that X 17 = u/ X X. 
Clearly, u/ is just the tangential component uj t of a;. Finally, n is a multiple of 
ai x &2 and n is a tangential vector, so it is a simple vector algebra calculation 
to verify that n = u T xn. ■ 

It follows that the tangential velocity version of (|20|l reads 

X r =u/'xX + X' 1 

and, in particular, the horizontal velocity of an m-triangle motion is always 
tangential, whereas the rotational velocity Ijl48jl in general has a tangential and 
normal component. 

Now, we turn to the angular momentum 

n = x x r = x x (w x x) 

which we assume is a fixed vector along the z-axis, say, and the expansion 

n = flk = n T + fi" = ( ffl ui + g 2 u 2 ) + g 3 n (149) 

defines its (time dependent) coordinate vector (<?i, #2, <?3) relative to the moving 
frame #(t) . The inner product of ft with a vector v may be written as 

n ■ v = (X x u>) x X • v = (X x w) • (X x v) = B x (u>, v) 

Hence, by letting v = Vj be any of the vectors from 5, 

g. t = ft ■ Vj = B x (w, = -Vj, i = 1, 2, 3, 

and we obtain the expansion 

tt> = w + u> ' = (— ui + — u 2 ) + — n, (150) 

Al A2 1 

In particular, the rotational kinetic energy can be expressed as 



~~ 2 



x- 



5 ^( W|W ) = 5 ( g + S + ? ) (151) 



with tangential and normal parts 



Now, let us also have a closer look at the individual velocities in <|147[) and 
their splitting, namely 

a* = a[ + a? = (wi(n x a,) + ^a^ + (w T x a*), cf. d (152) 
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where un is the scalar tangential angular velocity of a^. We claim that 

Wi = + y , cf. Remark [2 (153) 

where u- 1 is the scalar angular velocity in the case of vanishing angular mo- 
mentum, and a formula is given in Theorem CI, see (|44|l . The proof of Ijl53|l is 
really the same as in Section 3.2.1, if we only consider velocity components in 
the plane II (X) and replace fl by its normal component JT'. The identities in 
Section 3.2.1, in fact, expresses kinematic relationships valid at each moment i, 
and there is no need to assume $7 is a constant. 

Finally, according to (|152|l the individual kinetic energy terms expresses as 

T« = T[ + TV = ^-ujh + Jl) + l - mi |« T x ai | 2 (154) 

and hence they depend, in fact, only on the moduli curve f(t) of the motion 
and the moving frame coordinates gi of the angular momentum. 

5.2 Final proof of the Main Theorems D, B, E1,E2 

5.2.1 The Euler equations and proof of Theorem D 

Consider the horizontal (i.e. with vanishing angular momentum) m-triangle 
motion X h (i) = (bi, b 2 , bs), with the same moduli curve T(t) as X(i) and with 
moving eigenframe 

5*(t) = (u?(t),u5(t),n(to)), 
subject to the initial conditions 

X h (t ) = X(t ), $ h (to)=d(t ) (155) 

The existence of this motion is the statement of Theorem B in the simple case 
of vanishing angular momentum (cf. Section 2.2.2). 

Let T*(t) = (<p(t), 9 {t)) be the associated shape curve on M* = S 2 , expressed 
in the usual spherical coordinates, and consider the two nearby points 

p = fo>(f ), 6(t ), p' = ( v (to + At), 9{t + At)) 

with the longitude difference AO. The meridians through p and p' intersect the 
equator circle E* in the points e and e' respectively, and there is the piecewise 
geodesic closed path 

p - c e'— > p' - p 

enclosing the shaded region as indicated in Figure 7, whose area (on the unit 
sphere) is 

AA = cos Lp AO 

modulo higher orders of AO. 
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It follows from Theorem C2 and Theorem El applied to the above path, with 
a piecewise linear motion in the plane perpendicular to n(to) ; that 

uf {t + At) ee cos(A^)ui(t ) + sin(AVOu 2 (t ) (156) 
u£(t + At) ee -sin(A^)ui(t ) + cos(A^)u 2 (f ) 
modulo higher orders of A9, where 

A^^icos^Ae* (157) 

Consequently, we infer from (|156fl and (|157|l 

"i(*o) = (^cos^)u 2 |t , u%(t ) = (-~0cosc/>)ui|t o (158) 
Next, consider the following intrinsic frame version of i|2Ull 

iii = u> x Ui + u%, i = 1,2; n = tv x n, (159) 

By taking the inner product with f2 on both sides of these identities, using (114911 
and (|15Up . we perform the following calculations : 

n ■ Ui (t ) = X Ui + iif)| to 

= f2 ( (f^ui + f^-u 2 + ^-n) x u x + (i(9coscp)u 2 
\ Ai A 2 I I 

= Sl ^-|i n +(^ + ^cos^)u 2 
^7-^)53 + 5 cos 

n ■ u 2 (t ) = n ■(<■•> x u 2 + u§)| to 

= ft f (fi Ul + |^u 2 + ^n) x u 2 - (^cos^)u! 
\ Ai A 2 1 I 

= f2 (^j-n+ (-y - ^6 I C0S( ( 9)U1 



n • n = n • (w x n) = n- ( (f^m + ^-u 2 + ^n) x n 

\ Al ^2 I 



I to 



Ai A 2 / A 2 Ai 

Since CI is a constant vector and time to is arbitrary, the above three identities 
amount precisely to the ODE i|53|) . and this completes the proof of Theorem D. 
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Finally, we turn to the precession angle \{t) which records the motion of the 
normal vector n around the z-axis, that is, the fixed O-axis. For example, using 
spherical coordinates (ip, 9) on the unit sphere in Euclidean 3-space, with (p = 
at the north pole k and \(t) — 6 if), it is a simple exercise to deduce the first 
equality in (|51Jl . and by substituting n = u> x n the second expression in 
follows by calculating cross products in the frame 

5.2.2 The lifting problem and proof of Theorem B 

To complete the proof in the general case Jl ^ 0, let us also choose a horizontal 
lifting X h (t) of f (t). Then it follows from the identity that the lifting T(t) 
is the motion X(i) ~ (a^i), a 2 (t), a 3 (i)) determined by the following initial 
value problem 

-^X = (u; x X)+^-X'\ X(* ) = r(io) (160) 
at at 

Here the vector w is a function of X and the constant vector f2, namely for 
a fixed and nondegenerate X, u> is found by inverting the inertia operator on 
3-space : 

I x :w^Xx(wxX)=fi 
The matrix of this operator is 

B = |X| 2 Id - ADA 1 ~ diag(\ 1 ,\ 2 , A 3 ) 

where Id is the identity, the vectors a^ are the columns of A, and 

D = diag{mi 1 TO2, m^) 

The eigenvalues of B are listed in Section 2.2.4, and one of them vanishes 
when X ^ is degenerate. However, in that case it is easy to check that the 
indeterminacy of a; is a summand along the line n(X) and consequently the 
summand wxXin (|16(J[1 is still well defined as a function of X and Q. 

For another proof of Theorem B, more directly related to the construction 
of a position curve 7(i) in 50(3), we use either Theorem CI or D. The first the- 
orem applies to planary motions and calculates a position curve j(t) in 5*0(2), 
recording the rotation of the vectors a^. In the non-planary case the position 
of the m-triangle is represented by the position of its moving eigenframe 3S and 
the latter is determined by the coordinate vector (gi, (72,53) of f2 relative to $ 
together with the precession angle x (of the normal vector n) in the " invariant" 
plane perpendicular to $7. The four functions <?i and x are the solution of an 
initial value problem depending only on the moduli curve T(t), as explained by 
Theorem D and formula Ij54(l . 

5.2.3 Geometric reduction of Newton's equation and proof of The- 
orem El and E2 

To derive the reduced Newton's equations from the Newton's equations Q), 
we differentiate the kinematic quantities Ii up to second order, for example, 
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I\ = 2m\B.i ■ k\ and 

I\ = 2mi |ai| 2 + 2rmai • Si 

= 47\ + 2miai • ( ^(a 2 - ai) + ^(a 3 - (161) 

V»12 »13 / 

Then we use the Ceva-cosine law (|63[) . stated in the form 

&i ■ Sij = - — — C fc = — — — (m k I k - m l I l - rrijlj), 
2mimj 2mim,j 

to replace all inner products in (|161fl by linear combinations of the I^s. This 
procedure leads to the differential equations (J5SJ. 

In the above differential equations the individual kinetic energies 2j are cru- 
cial terms, and their actual splitting l|154|) distinguishes the two cases of planary 
and non-planary 3-body motions T(t). Of course, the actual case is also decided 
by the initial data r(i ),r(io)- However, it is not decided by r(t ) and the an- 
gular momentum vector, unless r(io) is nondegenerate. Clearly, the statements 
of Theorem El and E2 must be modified if they should also cover the case where 
the initial configuration T(to) is collinear. 

First, observe that a planary three-body motion is characterized by having 
all normal kinetic energies Tf = 0, and then its moduli curve T(t) is a solution 
of the 17-reduced equations l|55(l with Tj = T? given by the first summand in 

QUI) . 

Conversely, let ^(t) be a moduli curve which is a solution of this ODE, 
for a given value of Q. By Theorem B there is a unique lifting X(f), namely a 
virtual motion in the xy-planc, with angular momentum f2k and specified initial 
position X(£o)- From this knowledge we may calculate the initial velocity 

X(t ) = w(t«>) x X(to) + ± h (t ), w(*o) = y^yk, 

since the horizontal velocity ~X tl (to) is determined by gs7(io)- On the other 
hand, Newton's equations (JTJ also has a unique solution T{t) with the above 
initial conditions, namely T(to) = X(i ) and T(t ) = X(i ), and clearly the 
moduli curve of T(t) is a solution of the fi-reduced ODE. By uniqueness of the 
lifting we conclude that T(t) = X(i) for all f, and this completes the proof of 
Theorem El. 

Next, we turn to the general case, described by Theorem E2. The kinetic 
energies Tj in the fi-reduced ODE have the general form (|154|> . and we claim 
they depend only on the moduli curve and the functions gi. This clearly holds 
for the tangential summand TJ , whose expression is even independent of g\ and 
52- On the other hand, the normal summand is, say, for i = 1: 

2 



1 2 1 

7i = -mi |w T X ai| = -mi 



(-pUi + f^-u 2 ) x (cos^iUi + sini/ , iU 2 ) 
Ai A 2 



^mi(fi- sin^i - ^ cos^i) 2 (162) 
^ Ai A2 
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where ipi is the angle between Ui and ai, satisfying 

cos w 9 __ 

tanVi = —tan-, cf. Theorem d 163) 

1 + sin <f 2 

In this formula 9 is the longitude angle measured from the binary collision point 
bi and is increasing in the direction towards t>2. It follows (e.g. by symmetry) 
that one obtains the corresponding formula for T% and from (|162f) when 
ipi is replaced by the corresponding angle ip 2 an d determined by the same 
formula (|163Jl with 9 measured from 62 or S3, respectively. This proves the 
above claim. 

Problem 25 It is an interesting task to simplify the expression for the energies 
, and to express them in terms of coordinates Ij as in the case of T[ . But 
we leave the topic here. 

Thus, in the general case our "reduced" ODE actually consists of the reduced 
Newton's equations l|55[) together with the Euler equations (|53|l . The initial 
data will be a given nondegenerate configuration X(io) and angular momentum 
vector r2, and as before, this determines the initial velocity X(io) and hence the 
motion X(t) is unique (by Newton's equation Q). 

However, to avoid ambiguity in the initial value problem for the "reduced" 
ODE we must also specify the initial orientation (i.e. the normal vector n(to)) 
which decides whether the shape curve starts out on the upper or lower hemi- 
sphere of M*. With this choice the initial eigenframe {ui, u 2 , n} \t=t > an d 
hence also the initial values gi(fo), will be unique relative to a fixed convention, 
say, the angle ip l between Ui and ai is (initially) in the range (|114|l . Now, the 
remaining part of the proof of Theorem E2 is similar to the proof of Theorem 
El. 

5.3 Geometric reduction of the least action principles 

Recall from Section 1.2 the two classical least action principles with the action 
integral J\ and J2, respectively. The underlying geometric structure naturally 
associated to the former is Riemannian while that of the latter is symplectic. 
Thus the two types of least action principles are radically different in their 
basic geometric setting, although both of them characterize the same motion. 
Indeed, it is easy to verify that the Euler-Lagrange equations of both variational 
principles are equivalent to the Newton's equation of motion (|T|). 

However, the classical approach to the three-body problem is mainly based 
on J2 , namely the least action principle of Hamilton and the Hamilton- Jacobi 
theory, in the framework of canonical transformations and symplectic geom- 
etry. On the other hand, the kinematic geometry of m-triangles is, on the 
other hand, more naturally associated with the least action principle of Euler- 
Lagrange- Jacobi and the action integral J\ , and it is in this geometric framework 
that we have established many basic results, such as the universal sphericality, 
the kinematic Gauss-Bonnet formula (and geometric phase), kinematic moving 
frames and the generalized Euler equations. 
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5.3.1 Proof of Theorem F 



We consider virtual 3-body motions in the xy-plane, represented by motions 
S(t) of oriented m-trianglcs with k as their common normal vector, and hence 
a (nondegenerate) m-triangle S — (ai,a2,a 3 ) is positively oriented if a!xa 2 
points in the direction of k. The corresponding configuration space is R 4 , see 
and R 4 /SO(2) = M is the full moduli space. The SO(2)-orbit 5 of an 
oriented m-triangle S will be regarded both as a point in M and as a subset 
(congruence class) of R 4 . In the sequel, all motions in R 4 are also assumed to 
have a constant angular momentum. 

The set of differentiable (C^-smooth) curves in M from 6 to 6% is denoted 
tygg 1 , and similarly ^ s j t denotes the set of differentiable curves in R 4 which 
start at S and terminate at the orbit Si. For a fixed angular momentum Q the 
kinetic energy 

n 2 

T = T n =T + T u > =T+— , 

the potential function U, the Lagrange function Lq =Tq + U and total energy 
h = Tq — U, are functions which are also defined at the level of M. Therefore, 
for fixed value of f2, h or time interval [0, ti], we may consider corresponding 
subsets of g 

VjaW. Via (MJ) (164) 

with the obvious meaning, and similarly subsets of V[$g g 

fe(^), Vefrfah), qJ, tJl (n,[to,ti]) (165) 

Clearly, the solution curves of Newton's equation belong to sets of type (|165|) . 
We can also define the reduced action integrals 



Ji,a = j T n dt , J 2 , n = / L Q dt (166) 
acting on moduli curves, and there is the following commutative diagram 



An \ / J it a (167) 

K 

where the map ttq takes a curve T(t) to its moduli curve T(t) and J^q is the 
restriction of Jj to the subspace of curves with fixed angular momentum f2. 
Moreover, for S a nondegenerate m-triangle the map ttq is, in fact, a bijection 
due to the unique lifting property described by Theorem B. 

Now, let us turn to the proof of Theorem F, which we restate as follows : 

Theorem 26 The solution curves of the planary VL-reduced Newton's equation 
can be characterized as the extremal curves of Jx,n (resp. J2.f1) restricted to the 
sets 

yg- Sl (h), resp. fetMiD 
of moduli curves, with fixed energy h or time interval [to,tx], respectively. 
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Proof. First of all, extremal curves of Ji,q (resp. J2,n) are solutions of the 
associated Euler-Lagrange equations, and one checks that these are second order 
ODE whose solution curves are (as usual) uniquely determined by their initial 
position and velocity. 

On the other hand, the fi-reduced Newton's equation is also a second order 
ODE whose solution curves are uniquely determined by their initial position 
and velocity. Therefore, to show that the Euler-Lagrange equations and the 
^-reduced Newton's equation have the same solutions it suffices to verify this 
locally. More precisely, it suffices to show that any small segment of a solution 
curve of the Jl-reduced Newton's equation is also an extremal curve of J^n 
(resp. J2,n)- 

Let r be a small segment from S to S\ of a solution curve of the fi-reduced 
Newton's equation. We may assume the points are sufficiently close to ensure 
that r is the only segment (of a solution) linking them. 

Choose 6 in the orbit 8. Since the actions in H166(l are always nonnegative 
and the orbits 6 and Si are sufficiently close, there exists a Ji -minimizing (resp. 
^-minimizing) curve T in R 4 between S and the orbit Si, say Si is its end point. 
Then V is a solution of Newton's equation and it must, in fact, be the lifting of 
r. Moreover, it is the unique curve in tyj s g (Q) with minimal action integral of 
J i (resp. J 2 ). 

Consequently, T is, of course, also a small segment of an extremal curve of 
J^o and therefore by 1)1671) its moduli curve T is a small segment of an extremal 
curve of J^n- ■ 

Remark 27 Theorem F does not extend to the case of general three-body mo- 
tions. The reason is that the kinetic energy of a motion S(t) does not depend 
only on the moduli curve and the angular momentum vector S7, but also on the 
instantaneous configuration S(t). Hence, one cannot proceed as above, since it 
is not clear what should be the appropriate action J^o. at the moduli space level. 

6 The Newtonian potential function 

In this chapter our primary task is to analyze the Newtonian potential func- 
tion (J2J and its crucial dependence on the mass distribution. The function is 
naturally defined at moduli space level, 



where in each half-space M± — Kj_ the two triples (ii, l2,h) an d (si, 82,83), 
related by (|6~4*j) . are natural coordinate systems. Certainly, U has the simplest 
possible form when expressed by the mutual distances Sj. Even so, sometimes 
it is also convenient to use Euclidean coordinates or their associated spherical 
coordinates (I,ip,6), where I = p 2 = \/ x 1 + y 2 + z 2 , as explained in Section 




3 



- 3/2 



(168) 



4.5. 
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Let U* be the restriction of U to the "unit" sphere M* = S 2 = (I = 1). 
As a function of the coordinates Ij (resp. Sj) U is homogeneous of degree — \ 
(resp. —1), namely 

U(h,I 2 ,I 3 ) = lu*fe,%,h = -U*(S*) = -U*(^6) 
p I I I p p 

where 5* < — ► ((f, 9) represents the shape of an m-triangle. For the sake of 
convenience, the formula for U* in terms of spherical coordinates is 

3 • 3/2/ *\ — 1/2 

TT*=X2 , ^ K) (169) 
~[ sin <^ cos (0 - 6i) 

where 9\, 62, #3 are the longitude angles of the binary collision points 623, &31, &12, 
respectively. This follows by substituting the expression (I123[l with 6\ = 6 — 9i 
into In particular, with the convention that Q\ = made in Remark H5l 

we must use 

(01,03,03) = (0,/3 3 , 

where the angles f3 i; described in - (jlHrSfl . measure the longitude differences 
between the points bki- 

The analysis of U trivially reduces to that of the restriction U* = U |m* • In 
fact, by symmetry it suffices to investigate U* on the closed upper hemisphere, 
0<(y5<7r/2. U* has no maximum points since U* tends to 00 at the singular 
points bki (which are poles of U*). On the other hand, it is a classical result, 
dating (at least) back to Lagrange, that U* has a unique minimum value at the 
shape of a regular triangle, namely (for p = 1) 

si = s 2 = s 3 = — =, or li = l-mi — (170) 

Vm m 

This defines a unique point on each hemisphere M± which we also refer 
to as the physical center (as opposed to the poles which are the geometric 
center). It is easy to prove the above statement using the coordinates Sj and 
Lagrange's multiplier method subject to the constraint 1 = 1 = ^to*s|, cf. 
113UI) . Another proof follows from Lemma QUI below. The spherical coordinates 
of the shape (|170|l is worked out in Section 8.8, cf. 



6.1 Vector algebra analysis of the Newtonian function 

The vector algebra representation of M in Section 4.5 is also a convenient setting 
for the local analysis of U* , namely the function on the unit sphere of Euclidean 
3-space defined by 



i=l i=l IP ' 
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where 



P t» 



Si 



(cf. (mi) 



is the Euclidean distance from p to the binary collision point 6, . 

Fix a point p on the sphere with z > 0, and consider nearby points p' 
p + x on the sphere, that is, x is subject to the constraint 



2p • x + x • x = 0, 



(172) 



which in turn implies 



p'-b, 

Hence, there is the following expansion at p 

3 TI*(n) 3 oo 00 

U*( P + x) = J2 n _ u/2 = £ u:{p) £ c " z? = £ F " (p; x) (173) 

i=l V Z ^ -— 1 ™ = ° n 

where we use the notation 



P-t>i 



- 2b, • x 



2b, : • : 



,m*TOi. b, ; • X 



P bi 

1 • 3 • 5 • • • (2n- 1) 
2™n! ; 



c = 1 



Remark 28 It is easy to see that the convergence condition \zi\ < 1 for the 
series o/£/*(p + x) in jl73\f is equivalent to the condition 

( P + x) • b 4 > 2 P ■ b 4 - 1 

Geometrically, this means that p + x belongs to the hemispherical cap Di cen- 
tered at hi which is cut out by the plane parallel to the tangent plane at hi and 
separated by the distance Ri — 2(1 — p • bj). In particular, Di is the whole 
hemisphere ifp ■ h{ < 0, and by Lemma WR we know p • bi < holds for at least 
one i. The domain of convergence for the series ofU* is the "polygonal" region 

The zero order term of the expansion l|173|) is, of course, Fq = t/*(p), and 
the first order term is 



*i(p;x) = ( 



fcibi 



k 2 h 2 



^b 3 





3 




3 1 




P bi 




P b 2 




P b 3 



(174) 



which is essentially the gradient of U* at p. To make this precise, consider the 
following function from 5 2 (1) to the xy- plane 



p^B(p)=£. 



(175) 



»=i p - bj 
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Lemma 29 The gradient vector of U* at p is given by 

V*7*(p) = B(p)-(B(p)-p)p (176) 

Proof. Since by (|174|l 

Vf7*(p)-x = F 1 (p;x)=B(p)-x 

and x is tangential to p in the limit as x — > 0, it follows that 

VCT(p) t = B(p) t 

holds for all tangent vectors t. Hence, the tangent vector V?7*(p) is the orthog- 
onal projection of B(p) in the direction of p. ■ 

Lemma 30 The zero points of B are the two critical points of U* outside the 
equator circle z = 0, namely the physical center p^ (cf. on each hemi- 

sphere z > or z < 0. 

Proof. By Lemma |2~T1 B(p) = if and only if for some constant A > 0, 

3 

= A, 





3 




3 




P - bi 




P b 2 




P b 3 




— = m.-,- 




— = mt- 





fcl 

and by l|141|) . this is equivalent to sf = = s| = -|, namely si = S2 = S3. In 
particular, p is a point with z 7^ 0, cf. (|177c|l . On the other hand, for z ^ it 
is easy to see that B(p) = if and only if VJ7*(p) = 0. ■ 

The identity (|176fl also implies that the critical points of U* on the unit circle 
z = are the "eigenvectors" of B, in the sense that B(p) = Ap for some A, 
necessarily equal to B(p) • p ^ 0. Clearly, B has a pole at bi and B(p)/ |B(p)| 
tends to b^ as p tends to bY A simple analysis of B will show there are exactly 
one solution p between each pair hj, bj, of poles. These are the so-called Eulcr 
points = 1,2,3, and they are the saddle points of U*on the 2-sphere. We 
omit the proof of this well known fact. 

6.1.1 Series expansion of U* at its minimum point 

Henceforth, we shall focus attention on the expansion l|173|) at the physical 
center po = (x,y, z),z > 0, that is, po is the minimum point of U* on the 
hemisphere z > 0. The coordinates of po are the following mass dependent 
constants 

x = Po • bi = 1 - ^r—; — (177a) 

•mra\m\ 



-^( m2 ~ m3 ) (177b) 
m mi + HI3 



cos2ro = 4A -v/ TO i m 2'Ti3 = Vi^-r- (177c) 

771 
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where Ao is the area of the regular triangle with 1 = 1. These expressions follow 
from 1111(1 . (I140|) and 11701) . Note, for example, that z becomes arbitrarily small 
when some mass rrii tends to zero, and 2 = 1, that is, po is the north pole Af, 
precisely when the masses are equal. Concerning the sign of y = ±\/l — x 2 — 2 2 , 
we used l(14()|l to check, for example, that y > if m 2 > m 3 . 
The constant term of the series <|173[) is the minimum value 

F = U*(p Q ) = ]T U*(p ) = = ™ V2 > ( 178 ) 

and Fi = 0, of course. Moreover, by (|170[) . rhs 2 = 1 holds for all i, and therefore 
the n-th order term can be written as 

F n (p ;x) = — (m im *) n (b % ■ x)" (179) 

V m ) i= i mi 

Now, let us turn to the local analysis of the series, namely U* developed 
as a power series in suitable coordinates around po- To this end, consider a 
positively oriented orthonormal frame of the Euclidean 3-space of type 

{ti,t 2 ,p }, tiGlI fori = 1,2 (180) 

where II is the tangent space of the sphere at po. By condition l(172|l . the 
components of x = p — p must satisfy 

x = £ti+?7t 2 +Cpo, 2( + e 2 +^+C 2 =0, (181) 

and consequently the map 

where 

c = -l + \/i-(e+v 2 ) = -\{e + v 2 ) \(e + v 2 ? (m 

is a parametrization of the region of the sphere lying above the plane through 
the origin and parallel to II. Geometrically, the unit disk of II is projected down 
to the sphere in the direction of the normal vector p - 

It is natural to choose ti tangent to the meridian through po. Indeed, the 
intrinsic nature of this condition will lead to a series expansion whose coefficients 
are essentially symmetric functions of the masses rn^. The projection in the xy- 
plane of such an orthonormal basis {ti,t2} of II is (up to sign) given by 

vl — z VI — z 

In order to express the n-th term l|179ll of the ?7*-series in terms of £, T], we 
proceed as follows. Expand the "variables" in 11791) 

b i -x = a' i Z + b , i r 1 + c' i (, i= 1,2,3 (184) 
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where the three triples (a' i ,b' i ,c' i ) of coefficients are specific functions of the 
parameters rrii determined by 

a[ = b t ■ ti = b, • ti, 6^ = bi • t 2 , c- = b,-p (185) 

For example, by i fE&jl and (tT%3|) . 

/ xz , ~y , 

" ~ > h = n -o > c i " " 



z 



where x, y, z are the known functions in (|177a|l - (|177c|l . and the three triples in 
(|185|) permute covariantly with the parameters rrii. 

In view of 1)1791) . it is slightly more convenient to replace (|184H by 

(mim*)bi ■ x = + b t r\ + qC, i = 1,2,3 

and hence we calculate the modified coefficients, namely 

V3m^ , * m 



y/rh 2 - 3m m 



= (mi + i-m, +2 ) (i mod3) (186) 



2 v / m 2 - 3 

* m 
c, = m^, - 

m 

Thus, with the coordinate system (£,77), we finally arrive at the following sym- 
metrization of the power series in 1)1731 ) : 

OO 3 00 ✓ V 

[/*(£, 77) = E*»E — + m + c '0" = E E ( " ) ^.**V 

(187) 



n— 2—1 n— j+fc=n 



where 

1 • 3 • 5 • • • (2n - 1) m n +i 



Kn = 



2 n n\ m n - 1 

and Aj t k, k even, are symmetric functions of the masses rrii, whereas Aj t k is 
alternating symmetric when k is odd. 

Remark 31 The Newton sums Sk = ^2 m i are ' °f course, polynomials of rh 
and rh, cf. \15)) . For example, 

Si = 1 - 4m + 4m + 2m 2 

In terms of the Sk it is rather straightforward to obtain explicit expressions for 
the above symmetric functions A^k- For k odd, the alternating function Aj^ is 
a product of a symmetric function and the basic alternating function 

21 = (mi - m 2 )(m 2 - m 3 )(m 3 - mi) (188) 
= ^2 m,:(l - mi)(m i+1 - m i+2 ) 

i mod 3 
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6.1.2 The quadratic term of U* 

We shall work out explicitly the quadratic term of the function U* expanded at 
the physical center p , namely the n = 2 term of l|173[) or H187|) 

F 2 = ^^(Ae + 2B0l + Cr?) = «(Aif + A 2 ^ 2 ) (189) 
8m 

where (£, fj) is the coordinate system of a diagonalizing frame {ti,t 2 } of the 
tangent plane II at po- 

First, let us determine the coefficients A = A 2 .o, B = Ai^,C = A , 2 in (|187fl 
as symmetric or alternating functions of the symbols mi. By using l|188|l . the 
identities 

mj(l — rrii) 2 = m + 3m, rrii(mi+i — rrii + 2) 2 = m — 9m 
and the expressions (|186|l . we find 

\ - 1 9 1 3m ,m 2 + 3mm — 4m. 

A =y — a i = 7— ( — —o — tt- — ) 

z — ' m ? ' 4 m m z — 



„ \ - 1 , 9 1 m(m — 9m) 
G =) fe i = 7" 



m; 1 4 m 2 — 3m 



Hence, the eigenvalues of the quadratic in (|189fl are determined from the equa- 
tions 



which yield 



A 1+ A 2 = A + C=?, \ 1 \ 2 =AC-B 2 = ^, 
m 4 m 



nrjy . 

A* = — (l±Vl-3m) (190) 
Lm 



Consequently, 

f 2 = Am 3 / 2 ((l±, 1 )f + /i = VT^3^, (191) 

where (£, 77) are coordinates with respect to the eigenvectors 

ti = cos a ti + sin a t 2 , t 2 = — sin a ti + cos a t 2 

obtained by rotating the intrinsic frame {ti,t 2 } of the plane Ho. 
Similar to 1|117|) . the angle a is given by 

2B -\/3m(mi - m 2 )(m 2 - m 3 )(m 3 - mi) 

tan2a=- = ; — (192) 

A — C m + 9mm — om 

Then it also follows from (|119|l that the largest eigenvalue Ai ~ l+/i corresponds 
to the vector ti if and only if sin 2a and B have the same sign. 
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Remark 32 In the simplest case of uniform mass distribution, mi = 1/3, the 
analysis of U* is much simpler than in the general case. In this case, where the 
geometrical and physical center coincide, some of the above expressions such as 
^19iH) . are of indeterminate type. However, if mi — nij m^ then tan 2a = 0. 
and hence the frame {ti,t2} is already diagonalizing. 

7 A geometric setting for the study of triple col- 
lisions 

Recall the well known fact, stated by Weierstrass and proved by Sundman (cf. 
|14|. |15p . that three-body motions leading to triple collision must have vanish- 
ing angular momentum, = 0, and consequently they are also planary. Thus 
we shall focus attention on planary virtual motions X(t), namely curves in the 
configuration space M ~ R , with zero angular momentum, and we continue 
to use the vector algebra representation (cf. Section 4.5) of the moduli space 
M = Mq/SO(2) = R 3 , where the (equator) xy-plane E = M 2 represents con- 
gruence classes of eclipse (i.e. collinear) configurations. 

The total kinetic energy T — T can be expressed as a positive definite 
quadratic differential form on M \ {O}, namely the kinematic Riemannian 
metric 

o 2 

ds 2 = 2Tdt 2 = dp 2 + p 2 da 2 = dp 2 + tjds 2 , cf. JIHJ (193) 

which describes M as the Riemannian cone over the shape space M* — S 2 (l/2), 
and 

ds 2 = Ada 2 = dip 2 + (sin 2 ip)d9 2 (194) 

is the metric of the magnified sphere S 2 (l). 

Following Jacobi, we introduce the following conformal modification of the 
metric (|193|) for each energy level h, namely 

ds 2 h = (U + h)ds 2 , (195) 

which we refer to as the physical metric, and transform the action integral 
of 1|166|) into the arc-length integral in the Riemannian space 

(M h ,ds 2 ): M h = {peM;U(p) + h>0} (196) 

Consequently, the trajectories of three-body motions with total energy h are 
mapped to curves in M which are geodesies in the space (Mh, ds\). Notice that 
reflection in the equator plane E C M, that is, the transformation Lp — * n — tp, 
restricts to an involutive isometry of the Riemannian space 1196|) with the eclipse 
subspace E h — E C\M h as fixed point set, and hence this is a totally geodesic 
submanifold of Mh- In particular, a geodesic curve in (Mhjds^) is transversal 
to Eh, unless it lies entirely in Eh- Moreover, for a simple geometrical reason, 
a shortest geodesic in (M^jdsf) between a point outside Eh and the origin O 
cannot have any intermediate intersection with Eh- 
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Finally, we note that Newton's equation (JIJ has a 1-parameter group of 
space-time scaling symmetries {g s }, where g s sends a solution X(t) to a solution 

Y(t) = e 2s / 3 X(e- s t) (197) 

and changes the energy level from h to e~ 2s ^ 3 h. Hence, all the Riemannian 
structures in i|19tjfl with energy h of the same sign are mutually homothetic, 
and consequently there are essentially only three distinct cases, namely when 
the total energy h is negative, zero or positive. 



7.1 Geodesic rays and distance estimates 

Clearly, for h > the variety Mh is the whole moduli space M, whereas for 
h < it is the star-shaped union of all ray segments 



p e M* (198) 



from the origin O to the point where the ray through p intersects the boundary 
dMh, that is, the level surface U = —h. By definition, the physical metric (|195ll 
vanishes on the boundary, meaning that the distance between any two boundary 
points is zero. 

For h < the length of any ray segment (|198fl is 



L h (p) = I ds h = J J + h dp, (199) 
o ' 

and therefore there is a unique pair of shortest length 



L h = L h (pt) - j^ + hdp (200) 
o 

where the points p^ on the hemispheres z > and z < represent the shape 
of a regular triangle and hence U* has the minimal value 

^ = t/*(pj) = m 3 / 2 , cf. (UZHJ 

The following is a useful fact in Riemannian geometry which follows from 
general analysis of the first variation of arc-length. 

Lemma 33 Let ds 2 and ds 2 be two Riemannian metrics on a given manifold 
such that 

ds 2 = f 2 ds 2 

where f is a smooth and positive function, that is, ds 2 is a conformal modifica- 
tion of ds 2 . Let r be a C 2 -smooth curve and let n denote a normal vector at a 
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given point on T. Then the geodesic curvatures ofT in (the direction o/n) with 
respect to the two metrics are related by 

fc(n)=/C(n)-^ln/ 

We shall apply the lemma to the kinematic and physical metric, namely the 
metrics ds 2 and ds\ on M , cf. Ijl95(l . Thus, a moduli curve T is a geodesic with 
respect to ds^ if and only if its geodesic curvature AC(n) with respect to ds^ 
vanishes for all n, or equivalently 

K(n) = l-4-HU + h) (201) 
2 an 

where /C(n) is the geodesic curvature in the normal direction n, with respect to 
ds 2 . 

The simplest type of 3-body motions are the shape invariant ones, that is, 
the shape curve is a single point on M* = S 2 and hence the moduli curve 
is confined to a ray emanating from O in the cone M. Rays are, of course, 
geodesies with respect to the kinematic metric ds 2 , consequently a ray through 
p 6 S 2 (or a ray segment <|198|) if ft. < 0) is also a geodesic of the metric ds\ if 
and only if the normal derivative vanishes in all directions n normal to the ray, 
that is, 

■^-ln(-U* + h)=0 
an p 

This condition is independent of the radial coordinate p, and for p = 1 the 
vectors n span the tangent plane of S 2 at the point p. Consequently, the 
solutions are the five critical points p of U* , namely the three saddle points 
(called Euler points) on the equator circle p = 1 in the xy-plane, and the pair 
p of minimumspoints (also called Lagrange points). Thus, there are altogether 
exactly five geodesic rays (or ray segments) in (M/j,ds^). 



are the unique shortest 



Lemma 34 For h < 0, the two ray segments 

geodesic curves in (Mhjds^) linking a boundary point and the base point O 
(ignoring curve pieces of zero length along dMh). In particular, the distance 
from dMh to O is the number Lh in 

Proof. In (Af, ds 2 ), let Bh be the geodesic ball of radius jj^ centered at O. 
It lies inside Mh and touches dMh at the two points n^Po • If f is any curve 
between O and a point q on dM h , let f i be the portion of f between O and 
the first point qi on dMh- From the calculation 



L(T)>L(T 1 )= I Vh + Uds> I lh+^-ds> / * h+^-dp 



it is clear that the two ray segments of length Lh are, indeed, the shortest curves, 
and they are unique (modulo a portion along the boundary). ■ 
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Remark 35 For h > 0, the geodesic rays through p are still length minimiz- 
ing, whereas for any h this fails for the three geodesic rays ( or segments ) in the 
eclipse plane E. 

7.2 Existence of triple collision motions with minimal ac- 
tion 

We shall combine the above differential geometric setting, Theorem F and 
Hilbert's direct method to study the existence problem of three-body motions 
leading to triple collision, starting from a given non-degenerate m-triangle and 
with minimal action integral J7|), say. When this problem is pushed down to 
the level of M, at a given energy level h, it can be reduced to the problem of 
existence of a shortest geodesic, with respect to the metric ds h , between a given 
point in Mh — E and the base point O. 

Remark 36 For any energy h there are Newtonian motions, with the constant 
shape of a regular triangle or an Euler configuration, through which the m- 
triangle shrinks homothetically to a triple collision in finite time. To find the 
time parametrization of such a three-body motion, in fact, amounts to solve a 
two-body (or Kepler) problem, and this leads to the classical solutions found 
by Lagrange and Euler, see J^, J^, \VJ$ . For these motions minimal action 
is achieved for the Lagrange motions (regular triangle), but not for the Euler 
motions (which are collinear). 

More generally, let us first consider the case h < 0, and define the variety 

D h = {p£ M h ; d(p, O) < d(p, 0M h ) + L h ) (202) 

where d(p, O) (resp. d(p, dMh)) is the distance between O and p (resp. dMh) 
in Mh with the metric ds\. The interior Dh (resp. boundary dDh) of Dh is 
defined by strict inequality (resp. equality) in (|202H . 

Remark 37 The two surfaces dMh and dDh are interesting geometric objects 
in the study of triple collision orbits. They touch each other at the two points 
j^jP^ on dMh closest to the cone vertex O. 

We will prove the following existence result : 

Theorem 38 (cf. [S], Theorem 5) Starting from a given oriented m-triangle 
S whose congruence class S belongs to Dh, h < 0, there is a three-body motion 
with total energy h and minimal action integral which leads to a triple collision. 

Proof. As a consequence of Theorem F, the proof reduces to the existence of 
a curve with minimal length in Mh linking 5 to O. For 5 in Dh such a curve is 
necessarily a geodesic. 

Assume first 5 € Dh, and let {Pi} be a sequence of curves in Mh between 8 
and O whose lengths satisfy 

L(f i) < d{5, dM h ) + L h , lim L(f % ) = d(~S, O) 
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In particular, each f j is disjoint from the boundary dDh- 

Let us divide f j into ml 1 segments of equal length and replace each segment 
by the unique shortest geodesic between its end points. Then it is quite straight- 
forward to apply the direct method of Hilbert to find a suitable subsequence 
of {fi} with a limiting curve f , and this is necessarily a geodesic curve in Mh 
between 8 and O with minimal length L(Y) — d(S, O). 

Assume next 8 £ dDh, and let {8k} be as sequence of points in Dh with 8 
as its limit. Moreover, let be a sequence of curves, where f & is a shortest 

geodesic between 8k and O, that is, L(tk) = d(8k,0). It follows that 

lim d(8 k ,0) = d(S,0) (203) 

k — >oo 

and it is not difficult so see that there is a suitable subsequence of {Tfc} with a 
limiting curve Y whose length is the limit (|203|l . and moreover, T is a geodesic 
between 8 and O. ■ 

Finally, we consider the case h > 0, namely when Mh = M, and then we 
have the following analogue of the above theorem . 

Theorem 39 (cf. '31, Theorem 5') Starting from a given oriented m-triangle 
8, for a given energy level h > there is always a three-body motion leading to 
triple collision and with minimal action integral. 

Proof. This is similar to the case 8 6 Dh of the previous proof, and the 
application of Hilbert 's direct method will give the existence of the curve we 
seek. ■ 

Remark 40 The direct method of Hilbert can, of course, also be applied to 
study the existence problem of a geodesic curve Y realizing the minimal dis- 
tance between two given points 8\ and 82 of {Mh,ds\). Then, by Theorem B, 
there are liftings YofY which are planary three-body motions (with specified an- 
gular momentum) starting from a given oriented m-triangle 81 belonging to the 
congruence class 8\ . However, the end point configuration 82 ofY is already de- 
termined by Y and 81, according to Theorem C2. Consequently, only three-body 
motions with specific relative positions of their initial and terminal m-triangles 
Si can have minimal action integrals. In the above two theorems there is no 
such relative position constraint since the triple collision configuration 82 = O 
consists of a single congruence class. 

7.3 The uniqueness problem for triple collision motions 
with minimal action 

In view of the above existence theorems, it is natural to investigate the follow- 
ing uniqueness problem for motions starting from a given configuration at a 
sufficiently large energy level. 

Problem 41 To a given non-degenerate oriented m-triangle 8 and energy level 
h above a lower bound, say h > t)(S), is there a unique three-body motion from 
8 leading to a triple collision with minimal action integral ? 
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This problem requires a considerable amount of in-depth analysis of the 
geodesic equation of (Mh, ds 2 ), and here we shall leave it as an open problem. 

However, to facilitate future analytical studies of the above problem and 
related problems we shall discuss a geometric reduction technique which reflects 
some useful feature of the Riemannian structure of M = C(M*) as a cone over 
the subspace M* = (p = 1). In M the integral curves of the vector field -g- are 
the rays emanating from O, and they define the radial (i.e. a natural "vertical") 
direction at every point ^ O. 

The above problem is, indeed, simple and has an optimal solution in the 
special case mentioned in Remark namely for the shape invariant motions 
of Lagrange type. A 3-body motion is shape invariant if its moduli curve T is 
confined to a ray, that is, the associated shape curve T* is a single point. In 
fact, if a Newtonian motion is shape invariant over some time interval of length 
> 0, then for all time T* is a single point (necessarily a critical point of U*). 
However, in view of Remark 1361 the collinear solutions with the shape of an 
Euler point are not even action minimizing. 

In general, let T be a smooth curve in M — {O}. Since the associated shape 
curve F* is the radial projection of T onto the transversal subspace M*, T* will 
be smooth as long as f is transversal to the radial direction, whereas a cusp 
may occur at points where this fails. Hence, in the long run the typical shape 
curves of 3-body motions are rather piecewise smooth, but still they can be 
parametrized by arc-length. Moreover, unless T* is a single point, T may also 
be parametrized by the arc- length parameter of T* . 

Definition 42 Let T be a curve in the moduli space M and let T* be the as- 
sociated shape curve. The cone consisting of all rays emanating from O and 
passing through points on T (orT*) is called the cone surface ofT (or V*), and 
it is denoted either C(T) or C( T*). 

We assume T (and hence also T*) has a given orientation. Since the metric 
on M* = S 2 (l/2) is denoted by da 2 (cf. e.g. I38fl . a also denotes the arc-length 
parameter of T*. For a ranging over some interval [<7o,0'i], the corresponding 
surface C(T*) is immersed in (M,ds 2 ) with the induced kinematic metric 



and hence it is isometric to a flat Euclidean sector of angular width o\ — o~o, 
with (p,cr) as polar coordinates centered at the origin O. 

The moduli space M ss R 3 has the standard (right handed) orientation and, 
in particular, the 2-sphere M* has the induced orientation with as positive 
normal vector field. The surface C(T*) is naturally oriented with the positive 
orthonormal frame 



and we choose its normal vector field u so that v followed by the frame (|205|l 
is a positive orthonormal frame in M. 



^ s2 |c(r*) = dp 2 + p 2 da 2 ; a < $ < ai, 



(204) 




(205) 
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In H2U4|) the curve T* becomes the circular arc of radius p = 1, whereas the 
(original) moduli curve T appears as a radial deformation of T* . When s in (|204l 
is viewed as the arc-length parameter of 

f : s -> (p(s),a(s)), 

l|2U4|) becomes an identity along the curve. 

The extrinsic geometry of C(T*) c M is completely determined by the 
extrinsic geometry of T* C M*. Indeed, the lines of curvature are the two 
families of coordinate curves, namely the rays (a constant) and the "circles" (p 
constant) in C(T*). The principal curvature of C(T*) at a point 6 is zero in 
the ray direction and is equal to JC*/p in the direction of J^, where K* is the 
geodesic curvature of T* in M* at the corresponding point 8*. 

Along the curve T we will also consider the positive orthonormal moving 
frame {T,rj,i/}, where r is the unit tangent vector in the (chosen) positive 
direction of f , and hence the frame {r, rf\ of C(T*) differs from the station- 
ary frame l|205l) by a certain rotation angle a. Namely, we define the (radial) 
inclination angle a of V by wrtiting 

d . 1 d .8 19 ,„„„- 

T = cosa— — h since— — — , r) = — since— — hcosa-— — (206) 
op p oa Op p Oa 

dp da 1 dp d 

cos a — —, sin a = p— , cot a = — — = — — In p (207) 
ds ds p da da 

Briefly, a is the angle between the ray direction and the tangent direction, and 
< a < 7r since sin a in H207|l is not negative. The extreme values a = 0, n 
occur when T is not transversal to the radial direction, in which case a and T* 
(as functions of a or time) may encounter a singularity, namely T* encounters 
a cusp. 

Remark 43 The angle a{a) and the radial distance p{a) are mutually depen- 
dent according to \201\j . For example, we have for p(cro) 7^ 

p(a) = /9(er )exp( / cot a(a)da) (208) 

J (T 

The geodesic condition for a curve T in {Mh, dsfj is evidently equivalent to 
two identities of type H201I) , namely for two linearly independent normal vectors 
n to f. Thus, we shall consider the two cases 

(i) n = T] : tangential to C(f ) (209) 
(ii) n — v : perpendicular to C(f ) 

The first case amounts to the characterization of T as a geodesic in the 
(truncated) cone surface C(f ) D Mh with the metric dst, as follows : 

Lemma 44 Let u{a) be the restriction of the potential function U* along the 
shape curve T*(a). Then the geodesic equation for the moduli curve T in the 
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cone surface C'(T) with the metric 



P 

is equivalent to the equation 



dsl = {-U* + h){dp 2 + p 2 do- 2 ) 



da da If. u(a) u'(cr) , 

+ — = — sina , y -\- cos a y ' (210) 



ds ds 2p \ u(er) + hp W ( <J ) + hp 

Proof. Let r\ be the normal vector in (|2U6[) . Then the geodesic condition is 
by (UD 

where /C(t/) is the (geodesic) curvature of T in the Euclidean sector (|204|l . How- 
ever, in a Euclidean plane it is easy to see that /C(ry) can be expressed as d(/ds, 
where C is the angle between a fixed reference ray a — (say, the positive x-axis) 
and the tangent line, in fact, £ = a + a, see Figure 8. Finally, calculation of the 
normal derivative on the right side of the identity l|211|l . using the orthonormal 
frame (|205|l . leads to the formula (|210l) . ■ 

In the second case of (|209|l the geodesic condition is the identity i|201l) with 
n equal to the normal v of the surface. In this case K,[y) equals the normal 
sectional curvature of C(T) in the direction of T, namely the value II(t, t) of the 
second fundamental form. The latter has the frame (|205|l as eigenvectors, with 
eigenvalues and K-*/p respectively, and hence by l|2U6|) and Euler's classical 
formula for the decomposition of normal geodesic curvature 



/C!sin 2 a = p/C(i/) = ln(-[T + /i) (212) 
2 dv p 



As a summary we now state the following theorem, valid as long as the 
quantities involved are well defined. 

Theorem 45 (cf. Theorem 6) In the moduli space M with the kinematic 
Riemannian metric ds 2 , let T be the oriented moduli curve of a three-body mo- 
tion with total energy h and vanishing angular momentum, and let T* be the 
corresponding shape curve on the sphere M* = S 2 (l/2) with unit tangent (resp. 
normal) vector t* (resp. u* ) so that {t*, v *} is a positive frame on the sphere. 
Then T = (p,T*) can be characterized as a solution of the following system of 
ODE 



r-\ v* ■ 2 1 d , nr*\ 

(il) : /C.sm a = — [mil ) 

9 2u(a) + hp dv* 



where JC* is the geodesic curvature ofT* in M* , a is the arc-length parameter 
ofT*, U* is the restriction of U to M* and u{a) is its further restriction along 
r* , and a £ [0, 7r] is the angle between the (outgoing) ray direction and T in M. 
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Proof. By using the expression for sin a in (|2U7[) . equation H21U|) can be stated 
as 

da 1 u(a) da 1 ^i 17 ) 

ds 2 u(a) + hp ds 2p u(a) + hp 

When we replace the arc-length parameter s of T by a, using a formula from 
l|2U7|l , this equation reads 

da , „ 1 u(a) . 1 w'fc) 
3r = (- 1 +o ...^ \J + ^ c ota- 



da 2u(a) + hp 2 u(a) + hp' 

and by viewing u' (a) / u(a) as the tangential derivative of ln(C/*) we obtain the 
first equation l|213l) . 

The second equation of Ij213(l is merely a reformulation of Ij212(l , whose right 
side may be expressed as 

1 -jLln( 1 -U* + k)= ( X\ ±U* 
2du* p u(a) + hpdv* 

Here we use the fact that the normal vector i> along T may be identified with the 
scaling of v* by the factor 1/p, that is, v — v* / 'p, and moreover, differentiation 
in the direction of u commutes with the scaling. ■ 



Remark 46 The above system is easily seen to be scaling invariant. 

Namely, when the size function p is multiplied by a fixed constant k > 0, the 
energy level changes as h — > h/k and hence the product hp stays invariant. In 
particular, since the energy level h = is invariant with respect to scaling of so- 
lutions, the explicit dependence on p in \21 SI) disappears in this case. Moreover, 
the angle a, geometrically interpreted in ^2Ub)) as the inclination angle of the 
moduli curve T, is a neat scaling invariant which together with T* represents T 
uniquely up to scaling. In fact, p is generally obtained from a by "quadrature" 
along T* , cf. ^Uti\) . This explains the following result. 

Corollary 47 Let the pair (a,T*) represent the moduli curve T of a three-body 
motion with vanishing angular momentum and vanishing total energy, where the 
shape curve T* is not a single point and is viewed as a curve on the standard 
sphere S 2 of radius 1. Then (a,T*) is a solution of the following system of ODE 

t = -3 + >ta^ln(L/*) 
/qsin 2 a= |^ln(^) (2W) 

where s — 2a is the arc-length parameter of T* on S 2 and IC* is its geodesic 
curvature. Moreover, a solution (a,T*) can only encounter a singularity (cusp) 
when a — or tt, or when T* reaches a collision point. 

We are particularly interested in applying the system (|213J) or (|214(l to the 
study of triple collision motions. A triple collision is simply expressed by the 
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condition p = 0, but this singular event is not explicitly visible in l|214|) since the 
variable p is eliminated. However, the term hp in (|213|l also disappears when 
p — > 0, so the two systems should behave "similarly" in the limit. Hence, the 
system (|214[1 is likely to be significant also when h ^ 0. 

One of the major results of Sundman and Siegel in their work on the local 
analysis of triple collisions prove the existence of both a limiting shape, neces- 
sarily a critical point of U* , and a limiting position, cf. |14| . |15j . |12j . The 
existence of a limiting position is the statement that the 3-body motion T(t) 
has a "size normalized" limit, 




(215) 



at the configuration space level, and we shall express the statement concerning 
the limiting shape (due to Sundman) by saying the pair (a,T*) approaches a 
specific pair (a, 5*), namely 

ae {0,tt} =0[O,tt] , re^^e^p+PoJcS 2 (216) 

It is also known (cf. e.g. Sieeel-Moser 1 1 3j . p. 89) that an Euler point 
can only be the limiting shape of a triple collision motion confined to a fixed 
line. (However, e, may well be the limiting shape of a non-collincar motion as 
t -> ±oo ). 

The two "boundary" values of a in i|216l) actually distinguish between the 
two events triple explosion and triple collision, as follows : a — when T starts 
(or "explodes") out from the cone vertex O of M , and lima = tt when T is 
oriented towards O and terminates with a "total collapse". Anyhow, we are 
free to run a three-body motion in either directions, and the associated initial 
value problem for (|213|) or (|214|l is (a priori) of singular type in the above case 
since sin a = 0. In fact, a solution (a,T*) of (|214|l may also encounter another 
type of singularity (called cusp) when sin a = 0, but with p ^ 0. In Chapter 8 
these events and related problems will be further investigated in selected testing 
cases. 

8 Case study of triple collision motions with 
zero energy 

8.1 The basic setting and statement of Theorem G 

Due to the simplicity of the system (|214|l . the special case of vanishing total 
energy, h = 0, lends itself as the simplest testing case of three-body motions 
leading to a triple collision. We shall investigate this case more carefully, and for 
convenience, let us also restrict ourselves to the case of equal masses, TOj = 1/3, 
which largely simplifies the series expansions of the potential function and its 
derivatives. 
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We shall address the triple collision problem as an initial value problem, 
namely as a triple explosion, although we usually write "triple collision" mo- 
tions. Let us first recall the so-called Lagrange- J acobi equation which is the 
result of differentiating I — p 2 twice with respect to time t, using the homo- 
geneity of U and conservation of energy h = T — U , namely in our case 

d 2 

—I = 2(T + h) = 2T>0 (217) 
dt z 

It follows that / is a nonnegative convex function of time t, and starting from a 
triple collision (say, 1(0) — 0) it is strictly increasing and tends to oo as t — > oo. 

Following the setup from Section 7.3, we seek a description of the moduli 
curves of triple collision motions, valid for some appropriate time interval [0, t%\. 
For this purpose it is convenient to use the coordinates (p, ip, 9) in the cone M 
= C(M*), where as before p = y/l and (ip,0) are spherical coordinates on the 
unit sphere S 2 = M* centered at the north pole AT. In this setting a moduli 
curve r and the associated shape curve T* have coordinate representations 

f (s) = (p(s), <p( a ), 6(a)), T*( s ) = (<p(8),0( 8 )), (218) 

where s = s(t) is the arc-length parameter s = s(t) of T* and is an increasing 
function of time t. Here we must exclude, of course, the well understood shape 
invariant motions, namely the trivial case that T* is a single point (in which 
case p(t) is the solution of a 1-dimensional Kepler problem). 

Thus, we seek a description of all those moduli curves T(s) emanating from a 
triple collision, at s = say. According to Remark 0^] it suffices to consider the 
class of r modulo scaling, represented by the pair (a,T*) where a(s) E [0, 7r] 
is the (radial) inclination angle of f. In fact, recall from (|208|l that the size 
function p of T is recovered from (a, T*) by the general quadrature formula 

p(s)=p( So )e^°° cot{a)ds , p(s o )^0 (219) 

Remark 48 Using the parameter s rather than time t is, of course, crucial for 
our geometric approach below. The relationship between s and t is follows from 
the kinematic metric, using e.g. \204\ ), j207]) , \261$ . Namely, in the case 

of zero angular momentum there are the identities 

2 2 

2(U + h)dt 2 = 2Tdt 2 = ds 2 = dp 2 + ?-ds 2 = (cos 2 a)ds 2 + ?-ds 2 , 
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from which we deduce the relationship 



dt = = ds (220) 

2 3 / 2 sin a(s) ^u(s)/p{s) + h 

where u(s) — J7*(r*(s)). Moreover, by switching over to t it is, in fact, not 
difficult to see that p (at any time to) can be determined solely from the time 
parametrized shape curve T*(t) and the normal derivative of U* (near t = to). 
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Now, resuming the assumption h = 0, our approach is to determine the 
above pairs (a,T*) by solving the system 

ODE* : | /C;sin 2 a- i^ln([/*) (221) 
| l = (^ + ( si nV)(f) 2 

as an appropriate initial value problem which represents a triple collision, sec 
(JH21). The system <|2"2"Tj) is a copy of (I214|) since the third equation merely 
expresses the constraint that V* is a curve on the unit sphere. We will refer to 
the first and second equation of (|221|l as the inclination and curvature equation 
respectively. The first one relates the growth of the inclination angle a (of the 
moduli curve f ) with the tangential derivative of ln(C/*) along T*, whereas the 
second one - which is of order two- relates a to the geodesic curvature of T* and 
the normal derivative of \n(U*) along T*. 

The following theorem summarizes the main result of this chapter. It de- 
scribes the family 6(po) of all shape curves T*(s) representing triple collision 
motions, with the limiting shape of po at the collision. 

Theorem Gi In the case of uniform mass distribution and zero total energy, 
consider the family S(po) of arc-length parametrized shape curves T*(s), s > 0, 
which emanate from the north pole po = T*(0) of the 2-sphere S 2 and represent 
3-body motions with a triple collision at s = 0. This family has the following 
properties : 

(i) There is a unique curve Y* 6 for each initial longitude direction 9q. 

(ii) The family is invariant under the induced action of the dihedral isometry 
group D3 of S 2 which fixes po and permutes the three Euler points e%. In 
particular, r| o+27r / 3 is obtained from T* ea by rotating the sphere, 6 — » 9 + 2ir/3. 

(iii) Each curve stays within a sector of angular width 7r/3 and bounded by 
meridians representing the shape of isosceles triangles, at least until the first 
eclipse (i.e. crossing the equator circle). 

(iv) Each curve extends analytically through s = and Tg Q+K (s) = Tg Q (— s), 
and it has no singularity before the first eclipse. 

(v) For each 9o the associated inclination angle function ag (s) is the unique 
analytic solution of y221\) along Tq , with the singular initial condition a(0) = 
0, a'(O) > 0. Moreover, as + 7T {s) = — cte (— s). 

(vi) The sign of the curvature of the curves in 6(po) depends only on the 
sector, and in neigboring sectors the sign is opposite. 

Remark 49 In the case of rectilinear three-body motions, the corresponding 
sets S(ej) are obviously "congruent", each consisting of the pair of arcs 
of the equator circle, in opposite directions and starting at the Euler point e^ . 
The associated inclination angle function a±(s) is defined by a unique analytic 
function a(s) so that a±(s) = a(s) for s > 0, and a(—s) — —a(s). We refer to 
Section 8.4- 
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We also refer to Section 8.6.2 for more information concerning the geometric 
behavior of the curves in S(po). Indeed, we are actually close to a stronger 
version of Theorem Gi, but the proof needs more elaboration, so we formulate 
the following conjecture as an open problem. 

Conjecture 50 The different triple collison shape curves T*(s), s > 0, intersect 
the equator circle the first time at different points, and moreover, each point on 
the circle is reached by a unique curve. The curves do not intersect each other, 
except possibly after the first eclipse. 

Corollary 51 Under the current hypothesis of uniform mass distribution and 
zero total energy, consider the "moduli space" consisting of all three-body mo- 
tions in 3-space which start from a triple explosion at time t = 0, and moreover, 
the motion is neither rectilinear nor shape invariant (i.e. homographic) . This 
space can be naturally identified with the manifold 

50(3) x 50(2) x M+ 

In particular, the "moduli space" for those triple collision three-body motions 
confined to a fixed plane is 

0(2) x SO(2) x R+ 

Indeed, starting with the space 6(po) — SO(2) of shape curves described 
in Theorem Gi, 5*0(2) x R+ is the space of their associated curves f(s) in 
the moduli space M since the size function p(s) can be scaled by any positive 
number A S M + without affecting the shape curve. Next, the possible liftings 
T(s) of f(s) to three-body motions with T(0) = are distinguished by the 
normalized limit 8 in (|215|l , which is an oriented, regular m-triangle A of unit 
size. We may identify the various positions of such an oriented triangle with 
the rotation group 50(3) which measures its "deviation" from a fixed reference 
position. In terms of the fibration 

0(2) -> 50(3) -> 50(3)/0(2) ~ RP 2 

we can say that the projective plane RP 2 represents the choices of 2-planes 
containing A (and the motion), whereas 0(2) represents the possible positions 
of an oriented regular triangle in a given plane. However, we mention that there 
is no global "field" of reference positions (or gauge) for all the planes, since this 
would imply the fibration is trivial, which is certainly not true. 

8.2 Analysis of the potential function for equal masses 

In this chapter we shall choose the zero meridian 9 — for the polar coordinate 
system (<p,9) of 5 2 different from the convention in Remark El Namely, the 
three binary collision points hi,i — 1,2,3, which are now equally spaced along 
on the equator circle ip — 7r/2, will have the longitude angles 

0i = ™, 2 = I, 9 3 = tt (cf. Figure 9) (222) 
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Note, for example, the antipodal point of hi is the Euler point ej, and now 
9 = at e 3 . It is also convenient to use negative values of <p, with the usual 
interpretation so that (ip, 0) and (— ip, 9 + it) is the same point on the 2-sphcrc. 
This is consistent with our trigonometric formulae for U* (ip, 9) below, see (122411 
and j2jg) . 

For convenience, let us normalize U* by a constant factor to make its mini- 
mum value U*(f>o) = 1. In fact, scaling of U* has no effect on the system <|221[) . 
Thus, for peS 2 

i—1 



i=i p - b 

where (for any mass distribution, indeed) 



P b, 



= %/2(l - z,) 1/2 , with *i = sin <£cos(0 - 6i), (224) 



is the usual Euclidean distance from p to the binary collision point hi and 9i is 
the longitude angle of h t , cf. (TIT)) . (IS5jl . P7T)I.I555I 

Consequences of the invariance of [/* with respect to permutation of the 
points hi will be analyzed and exploited later (cf. Section 8.3.1). At the alge- 
braic level, however, the following symmctrization technique will facilitate the 
analysis of U* and related series expansions. For each integer k > 0, define 

3 

S k = ^cos k (6 -9t) 

i=l 

and write 

f(x) = f[(x - cos(6> - 9{)) = x 3 - ^x + i cos(36>) 

i=l 

Logarithmic differentiation of f(x) leads to the formal identity 



2 



from which we deduce the recursive formula 

Sfc+a = |s fc+ i - i cos(30)S fc , fc> (225) 
For convenience, the first few Sfc are listed as follows : 

So = 3, S 1 = 0, S 2 = ^, S 3 = -~cos(30), S 4 = ^ 

15 27 3 63 

^ = --008(30), S e = - + -cos 2 (39), S 7 = --cos(W) 
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It follows, for example, that Sk as a polynomial in cos(3#) has all its nonzero 
coefficients positive (resp. negative) when k is even (resp. k is odd). By 
expanding U* as a sum of binomial series in the variables z, and using the 
identity 

3 

8=1 

we arrive at the following trigonometric series 

1 00 /-IN 

U* = U*(<p,e) = -J2( 5 (-l) fc (sin fc ^)S & (226) 

3 5 
= 1 + — sin 2 tp - — (sin 3 (cos 30) + ...+ 
16 64 

For later use we also introduce the following functions on the sphere 

F{P) = ^±f<^, G(p) = ^±s^^ (227) 

i=l p-bj i=l p-bj 

It follows that 

-jfir = —Fsimp, ^- = Gcosip (228) 

ad dip 

Remark 52 It is easy to check that F vanishes precisely along the six meridians 

r* : 9 = fe|, < k < 5 (229) 

passing through either a binary collision point hi or an Euler point (= — b,-). 
Hence, F changes its sign by crossing these meridians, but on the other hand, 
the function G is positive (except undefined at poj. For example, F is negative 
for < 6 < 7r/3 ; and in this sector the gradient flow of U* is depicted in Figure 
9. 



8.3 Reduction, regularity and singularity 

Here we shall describe a finite group acting on moduli curves and, in particular, 
it is a symmetry group of the space of solutions s — * (a(s), T*(s)) of l|221|) . 
Moreover, regularity and singularity aspects of the solutions we seek are also 
briefly discussed. 

8.3.1 Discrete symmetries and reduction 

In addition to time translation and space-time scaling symmetries which trans- 
form solutions of the general 3-body problem as in H197|) ~). there is also an 
additional symmetry group of order 4 which we denote by 

2ixZ 2 = {a, f ) ~ Z 2 x Z 2 (230) 
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The involution r represents reversal of time, t — ► — t , and its induced action on 
oriented curves in the moduli space is expressed by 

f : (s, a) — > (s, a) = (— s, 7r — a) (reversal of direction) (231) 

which takes a solution (a(s), T*(s)) of the system (|221[) on the interval si < s < 
S2 to the "reverse" solution in the opposite direction and defined on the interval 
— S2 < s < — si(or any translation of this interval). The other involution a 
is a purely geometric symmetry, arising from the reversal of orientation of m- 
triangles. At the moduli space level, the latter is the reflection of M ~ M 3 in 
the (equator) xy-plane, that is, the map ip — > tt — ip in the coordinates (p, ip, 8). 

On the other hand, under the present assumption of equal masses, there is 
also the (order 6) dihedral isometry group D3 C 0(2) of S 2 which fixes the poles 
p and permutes the Euler points e,. It is a symmetry group of the 3-body 
problem since it leaves U* invariant, cf. Section 8.2. The action is generated by 
the rotation 8^8 + 2tt/3 and the reflection 8 — > —9, and altogether we have a 
symmetry group of order 24, 

© = D 6 x Z 2 = (£> 3 x £>i) x Z 2 (232) 

where we also regard De C 0(3) as a dihedral isometry group of S" 2 generated by 
reflections. Thus, the action of 'Sq divides the sphere into 12 congruent spherical 
triangles called chambers, and we choose one of them to be our fundamental 
chamber, namely the following geodesic triangle on the upper hemisphere 

€ = {(p,8)eS 2 ; 0<p<^, 0<8<^} (233) 

(cf. Figure 9) with the vertices 

po = (ip = 0), e 3 = O = tt/2, 8 = 0), b 2 = [p = tt/2, 8 = tt/3) (234) 

In particular, the action of on solutions (a, T*) of (|221|l reduces the study 
of solutions to the study of "solution segments" T* inside Co- ln particular, we 
may restrict the study of triple collision solutions T* to those emanating from the 
vertex po or 63 with initial direction leading into the chamber Co- The natural 
first step of this program is to look for solutions whose curvature equation in 
H221|) is trivially satisfied, and Section 8.4 is devoted to this preliminary study. 

8.3.2 Cusps and other singularities 

It is well known that a 3-body motion T(t) can be "regularized" through a 
binary collision (cf. ^31)i an d the only "real" singularity must be a triple 
collision. Away from collisions the moduli curve T and the shape curve T* are 
also analytic functions when we parametrize by time t or the arc-length s of 
r. However, we also want to parametrize by the arc-length s of T*, and then 
singularities may occur at specific instants where the time derivative s(t) > 
vanishes. Although this type of " singularity" is rather artificial, it has geometric 
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significance which explains the possible cusps of the embedded curve T* on the 
2-sphere. These are also singularities of the system (|221ll . and in this subsection 
we shall discuss them in some detail. 

Since 4| > for all t, the event s(ti) = is equivalent to the condition 
a(ti) — or 7r. In this case, -^T*(ti) — and we say p = T*(ti) is a halting 
point for the shape curve. Geometrically, this can be a singular point for T*, 
namely it is the type of singularity that may occur when a regular curve ip(t) 
in 3-space with vertical tangent at t = t\ is projected to a curve i^*{t) in the 
xy-plane, cf. also (|2 18|> . On the other hand, when T* is parametrized by s and 
si = s(ti), the unit tangent vector 4-T*(si) still exists (as a one-sided limit) at 
the halting point p. Moreover, the (one-sided) geodesic curvature K.* of T* will 
be bounded near p. In fact, for a "thin" region < p 1 < p < p 1 + Sp of the 
moduli space (M,ds 2 ) with the kinematic metric l|193|) . the projection to the 
shape space (M*,ds 2 ) may be viewed as a Riemannian submersion modulo an 
almost constant scaling. Restricting to the above region, the curvature of the 
moduli curve T is certainly bounded, and its image curve in M* will also have 
bounded curvature. 

Now, consider a pair (a(s), T*(s)) which is a solution of the system l|221|) . 
The pair is regular on the interval (si, S2) if the three functions a(s), y>{s), 9(s) 
are analytic and a(s) ^ 0, 7r, and a singularity is encountered at s — s, ; if we 
cannot extend the functions regularly beyond this point. In that case it follows 
from l|221|l that either sin a(si) = 0, in which case we call p = r*(s,) a cusp, 
or p is a binary collision point bj (in which case a(si) = n/2 and a'(si) = 00). 
On the interval (si, S2) the growth of a(s) is governed by the inclination angle 
equation (cf. (|221|l ) 

da 11 . . , . . 

^ = -^ + 5 cot <*(«)£>(«) (235) 

whose dependence on the shape curve is solely through the tangential logarith- 
mic derivative 

u'(s) 

£>(s)=V(lnt/*)-T* = -^, u(s) = U*(T*(s)) (236) 
u{s) 

Assume there is a halting point at s = Sj, say a(si) = and write p = r*(sj). 
If D(si) ^ 0, then by (j235(l a'(si) = ±00 and p is a cusp. On the other hand, 
if D(si) = and p is not a critical point of U*, then the curvature (i.e. second) 
equation of 1)22 l[l . whose right side is nonzero at s = s,, would force the geodesic 
curvature of T* to become infinitely large towards p. However, as observed 
above, such a behavior of the shape curve is not possible. Hence, D(si) = is 
only possible when the halting point p of T* is also a critical point of U* . 

Finally, assume the halting point p is also a critical point of U* . For s close 
to Si we have cot a ~ 1/a, and the local behavior of a(s) near s = Sj is largely 
governed by equation l|235|l . from which we can show a'(s) is bounded in a 
neighborhood of s,. Choose some sq so that si < sq < S2 and consider the two 
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cases depending on whether s is approaching Si from above or below 



ds 

i) Si — si : I cot a(s) ds ~ I — — = — oo (237) 
J so J so a i s ) 

fS2 

ii) St — S2 '■ / cota(s) ds — oo (238) 

J so 

In particular, in case i) an associated 3-body motion must necessarily encounter 
a triple collision at s = s\ since formula (|219|l implies p{s\) = 0, whereas in 
case ii) p(s2) = oo and p = T*(s2) is the limiting shape of the 3-body motion as 
t — > oo. In contrast to this, for a cusp singularity with < p{si) < oo, we have 
a(s) — ► and a'(s) — > ±oo as s — + Si, in such a way that the integrals (|237|l or 
1238(1 will converge. For example, this would be the case if a(s) <~ k(s — Si) p for 
some constant k and p < 1. 

We claim, in fact, that (a(s), T*(s)) is regularizable at the above halting 
point p = r*(sj) which is also a critical point of U*, that is, as a solution of 
(12211) the functions can be extended analytically beyond S{. Indeed, once the 
derivative a'(si) exists, the power series developments at s — Sj of the functions 
a(s), y(s), 9(s) are recursively determined from the system (|221|l . The recursive 
scheme is worked out in detail in Section 8.4 for the case (|237l) with a triple 
collision at the north pole po = T*(si), and case ii) is similar. 

The collinear type of triple collision is at an Euler point such as p = §3, 
in which case T*(s) moves along the equator. In particular, u(s) in (|236l) has 
a minimum at s = si. In fact, approaching §3 from other directions (not 
tangential to the equator) would force a'(s) to become complex valued. It is 
remarkable that a'(si) — clq (resp. 60) turns out to be the same constant for all 
possible triple collision curves emanating from pjj" (resp. e,), see (|25U|I . 

In order to describe analytically the regularization of the triple collision 
motions it is convenient to extend the domain of the angle a to negative values 
as well. Indeed, —a should be identified with tt — a, and therefore we introduce 
the a-circle 

M ~ 51 c C : a -» z a = e 2M (239) 

(0~7T) 

as the new domain for a. Then the continuous motion z a (s) on the circle 
illustrates the qualitative behavior of the solution (a(s),T*(s)) and hence also 
the moduli curve T(s) of an associated 3-body motion. For example, p(s) is 
increasing (resp. decreasing) with s when z a lies on the upper (resp. lower) 
semicircle. A halting point is characterized by z a = 1, and it represents either 
a cusp, a triple collision (p — * 0), or an escape (p — * 00). 

The only way z a (s) enters the other semicircle is at z a = 1, but for p — > 
or 00, since at a cusp z a "bounces back" on the same semicircle. Similarly, z a 
reaches the value —1 at a binary collision point hi, but again z a "bounces back" 
(if the moduli curve is continued via regularization with p increasing ). 

Summary 53 The solution (z a (s),T*(s)) has a singularity at s — Si if either 
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i) the pair takes the value (l,p) where p does not belong to the set 

{±ft),Si,§a,63}u{bi,b3,b 3 } ) (240) 

or ii) it takes the value (— l,bj),i > 0. The singularity is a cusp (resp. a binary 
collision) in the first (resp. second) case. 

8.3.3 The initial value problem for triple collision solutions 

For later reference we introduce the set S(p) of all shape curves T*(s), s > 0, 
representing a triple collision motion with the limiting shape p = T*(0) at the 
collision. Here p can be any of the five points of the first subset in l|24(J|) . In fact, 
each r*(s) is associated with a unique (inclination angle) function a(s) so that 
the pair (a,T*) is a solution of the system ODE* in (|221l) with the additional 
and singular initial condition 

i) a(0) = 0, ii) o/(0) > (241) 

In particular, for a fixed T*(s), a(s) is the unique solution of (|221|l and l|241|) . 
Thus we may as well consider the totality of pairs (a,T*) and define 

6(p) = {(«(*), r»); s > 0,r*(0) = p, a(0) = 0,a'(0) > 0} (242) 

as the solutions of a specific initial value problem for the system ODE*, as 
indicated. 

For the calculation of the sets (|242|l we may assume (by symmetry) the initial 
shape p belongs to the fundamental chamber £ , namely p is either its vertex 
S3 or po, see (|234H . Observe that the set (|242|l has the induced symmetry group 
(& p which is the "isotropy" subgroup at p of © (cf. (1232(1 1 

(S e3 ={l,f}x{l,f}~2) 1 xZ 2) ^=3 3 xZ 2 (243) 

where r £ D3 is the reflection 8^—9. 

Remark 54 Contrary to the above, the initial value problem for \221)) at a cusp 
or binary collision is not well defined since we would have a'(-Si) = ±00. There- 
fore, one cannot continue a solution (a(s),T* (s)) across the singularity using 
only the system ODE* . To circumvent the problem one should, for example, 
turn to the moduli curve T itself, which is regular in any case. In Section 8.4-2 
below we shall use Newton's equation of motion more directly to develop in time 
a specific solution through several cusps. 

8.4 Isosceles and collinear triple collision motions 

In this section we take the opportunity to illustrate the above approach applied 
to the " simplest" type of triple collision motions apart from the shape invariant 
ones. We shall also supply with numerical calculations, for comparison reasons 
and illustration of examples only. 
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Namely, consider the possibility that the shape curve T* of a triple collision 
motion is confined to a great circle on the sphere, that is, JC*(s) — for all s. 
Then the curvature equation in l|221|l forces the normal derivative of U* along 
T* to vanish. Hence, by l|228|) and Remark|221 T* "moves" either on the equator 
circle (<p — it/2) or on one of the six meridians (|229|l passing through an Euler 
point e, or a binary collision point hi. These meridians represent the shape of 
isosceles triangles, so the associated 3-body motions are either of collinear type 
or isosceles triangle type. 

By the symmetry reduction explained in Section 8.3.1 it suffices to consider 
three separate cases, namely T* is (initially) confined to one of the boundary 
arcs of the fundamental chamber (|233[) . We list the starting point (at the triple 
collision) and a choice of arcdength parameter s to be used (up to the first cusp 
or binary collision point) in each case : 

(1) po : 6 = 0, s = ip > 

(2) e 3 : ip = tt/2, s = 9 > (244) 

(3) po : 9 = vr/3, s = Lp > 



8.4.1 The inclination angle a and its ODE 

We shall investigate the first equation of 1221|) 

^ = -~ + ~cot(a)A(s) , * = 1,2,3 (245) 
as 4 2 

for each of the three cases (|244|) . where -Dj(s) is calculated from the appropriate 
expression of U* = U*(Lp, 9) = U*(lp, 9 + 2tt/3), namely (cf. 

U*{tp,Q)= \{ . 2 + w ^=. ), (246) 

U*(n/2,6) = -( , 1 + 1 — + 1 ) 

3 v v /l-cos(6' - f) y/1 - cos(<9 + f ) VI + cos 9 

U*(<p,ir/3) = U*(-cp,0) = U*(<p,ir) 

As will be demonstrated below, there is a unique solution of the initial value 
problem l|24ip . 

The derivatives Dj are the following analytic functions expanded at the point 
po in case (1) and (3), and e 3 in case (2) : 

n ( s ,3 15 23 2 

DM = a^'frAf) = 6{ ™ + + + ...) (247) 

y ' U*{n/2,9) v 20 1200 12000 ' K ' 

n . , j^W/g) 3 15 23 

D ^ = UH^j3) = ^ ( 8 + 64^ + 256^ + ^ 
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It follows that Di(ip) = —Ds(—<p) and hence case (3) of (I244[l can be subsumed 
under case (1) by using the range tp < and moreover, a < interpreted in 
accordance with (I239[l . In fact, if ai(ip) is the solution of i|245|) for i = 1,3, then 
a>3(ip) = —ai(—ifi). Hence, we need only consider the cases (1) and (2) of l|245|) . 

To investigate the nature of the singularity a = 0, let us first approximate 
the functions D\{ip), D 2 {0) and 

-i 1 1 ^ 

cota = a a or + 

3 45 

by their first term |</?, ^8 and a -1 respectively. Then the initial value problem 
a(0) = for the simplified version of l|245|) has the two straight line solutions 

(1) : a = a ip, a Q = — 1 (248) 



(2) : a = b 6, b = -tL- 

found by solving a second order polynomial, namely 

1 3 1 1 29 1 

4 16 ao 4 40 b 

However, the extended initial value condition (|241|l demands a to be initially in- 
creasing and hence selects the positive solution in (|248|l as the leading coefficient 
for the two types of triple collision, namely 

^/^3 \ 

Lagrange type : a = « 0.326 (250) 

8 

±vTl85- 1 
Euler type : b Q = ^ w 0.736 



Returning to the original equation (|245|l we can determine recursively the 
power series expansion of a, 



(1) a — a ip(l + ai<p + a 2 ip 2 + ...) (251) 

(2) a = b o 6(l + b 1 + b 2 2 + ...) 

which depends solely on the initial oefficients eto,6o m l|250l) . For convenience 
we list (with a few decimals only) the first terms of the expansions : 

(1) a « (^(0.3257 - 0.0955<^ + 0.0129(^ 2 - 0.0233</? 3 + ...) 

(2) a = 0(0.7356 + 0.19416* 2 + 0.04876* 4 + ...) 

As indicated, in case (1) the series is alternating and in case (2) a(s) is an odd 
function since bi = 63 = 65 = ... = 0. 
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The above series can be easily developed and used with high accuracy for 
small ip (or 9). However, we shall only use it to calculate an initial value of a 
for some small tp (or 9) and then solve equation <|245|) by standard numerical 
procedures, on the maximal interval bounded by the first singularity in each 
direction. Namely, in case (1) we find 

lim a(tp) = -it/2, lim a'(tp) = oo (252) 

if — > — 7T/2 ip—> — 7rf'2 

lim a(<p) = 0, lim a'{<p) = -oo, ip t w 1.876 w 107.5° 

On the interval [0,^], a(ip) increases up to its maximum at ip w 1.2 and is 
thereafter decreasing, see Figure 10. 

On the other hand, if we rotate by 180° the graph of a(ip) over the interval 
[— 7r/2,0], then we obtain the graph of a(cp) on [0, ir/2] for case (3) of l|244|l . 
Finally, the graph of a(9) in case (2) on the interval [0,7r/3] is quite similar to 
that of case (3), see Figure 11. Its graph over [— 7r/3,7r/3] is symmetric with 
respect to the origin. 

As calculated above, the shape curve T* along the meridian 9 — reaches 
the first cusp at <p = (p±, which is beyond §3. At the cusp the motion T*(s) 
changes its direction and continues northward and across e 3 . See Section 8.4.2 
and the time development of this motion. 

Here is a brief summary of the analysis of the preliminary cases l|244f) and 
the solution sets (|242|l to which they belong : 

• The set 6(63) has only two solutions (a+,ri), (a_,r^), and they are 
equivalent modulo the isometric reflection f : 9 — > —9 belonging to Q3 D 
C3e 3 - Let r*(0), 9 £ (— 7r/3, 7r/3), be the arc-length parametrization of the 
equator circle from bi to B2. There is an analytic function a(9) so that 
the pair (a, T*) is a solution of (|221|l and moreover, 

at-(s) = — a(— s) for s = —9 > 0; oq-(s) = a(s) for s = 8 > 
T*_(s) = r*(-s) for s = -9 > 0; T* + {s) = T*{s) for s = 6 > 

• The set S(po) has a unique solution (ao, TJ) and T* / 3 ) in case (1) 
and (3) of (|244|l . respectively. Consider also the solution 

MK /3 ,r; /3 ) = K,r;) e 6(p ) 

along the meridian 9 = tt, obtained by applying the rotation p, £ D 3 : 
9 — > 9 + 27r/3, as indicated, and let r*(<^),<^ e (— 7r/2, tt/2), be the arc- 
length parametrization of the half-circle (b 3 — * p — * e 3 ). There is an 
analytic function a(ip) so that the pair (a,T*) is a solution of (|221ll and 
moreover, 

a T (s) = — a{— s) for s = —ip > 0; ao(s) = a(s) for s = ip > 

r;(s) = r*(- s ) for s = > o ; rg(s) = r*(s) for s = <^ > o 
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8.4.2 Time dependence and Newton's equation 

We shall choose case (1) of l|244[) and compare the above approach using the 
system (|221H with the time parametrized motion T(t) using Newton's equation 
Namely, in the xy-plane we consider the motion of three point masses of 
mass 1/3, symmetric with respect to the y-axis and with position vectors 

ai = (-x,y),a 2 = (x,y),a 3 = (0,-2y) 

Newton's equation (JIJ reads 

■ ±= 1 N 1 jj = y (253) 

12 x 3 3 {x 2 + 9y 2 ) 3 /2 ' y ( x 2 + 92/ 2)3/2 V ; 

As initial condition at time t — 0, assume y = (i.e. T* is at the Euler 
point §3) and moment of inertia I = p 2 = 1. Moreover, let (3 denote the oriented 
angle from the positive x-axis to the initial velocity vector a.2, whose length is 
denoted v. Assuming the total energy h = T —U vanishes, the initial condition 
now reads 



(•'■• </)|/-o = (\/|,0), (x,y)\ t=0 = v(cos/3, sin/3), z; = ^ 



5 2 

" V1 (254) 



1 + 2 sin^ f3 



Finally, let us assume < (3 < it, which means the shape curve T* is heading 
southwards from 63 for small t > 0. 

The above data specify a 1-parameter family (parametrized by (3) of isosceles 
3-body motions T(t) — (x(t),y(t)) with total energy h = 0, with normalized 
size 1=1 and collinear shape S3 at time t = 0. We shall use the equation 
Ij253(l to investigate the time dependence of the various geometric and kinematic 
quantities of the motion, such as /, T, (p, a, where 

I=-x 2 + 2y 2 , C os<p = ±^=- = ^L 9 f V - (255) 
3 W ^ 3\/3 / V3 |.x 2 + 2y 2 V ; 

r = iv = l ( A + 7 2 2 § 2 ) = \tf + &) 

By definition of the inclination angle a, 

coca- ^ r(t) - g = M+M _ 
||r(t)| V2T V2T(7 ' 

and denoting the angle at t = by ao we deduce 

cosao = ±(H)V4 WC0S/ 9 = - C0S/? (257) 
V5 3 v/l + 2sin 2 /3 

Thus, the correspondence (3 < — ► ao is a bijection of the interval (0, 7r) such that 
n — p corresponds to 7r — ao- 



81 



We claim there are exactly two values of (3 leading to a triple collision motion, 
namely the pair (3 and n — /3 for some /3 < n/2. This choice of f3 = f3 yields 
p(0) > 0, and hence the triple collision occurred in the past, namely at some 
negative time to < 0, with the shape curve T* at the north pole po- Similarly, 
using (3 = it — /3 the triple collision is reached at time —to > with T* at the 
south pole. 

The angle f3 is calculated using the formula (|257|l . where ao = a(ir/2) and 
a(ip) is the solution of the equation 

da ip , ,,3 15 23 9 , ,„_„n 

^-- 1/4+ f COt(a)( 8-64^ + 256^ -•- ) (258) 

with initial condition a(0) = 0, cf. case (1) of I245|l and l|247|) . The solution 
found by the approach in Section 8.4.1, is approximately 

ao = a{n/2) » 0.18673..., f3 = 0.10865.... 

Thus, we know the initial data l|254|) corresponding to triple collision motions. 
As a test, by running the system l|253|) backwards in time one will find that the 
triple collision occurs approximately at to = —1.0228... 

It is also interesting to follow the shape curve T* of T(t) for t > 0, for 
example, using the ratio y(t)/x(t) or calculating (p(t) directly from l|255[) . The 
solution r(i) can be continued in time t through the cusps since they are not 
singularities for Newton's equation, and only truncation errors or numerical 
instability may invalidate the calculation in the long run. At t\ ~ 10.4, y/x is 
maximal and (p{t\) = (p 1 (cf. (|252|l ) is the colatitude of the first cusp - here a 
= and r* (t) turns northward. After passing 63 there is a second cusp where 
r* turns southward again and crosses 63, but the next cusp is closer to §3, and 
so on. Concerning the possible asymptotic behavior of T* (t) , at triple collision 
or as t — > 00, see also Section 8.6.3. 

Finally, we consider the time dependent size function p(t) — \J I(t) and 
compare it with the integral formula l|219() . We have, by assumption, p = 1 at 
time t = 0, and Newton's equation (|253J) yields, for example, 

p(-0.5) w 0.634726, p(-0.5) = (p « 1.381793 

On the other hand, numerical integration of the solution a(ip) of Ij258(l yields 

p\ ^ = e i / * cot adip sa 0.634725 (259) 

Alternatively, let us also evaluate this integral using the time parametrized 
function a(t) calculated by (|256|l and developed via Newton's equation. Thus 
we change the variable ip in 1259|) to t using (|255|l . namely 

dip = tp'(t)dt = (4arccos( ^{^vi 1 ) ^ 

and then numerical integration similar to the case Ij259(l yields 
p(-0.5) = e Ho"°- 5 ^(a) v '(t)dt „ 0.634726 
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8.5 Analytic uniqueness of triple collision motions 

In this section we turn to the full family 6(po) of triple collision solutions (a, T*) 
of the system ODE* , with T* starting out from p . Due to the symmetry group 
2)3 we may assume the shape curve T* enters the fundamental chamber C > 
which limits the initial direction 9q to the range [0, 7r/3]. Therefore, in terms of 
spherical coordinates (if, 0) the appropriate and complete initial condition (|242ll 
now reads 

a(0) = 0, a'(0) > 0; ip{0) = 0, 0(0) = 6 , < 9 < tt/3 

The border cases 0q = 0,7r/3 are the cases (1) and (3) of (|244[) already 
investigated in Section 8.4, and now it is natural to generalize the procedure 
used there to the whole range of angles 8q. This time, however, the strength 
of the curvature equation of (|221[) must be fully utilized. At this point, we 
assume (tentatively) that the functions a,tp,9 have power series expansions at 
Po, necessarily of type 

Ifi = s(c + CiS + c 2 s 2 + c 3 s 3 + ... + ) 

9 = 9 + s(d + d lS + d 2 s 2 + ...+), 0<6» <^, (260) 

a = a a s(l + ais + a 2 s 2 + •••+) 

Recall from Section 8.1, we have excluded the trivial case of constant shape, 
and hence a does not vanish identically. To justify the notation in the third line, 
it will be demonstrated below that the leading term is oos with a>o =/= 0. 

By considering the leading coefficients of the series for ip and 9 the third 
equation of ODE* implies 

co = l, ci = 0, c 2 = -^d 2 <0 (261) 

and clearly 

(sin^)0' = e(l-^' 2 ) 1 / 2 = e(-6c 2S 2 + ...+) 1 / 2 , e = ±1 (262) 

The calculation of ciq in the expansion a — oqs + ... is really the same as in 
Section 8.4.1 and gives the same value (|250|l independent of 9q. To see this, we 
write for clarity the first equation of ODE* as 

a(2a' + -) = (a cot a)D(s), (263) 

where 

112 1 

a cot a = 1 a 2 a 4 a 6 a 8 + 0(a 10 ), 

3 45 945 4725 v h 

and the potential function (|226|l and its logarithmic derivative along Y* have 

the expansions 

00 

u(s) = U*(T*(s)) = VW = 1 + —s 2 - — (cos30 o )s 3 + •- (264) 
^— ' 16 64 

i=0 
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d , , , u'(s) 1 ,dU* . dU* .. 

D(s) = - In (u) = -^f = -(-o-v' + -57-0 
as u(s) u at/3 aw 

00 3 1 5 

= S £ ^ = s{~ (- cos30 o )s + -..) (265) 

«=o 

In particular, the first order term in (|265|l is independent of 9q and the leading 
terms of l|263f) yield the single condition 

2 / V13- 1 , . 

4a + ao — 3/4 = , with positive root : ao = (266) 

8 

The identity l|263|) also provides recursive relations for the calculation of cifc, k > 
0, expressed in terms of the coefficients n i; i < k, see (|277|l below. 

8.5.1 The method of undetermined coefficients 

The proof of Theorem Gi , concerning the existence and uniqueness of the curves, 
is based upon formal power series substitution for the three functions (|26UI) in- 
volved in ODE*. This leads to a recursive procedure - the method of undeter- 
mined coefficients - which is consistent and determines successively the higher 
order coefficients in terms of C2, do and 6q. By (|261|l c-i is already determined 
by do, and at the final stage we shall find that do is actually determined by 
#o- Consequently, the expansions in 126U|) are, indeed, determined by the initial 
angle 9 alone and hence 9 parametrizes the whole solution set 6(p ). 

On the 2-sphere there is the positive, orthonormal frame ~^p'§g} as ~ 
sociated with the coordinates tp, 9. Along the oriented shape curve T* we also 
have the positive, orthonormal moving frame {t* , v*}, where t* is the tangent 
vector. The latter frame differs from the stationary frame by a rotation angle 
/3, namely in analogy with (|206jl - (|207jl . 

d sin/3 d „ . d cos/3 d 

t = cos/3— + - — , v =-sm/3— + - — (267) 

op sin ip 00 op sm ip o9 

cos/3 = pi', sin/3 = (smcp)6' = e(l - p>' 2 ) 1/2 , cf. G53 (268) 

Now, we turn to the second equation of Ij221|l , namely the curvature equation 
written as an identity 

LHS = RHS (269) 
between the left hand and right hand side 

LHS = (1 - cos2a)K* g U* (270) 

dip sm ip 09 
In LHS the geodesic curvature term decomposes as 

K* = ^ + cos p^- = -e(l - p 12 )- 1 ' 2 ^' + (cos <p)6' (271) 
as as 
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where is calculated using (|267|) . We have actually e = 1, see the Remark 
below. 

By substituting the first expression for sin [3 in i|267|) into RHS and compar- 
ing the leading terms of the expansions of LHS and RHS, it follows that either 
c-2 or all Ci — for i > 1, and that C2 = implies sin3#o = 0. Thus, C2 = 
means 6*o = or 7r/3, that is, the two meridian solutions of isosceles triangle 
type already discussed in Section 8.4. 

Henceforth, we shall assume C2 ^ 0. Comparison of the leading terms (of 
order 2) in yields 

3 15 
4agd = ~-d + — sin 36» 

o 04 

which combined with H261[) gives 

15 sin 3£?o 75 sin 2 39q 



~ u 16 (16a 2 + |) ' C2 512 (16a 2 + |) 2 (2?2) 

with the approximate values 

d w 0.293 sin 36» , c 2 ~ -0.014 sin 2 3<9 

Remark 55 In particular, in $262\) we have s = sgn(do) = 1, and the expres- 
sions in \272\j are, in fact, valid in the whole closed chamber Co : < 9q < ir/3. 
However, e = ±1 actually changes sign across the border meridian of two neig- 
boring chambers. 

In view of l|267|l . I|271|l we also need the expansions 

(1 - if' 2 ) 1 / 2 = d s(l + b lS + b 2 s 2 + ...) 

(l-v/ 2 )- 1 /^ J- {1 + - blS + b 2S 2 + -) 

d s 

where by writing Ck = Ck/c 2 for k > 3 we have by simple inspection 
n + 3 

b„= - c n+2 + B n (c 3 , ...,c n+ i), n>\ (273) 
n + 3 

&n = ^ — c n+2 + B n (c 3 ,...,c n+ i), n>l 

6 

where -B„ and B n are polynomials and B\ = B\ = 0. 

For simplicity, let P(yi, y%, ...) denote any (unspecified) polynomial in the 
variables j/i, except that yi = #o means it is polynomial in sin3#o and cos3#o- 
Using the notation 

sin^ = s(l + 52 s 2 +.93S 3 + ...), g k =Ck + P(c 2 ,...,c k -i) (274) 

cos<^ = 1 - is 2 + /i 4 s 4 + /i 5 s 5 + h k = P{c 2 , c fe _ 2 ) 
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we derive from (|262[1 the following formula for d n 

(n + l)d n = d b n - (dog n + 2dig n -i + ■■■ + (ra - 1)^-292), n > 1, (275) 
and from (|271|l we calculate the curvature expansion 
/C* = 2d (l + fcis + fc 2 s 2 + ...) 

where 

2d k n = (d b n + (n + l)d n - (n + 3)(n + 2)^J (276) 

^ ri+l n — 1 

- — ^(fc + l)kc k b n+2 -k + + l)d k h n ^ k 

fc=3 fc=0 

Next, let us have a closer look at the coefficients of u(s) and u'(s)/u(s), 
using (|226|) , l|264|) , l|265|> , and also at the recursive generation of the coefficients 
of a(s) using l|263|) . It follows that they are of type 

u n = P(0o, c 2 , c„_ 2 , d , ...,d n -i), n > 4 

fi n =P(9 ,c 3 ,... ,Cn,d ,...,d n -2), n> 2 (277) 

a„ = P(ai,a 2 , a„_i, ...,/i„), n > 1 

For example, 

gMi 8/fi 15/8 

12a 2 + a 10 - 10 - \/l3 

For convenience, write 

1 - cos 2a = 2a 2 1 s 2 (l + A\s + A 2 s 2 + ...), (278) 

where 

^i=2ai, A-2 = 2a 2 + a 2 - ia„, . . . , A n = P(a x , a 2 , a n ), 
and thus we arrive at the following presentation of LHS as a product of series 
LHS = 4a 2 ,d s 2 (l + A v s + ....)(! + k x s + ....)(! + u lS + ....) (279) 



From the structure of U* as a trigonometric series (|226|) in the variables 
sin<^ and cos 30, we can write 

dU* dU* 

— — = (sin w cos iii, — — = (sin 3 <£>sin30)i?2 (280) 
dip ' 06 

where 

3 I 5 • 0/1 I 5 945 • 2 
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are again polynomial scries in the variables sinip,cos39. Hence, by substituting 
the series of sin tp, cos tp, 9' , <p' , cos 39, sin 3$ into the expression 

RHS = (sin 2 tp) (-6' cos ipR x + tp' sin 39R 2 ) (281) 

equation l|269|) renders a recursive procedure for the calculation of c n and d n -2, 
n > 2, starting from c 2 and do (|272|l . 



Lemma 56 TTie coefficients c m and G? m _ 2 can be expressed as 

c m =P(0 )c 2 , d m -2 = P(e )d ; m> 2 (282) 

where P(0q) denotes some polynomial of sin 39 'o and cos3#o (generally different 
for each coefficient). 

Proof. The first step is to compare terms of order 3 in (|269|l . which renders 
the identity 

3 15 

4ald (Ai + fci + ui) = —-^d\ + —d Q cos3# 



where 



Consequently, 



a i ; C3 n j d ° C3 

Ai = 2ai, ki = — , ui = 0, di = - — 

c 2 oc 2 



15(16a + l) 
4(16a§ + l)(12a + 1) ' 



c 3 = ^/r.^r , ^ ( c ° s ^2 (283) 



15(16a + l) . 

"1 = ; o w r(cos30o)an 

12(16a2 + l)(l2a + l) v 0J 

We proceed by induction and assume that Q282JI holds for m in the range 
2 < m < n. By (|273|l - (|278j). we infer that b m - 2 ,b m -2, k m ~2, g m h m ,u m , [i m 
and a m are all of type P(9q) for m < n. Furthermore, 

nd„_ x = (^p c «+i + p (#o)) d Q (284) 
(n + l)(n + 2) 

Kn-l = T~ Cn+i + P{a ) 

12c 2 

Consider the terms of order n + 1 in equation (|269|l . By equating the coef- 
ficients of s n+1 in LHS and RHS we deduce 

4a 2 d fc„-i + P(# )do - - \nd n -! + P(9 )d (285) 

o 

Here, k n -\ and d„_i are the only coefheients depending on c„ + i, and by sub- 
stituting their expressions from l|284|) into the identity l|285|) . we deduce the 
identity 

'4 ag(n + l)(n + 2) , n + 2\ 

12 + j C " +1 = m)C2 ' 
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and consequently, 

c„+i = P{0 o )c 2 , d„_i = P(0 o )d o 

This completes the induction step, and hence l(282[) holds for all m > 2. ■ 

This settles the existence and uniqueness question for the series expansions 
1)260(1 . for each initial longitude angle 9q. Their radius of convergence is certainly 
positive (e.g. by an inductive argument showing the coefficients are bounded), 
but we shall not try to estimate the radius here. Clearly, for 9q = kir/3 the 
radius is at most w/2. 

8.5.2 Symmetries of the solution set S(po) 

In the previous subsection it was established that the triple collision solution 
set (124211 for p = po is naturally parametrized by angles 8o, namely 

6(po)={(a eo ,r; o );O<0 o <27r} (286) 

is in 1-1 correspondence with points on a circle and therefore inherits "sym- 
metries" of a circle. However, the actual symmetry group should act with or- 
bits representing the various "species' or "congruence" classes of solutions, and 
moreover, knowledge of each class suffices to generate all solutions by a straight- 
forward transformation procedure. We contend that <&p defined in (|243fl is, in 
fact, the appropriate group. 

First of all, (&p a contains the group D 3 which acts on the 2-sphere and 
represents the purely geometric symmetries. On the other hand, we have also 
seen that each solution curve T* = Tg Q corresponds to three analytic functions 
(a(s),(p(s),8(s)) in a neighborhood of s = 0, and for s < these functions 
also describe a motion approaching a triple collision as s — > CP. Hence, by 
inverting its direction we should obtain a triple collision motion emanating with 
initial longitude angle Oq + tt, which by uniqueness must be the solution in (|286fl 
labelled by 9q + ir. Consequently, in agreement with the summary of Section 
8.4.1, for each "antipodal" pair (ctg , Tg ), (ag 0+7T , Tg o+7T ) in (|286fl there is an 
analytic curve r*(s) passing through p and an analytic function a(s) so that 
(a(s),T*(s)) is a solution of the system (|221fl . and moreover, 

ag a + 7T {s) — —a(—s) and ag (s) = a(s) for s > 
n 0+ M = r*(-a) and r; o (*) = T*( S ) for s > 

In particular, the inversion operator f applied to (a(s),T*(s)) induces an 
involution 

f : (ae ,r*g o ) -> (ae 0+7T ,Tg 0+7T ) 

of the set 6(po) which commutes with the action of D3, and together they 
generate the dihedral symmetry group 

e ftj =3) 3 x{l,f} ~S) 6 (287) 
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which may be viewed as an "isotropy" subgroup of © in (|232|) . In effect, this 
divides the fundamental chamber Co (|233() in two sectors of angular width 7r/6, 
say 

C : < 6»o < tt/6 

is our reduced fundamental chamber. However, we remark that ®p does not act 
on S 2 , so Co is not a fundamental domain in the geometric sense. But solutions 
starting out in this region suffice to generate the whole solution set (|286ll using 
analytic continuation. 

The power series developments <|26U|) of the three functions a(s),ip(s),6(s), 
where Tg o (s) = (tp(s),9(s)) is the shape curve with initial direction 9q at the 
north pole, also exhibit a specific symmetry pattern which reflects the Dq- 
symmetry of their coefficients. For example, consider the "reflection" in <3p 

9 a tt/3 - 9 

which divides Co into two reduced chambers, and write 

X = cos 36*o, Y = sin 300 

a k = A k (X,Y), c k = C k (X,Y)c 2l d k = D k (X,Y)d 

where A k ,C k ,T) k are polynomials of two variables (not unique, of course, since 
X, Y are algebraic dependent). Let a k ,c k ,d k and a k ,c k ,d k be the coefficients of 
the solutions (a, ip, 9), (a, ip, 9) corresponding to initial angles 9q and 7r/3 — 9q, 
respectively. Then we have 

f a k = a kl c fe = 5fe, d k = d kl for k even 

\ a k = -ak, c k = —Ck, d k = -d k , for k odd 

Equivalently, as functions of X the polynomials A kl C k ,T) k are odd (resp. even) 
functions for k odd (resp. even). This is due to the fact that Y is invariant 
whereas X changes sign under the substitution 9q — > ir/3 — 9q. 

8.5.3 Symbolic manipulations and numerical calculation of power 
series 

The recursive scheme used in Section 8.5.1 will generate all higher order coef- 
ficients as polynomials of V = sin36>o and X = cos3#0j but the explicit calcu- 
lations involve a substantial amount of symbolic manipulations. For example, 
various types of algebraic operations, together with composition, are applied to 
power series. 

In principle, calculations involving elementary functions of power series, such 
as sin(^pi2; 1 ), can be reduced to symbolic manipulations on power series of the 
type 

OO 

(po+Pix+p 2 x 2 + ...) n = Y. P t n)xk 

k=0 
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where the n-th multinomial polynomial pj^ records the n-partitions and asso- 
ciated multinomial coefficients which can be calculated recursively with some 
effort. 

On the other hand, available computer software developed for symbolic com- 
putation have built-in procedures which effectively generate the intermediate 
power series expansions as well as recursive formulas. We have employed such 
symbolic software for the calculation 1 of ak,Ck,dk in (|260JI . for small fc, see also 
Section 8.7. For convenience, we list the first of them below (omitting the al- 
ready known C2, do, do), and we remark that the exact (or symbolic) expressions 
are growing fast in complexity as fc increases : 



C3 = C2 1^10W13) 
lib 

/ 2040762505 + 136353812\/l3 20984375(113 + 20\/l3) . „„ ,\ 

c 4 = c 2 1 — -(cos66> ) 

^ 2891425280 2891425280 v 'J 



5(10 + ^/13) , 

"i = «o 775 (cos36>o) 

lib 

, / 2(19733316 + 84347^) 302175(113 + 20^), „„ ^ 

d 2 = d \ — — - H i — -(cos6<9 ) 

I 216856896 216856896 v ' 



-5(10 + V13) , 
«i = 232 (cos36» ) 



3(28004 + 4175V13) -25(644 + 47V13) . 

a 2 = — — - H — Mcos66» ) 

2745024 2745024 v ' 

8.6 Global behavior of the shape of triple collision motions 

First we shall investigate the curvature properties of the flow consisting of the 
gradient lines of the potential function U* on the sphere S 2 (l). This information 
will be related to the curvature properties of the " flow" consisting of those curves 
r* belonging to the set l|286(l . that is, the triple collision shape curves emanating 
from the north pole Pq. 



8.6.1 Differential geometry of the gradient flow of U* 

We start with the following elementary result about curves on the unit sphere 
S 2 in Euclidean 3-space. 

Lemma 57 Let t — > p(t) be a parametrized curve on S 2 . Then its geodesic 
curvature is given by the following triple product 

^) = (| fpxp-p = pxp'.p" (289) 

lr The symbolic and numerical calculations were performed in 1995 by Chee-Whye Chin, an 
undergraduate student at U.C. Berkeley, using Mathematica software. 
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where (as usual) s is arc-length, p —-^p and p' = ^P- 

Proof. The unit tangent vector t* = p' points in the positive direction of the 
curve, and p" is the curvature vector in 3-space. With u* = p x t* as the 
normal vector field along the curve, {t*,u*} is a positively oriented frame of 
the sphere. By definition, the geodesic curvature vector in the sphere is the 
orthogonal projection of p" into the tangent plane, namely 

p" - (P" ' P)P - K g u* 

where the coefficient K g is the (scalar) geodesic curvature. Clearly, K g equals 
the triple product in (|289f) . ■ 

Next, we turn to the gradient field VC/* on the sphere, whose integral curves 
will be referred to as the gradient lines of U* . Their geodesic curvature will 
be denoted by K g . Until further notice there is no restriction on the mass 
distribution, and we use the expressions i|171|) . (|175|l and (|176|l for U*, the 
vector function B(p) and the gradient Vf* field, respectively. 

Lemma 58 The geodesic curvature of the gradient line of U* passing through 
p is given by the following triple product 

K g (p)= * / xC-p (290) 

where B = B(p) and 

3 , 

C = C(p) = ^ i—g \p x B ■ (p x bjj b, (291) 

»=i p - hi 

Proof. Consider a gradient line parametrized by arc-length, s — * p(s). We can 
write 

VU*(p) - /(s)p'(s) 
where f(s) = |VJ7*(p(s))|, and hence by lfT7^|) 

, = B - (B p)p 

P |VC/*| 

On the other hand, = pxp'-V?7*=pxp'-B and by differentiation 

p x p" B = p x p' • 4-B (292) 
ds 

Hence, by substituting the expression 

B = |V(7*|p' + (B-p)p 
into the left side of QSGfy . we obtain 

IVP-lpxpr-pP-PXp-.^B-E^.^B (293) 
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Finally, we calculate from (|175fl 



ds 



-B 



^^(MpHBxp)^ 



and by substituting this into (|293[l it follows from Lemma 1571 

3 v-, (pxB).(pxb.) 



K g = pxp'p" = 



(pxB-b,) 



This expression can be rewritten as (|290|l . ■ 

Problem 59 Regarding K g as a function on S 2 , determine the curves defined 
by the condition K g {p) = 0. 

Remark 60 Note that B(p) and C(p) are vectors in the xy-plane in the Eu- 
clidean model M = M , and the function K g on the sphere S 2 : x 2 + y 2 + z 2 = 1 
is undefined precisely at the critical points of U* , namely for z > these are the 
points £>i,0 < i < 3, and the minimumspoint (physical center) pQ . From the 
triple product formula 129 (J^ it follows that K g (p) vanishes on the eclipse circle 
(z = 0), whereas for z > K g (jp) vanishes if and only i/B(p) and C(p) are 
linearly dependent. 

Henceforth, we shall retain our assumption of uniform mass distribution, 
and a deeper understanding of the function K g will be achieved. In fact, in this 
case the above problem has a simple solution, as explained at the end of this 
subsection. 

By assumption, 



3' h 



V2 
3 ' 



hi ■ b , 



- (i, j > 1 and i ^ j) 



and we introduce the three distance functions 
5 { = Si{p) = 



p-b 

which are algebraically related by 

<5l + <$2 + ^3 = 6 > 
due to the identity ^ p • b^ = p- 53 bj ~ anc ^ 



V2y 1 — p • bj, i > 1, cf. Section 6.1 



(294) 



(295) 



P • b, = 1 



(296) 
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Let us write 



B( P ) 
C(p) 



Y = ~y~ (( ei ~~ e3 )^ x + ( e2 ~~ e ^ 2 



V2 



1 



1 



()Z ^/ibi = -((/i-/ 3 )bi + (/ 2 

2= 1 



/s)b a ) 



where 



1 3\/2 

73, /< = -s-(p x B(p)) • (p x bj 

0; 0,- 



(297) 
(298) 

(299) 



To simplify our notation we denote products (monomials) of the functions 



Si by 



S a ,b,c = S*6 2 6l, 



(300) 



where a, b, c are nonnegative integers, and the alternating polynomial generated 
by the monomial (|300|l is 



A a ,h,c 



rd rd rfl 

d l d 2 d 3 

cfa rfc 

dj. d 2 d 3 

5£ S c 2 S c 3 



<r£S 3 



sgn{(T)5 a ab 



(301) 



where S3 is the permutation group of {^1,^2,^3} acting on monomials in the 
obvious way. In particular, the basic alternating polynomial is 



A = A 1 2 



1 

Si 



1 

S 2 



s i s 2 



1 

s 3 
si 



(S 1 -S 2 )(S 2 -S 3 )(S 3 -S 1 ) 



(302) 



On the other hand, the symmetric function generated by <5 a .fc,c, where we may 
assume a > b > c, is the smallest ^-invariant sum 



S a , b ,c = S i S j S k = SI 



b,c 



(303) 



containing <5 a ,&,c- In particular, 6*2,0,0 = 6, by Ij295|l . 

Note that £ a ,&, c is unchanged, whereas A a ,f,, c may change sign, when a, b, c 
are permuted. The alternating polynomials can be decomposed as a product of 
the basic alternating function H3U2|) and a symmetric function, for example 



A), 2, 4 = (£2,1,0 + 25i,i,i)A 

A-0,1,5 = (83,0,0 + 82,1,0 + Si,i,i)A 

A), 1,6 = (£4,0,0 + 63,1,0 + £2,2,0 + S , 2,i,i)A 

Aq,3,6 = (£4,2,0 + £4,1,1 + 263,2,1 + 35 2 ,2,2)A 



(304) 



The induced action of 6 3 on the triples {ei, e 2 , 63} and {/1, / 2 , fy} is covari- 
ant with the action on {5i, 62, S3}. Certainly, the above coefficients in (|299|l 
are simple rational functions of the S'jS, and now we prove a similar statement 
for the /[s, as follows. 
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Lemma 61 As a rational function 0/61,62,63 

fl = T (—351,3,6 - 3<5i, 6 ,3 + 2# .6,6 + £3. 6j3 + £3.3,6 + #1,5,6 

06,6,6 

+#1,6,5 — 2^ 2 ' 6 ' 6 — 2^ 3,5 ' 6 _ 2^ 3 ' 6 '' 5 ^ 

and fa (resp. fa) is obtained /rom fa (resp. fa) by cyclic permutation, a — > b — > 
c — * a, 0/ the indices 0/ each monomial # ,&, c . 

Proof. Use the identities H296[) and substitute the expression H297|) for B into 
the formula (|299|l for fx, namely 



/i = ^(B-bi-( P -b 1 )(B. P ; 

#1 v 



Then one obtains the above rational expression for fa by straightforward calcu- 
lations, and by symmetry it is also clear that fa and fa are obtained from fa as 
claimed. ■ 

Now, turning to the formula (|290|) for the curvature function K g and insert- 
ing the expressions 1)2971) . 12981) . we write 

B(p)xC(p).k=^U 

where k is the unit normal vector of the xy-plane and A is, by definition, the 
function 

i(#i,# 2 ,# 3 ) = {ex fa - e 2 fa) + (e 2 fa - e 3 fa) + (e 3 fa - ex fa) (305) 
When the products eifj are calculated using the above lemma, for example 

61/2 = 7 ( — 3(53,1,3 — 3#o,l,6 + 2^3,0,6 + #3,3,3 + #3,1,5 + #0,3,6 

06,6,6 

x 1 x 1 x 1 s \ 

+ #2,1,6 — 2#3,2,6 — 2#3,3,5 — ^ 2' 3 ' 6 >> 

the expression l)305|) may be written as 

A = -T— (3Ao,l,6 + Ao.3,6 + ^1,3.5 + A 1)2 , 6 ) (306) 
06,6,6 

Finally, substitution of expressions from l|3U4|) into (|306[) yields 

l(#l,#2,#3) = -5(#i,# 2 ,# 3 )A 

where 

S{6x, #2, #3) = 775 rs- [3(5 < 4.o,o + 5*3, 1.0 + 52.2,0 + 62,1,1) + £4,2,0 + 64.1,1 

^1,1,1) 

+ 253,2,1 + 3(51, 1,1) + 25i, 1,1^2, 1,0 + 3(51, 1,1) + 5*1, 1,163,0,0] 

In summary, we have established the following proposition, where the factor 5 
of K„ is always positive ! 
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Proposition 62 The geodesic curvature function of the gradient lines ofU* on 
the unit sphere S 2 is given by the product 

\/Q S 

K 9 (P) = -— — p=(x,y,z), z> (307) 
18 | vU*(p)\ 

Observe that the 2)3-chambers of the (upper) hemisphere of S 2 are defined 
by inequalities Si < Sj < 6k, and therefore A, and hence K g as well, has constant 
sign in each chamber. For example, the fundamental chamber l|233|) is given by 

C : S 2 < 8 1 < S 3 

and here K g > 0. The meridians which are the walls of the 2)3-chambers are 
defined by relations of type Si = Sj; these are the zero set of the function A. 
Together with the equator circle they are the curves where K g vanishes (or 
is undefined), and this also solves Problem (in the case of uniform mass 
distribution). 

It is easy to visualize the gradient flow on the 2-sphere. For example, in the 
interior of the spherical triangle Co the flow has the vertex po as source and 
converges towards the vertex b2, with positive curvature everywhere. The flow 
in C is illustrated in Figure 9. 



8.6.2 Geometry of the triple collision shape curves 

It is possible to draw qualitative information about the family 6(po) of shape 
curves by relating it with the geometry of the gradient flow of U*. By symmetry 
it suffices to consider those curves Tg g starting out in the chamber Co, that is, 
< #o < § > and we observe that the boundary meridians 9q — and 9q = 7r/3, 
emanating from the north pole po towards the equator, are themselves both 
shape curves and gradient lines. So, the question is what one can say about the 
shape curves in the interior of Co? 

A rough description goes as follows. In Co there are three different "flows" 
of curves emanating from po, namely the shape curves, the gradient lines and 
the meridians (9 constant). At each point p = r*(si),si > 0, the shape curve 
r*(s), s > si, is "trapped" between the gradient line and the meridian through 
p. Being positively curved one may imagine the shape curves arising by grad- 
ually bending the meridians towards the gradient lines by means of a "force" 
field directed eastward, see Figure 12. 

To be more precise, we shall focus on five properties as stated below. For 
this purpose we introduce two angular functions (3{s), 7(5) as follows. Namely, /3 
is the angle between the meridian and T* , as defined in (|267|1 . It is the oriented 
angle from ^ to the velocity vector t* = ^T*, and 7 is the angle from ^ to 
the gradient vector 

„„, dU* d 1 dU* d 

VZ7 = 1 

dip dip sin 2 <p 09 39 

lr7ml , d . Id. 

= VZ7 (cos7— + SU17- — 

dip sin <p o9 
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We also recall the role of the inclination angle a(s), which is not directly related 
to the geometry of the spherical curve T* itself, but to the associated moduli 
curve f . However, the curvature equation from 1|221[1 can now be stated as 

2CT7C* sin 2 a = |V(t/*)| cos(^ + f3 - 7) (308) 
y 2 

and hence relates all three angles with the curvature of T*. 

Now, we contend that the following properties are valid for T* = Tg o , < 
#o < § j at least until T* leaves £q the first time (but not necessarily later) : 



(i) < a < J, for s > 0. In particular, the curve s — > T*(s) has no cusp 
singularity for s < tt/2. 

(ii) The spherical coordinates tp(s),8(s) of T*(s) are strictly increasing 
functions of s. 

(iii) < 13 < 7 < f , for s > 0. 

(iv) The geodesic curvature K.* of T*(s), s > 0, is nonnegative. 

(v) T* leaves the chamber £ by crossing its boundary arc (e 3 ,b 2 ) on the 
equator circle. 



In order to verify these statements one may proceed as follows. First, note 
that < 7 < 7r/2 follows from the fact, due to (|228|l and Remark 021 that 
^g— > and ^ijg- > inside £0, and moreover, the gradient lines emanate from 

Po with 7 = and approach the binary collision point t>2 with 7 = 7r/2 in the 
limit. This also explains why property (v) follows from (iii), and using (|2"H7l and 
H3()8(l we also readily deduce properties (ii) and (iv) from (iii). Thus, we are left 
with the statements (i) and (iii), and let us first establish property (iii) (using 
property (i) if necessary). 

Observe that (3 > by l|267l) and Remark [S5j But (3 = for some s would 
imply 6 = 0, ip' = 1, and hence by (|271(l IC*(s) would vanish. However, with p 

in the interior of £0 we also have (p) 7^ and then 7 > in the right side of 
the identity l|3U8|) . Consequently, A^(s) 7^ and this contradiction shows (3 > 
must hold. Again by l|308() . JC* is positive as long as (3 < 7, and this certainly 
holds for small s since 

IC (0) = d 



16 (16a 2 + f ) 

We claim that (3 < 7 holds (at least) until T*(s) leaves the chamber. To see this, 
suppose we had = 7 for s = si, j3 < 7 (resp. f3 > 7) for s < si (resp. s > s\) 
and s close to s%. Then T* would be tangent to the gradient line at p = T*(si) 
and is (locally) lying on the "upper" side of it, hence K,*{si) > K g (p) > 0. This 
contradicts the fact that /C*(si) = 0, by (|308|l . 
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Remark 63 Property (Hi) implies that U* increases along the curve T*(s),s > 
0. In general, the event = 7 means T* is perpendicular to the level curve of 
U* , and by ISUWll this can happen for two reasons, namely i) K* vanishes or ii) 
r* reaches a cusp. We can rule out the second case due to property (i), but only 
up to the first crossing of the equator. 

Finally, we turn to property (i). From the relations 

?^ = 17 >0 ' ^f = p(s)cota, cf. (EH, 
as at at as 

we deduce < a < ir/2. In fact, a = ir/2 (and ds/dt = 00) only at the vertex 
h>2, and our claim is that a = only holds at po = T*(0). 

It is certainly evident from the numerical analysis of the shape curves (cf. 
Table 1) that a > for s > 0, at least until T* leaves Co- Indeed, using numerical 
data and a continuity argument we can establish the uniform lower bound a > 
0.18 for 7r/4 < ip < 7r/2. However, we shall also explain an alternative and more 
qualitative approach to settle the problem. 

To show a > holds, let assume the contrary and recall the geometric 
arguments in the setting in Section 7.3, where we regarded T*(o~), a = s/2, as 
a curve on the sphere S 2 (l/2) and C(T*) C M denotes the cone surface with 
the induced Euclidean metric (|204|l . The associated moduli curve s — ■> T(s) 
lies in this surface, and assuming the first cusp occurs at a = o\ we consider 
the Euclidean sector < a < o~\ bounded by the rays cr = and a = u\. By 
our assumptions, there is a bijective correspondence [0, s{\ < — > [0, o~\] between 
the arc-length parametrizations of the moduli curve T and T*, and a > for 
< s < St. 

The curve T(s) starts out from the origin and its radial distance p is increas- 
ing. It has the positively oriented moving frame {t, t]} of (|206|l . where r (resp. 
77) is the unit tangent (resp. normal) vector. By 1)2111) and a well known for- 
mula for the curvature of curves in the Euclidean plane, the curvature of T in 
the above sector can be expressed as 

^ = --^-(\nU), ( = o + a , cf. Figure 8 (309) 
ds 2 dr) 

Towards the point r(si) the curve T becomes tangential to the boundary ray 
a = (7i, that is, the angle a decreases to zero. Therefore ^| vanishes, by l|2U7|) . 
and hence 

d( 

— < as s — > si 
ds 

Moreover, the frame {t, 77} approaches as s — > si, so we also deduce 

arj u(a) 

and hence u(a) = U*(T*(a)) is decreasing at u\, or possibly u'(ai) = and 
hence T* is tangential to the level curve of U* on the sphere. However, U* is 
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actually increasing towards the point T*(ai) by Remark Ifi3l f which applies here 
since there is no cusp for a < a\). This contradiction rules out any occurrence 
of cusps inside Co- 



8.6.3 Final escape limiting behavior of the shape curves 

As an interesting example, recall the time parametrized meridian solution r*|(i) 
from Section 8.4.2. The first cusp appears at t\ sa 10.4, with colatitude tp 1 ~ 
107.5°. The next cusps occur roughly at times t 2 pa 435, t 3 pa 162400, t 4 pa 5T0 7 , 
is pa 1.24T0 9 (with due regard to numerical instability) and they appear to 
be approaching §3 as a final limit of the shape curve. The behavior of the 
curve resembles a damped oscillation converging to its "stability" point §3 as 
t — > 00. Let us refer to this limiting behavior as irregular ; namely the limit 
shape p is reached through converging cusps and a'(s) has no limit, as in the 
above example. Such a behavior at the final escape at infinity is, however, not 
necessarily related to the fact that the shape curve is a triple collision curve in 
the other direction. 

Thus, one may consider more generally the limiting behavior of moduli 
curves T(t) of three-body motions as t — > 00, assuming the limit shape p exists. 
In the irregular case, however, we do not claim that p is necessarily a central 
configuration (that is, a critical point of {/*), although this is rather likely. On 
the other hand, if p is a central configuration, we claim that it is an Euler points 
dj, and moreover, T* is not confined to the equator circle. 

We define the final limiting behavior to be regular if the final shape p = 
r*(s2) is a central configuration, where a(s) > for s 2 — e < s < s 2 and 
a(s 2 ) — 0. Moreover, T* is confined to the equator circle (and hence the 3- 
body motion is collinear) if p = e^ As in the case of triple collisions, a useful 
tool in the study of such limiting behavior is again the system ODE* (|221[1 . 
whose solutions are pairs (a(s), T*(s)). Then the fact that p — ► 00 as s — ► S2 
is expressed by the divergence of the integral l|238|) . and moreover, the system 
<|221|) itself imposes the condition that p must be a critical point of U*. 

A closer study of the above regular solutions (a(s),T*(s)) near the final 
limit may proceed in the same way as we studied triple collision shape curves 
in Section 8.5. For convenience, let us translate the arc-length parameter, s — > 
s — S2, and consider the power series expansions of a(s) and T*(s) at s = 0. The 
calculations are similar to the triple collision case worked out in Section 8.4.1 
and 8.5.1, but this time ct'(s) converges to the negative root of the polynomial 
in (|249jl . namely 



VT3+1 WH85 + 1 , . 

a = - pa -0.575 69, b Q = — pa -0.9856 (310) 

8 8 

Thus, in the collinear (Euler) case with the limit shape T*(0) = §3, with r*(s) an 
arc-length parametrization of the equator circle near S3, the differential equation 
l|245(l in the case i = 2 has a unique solution a(s) with a(0) = 0, a'(0) = b . 
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Next, for the final limit shape of Lagrange type we consider the following 
(singular) initial conditions 

r*(0) = p ; a(Q) = 0, a'(0) = a (hence a(s) > for s < 0) 

which define a family of analytic solutions (a(s),r*(s)) of the system (|221f) . 
The calculations are similar to the triple collision case, with ao equal to the 
positive number in l|25U|) . worked out in Section 8.5.1. Therefore, we leave it 
to the reader to modify these calculations and perhaps establish the same kind 
of analytic uniqueness, namely that solutions are parametrized by the terminal 
angular (longitude) direction 6 of T*(s) at po, cf. l|2tj(J[l . In particular, by 
reversing the direction of the curve segment {r*(s),s > 0} one obtains the 
shape curve with terminal direction 9q + it. 



8.6.4 More about the asymtotic behavior at triple collision 

Finally, we turn to the asymtotic behavior, in terms of the time parameter t, of 
a triple collision taking place at t = 0. Let T(t) — (p(t),T* (t)) be the moduli 
curve of such a motion, with F*(0) = po or ej, and write p, = U*(T*(0)). Then 
it is a classical result, dating back to the work of Sundman and Siegel, that (for 
any energy level h and mass distribution) 

p(t) ~ Kt 2/3 as t -> 0, K =(^/i) 1/3 (311) 

and moreover, the total kinetic energy is asymtotically dominated by the "change 
of size" , in the sense that 

T(t) ~T? = -p{tf - -« 2 t~ 2/3 ~ ~~~7~\ (312) 
W 2 Fy ' 9 p(t) y ' 

In particular, the residual kinetic energy 

1 r/s 

T° = -p{t) 2 {-Y=T-T? (313) 

due to the "change of shape" must be of lower order in t than that of T p , in 
the sense that T a /T p — * 0. However, this does not exclude the possibility that 
T a — > oo as t — > 0. We claim, however, that T a 0, and we shall present the 
following (rather heuristic) argument for this. 

Using the relationship (|219J) between p(s) and a(s) and the series expansion 
of a(s) we obtain an expansion of type 

p(s) = p s p (l + ns + r 2 s 2 + ...), p = 1/20/(0) 

Hence, for suitable nonzero constants kq and k s 

s(t) ~ KoP (t) 1/p = Ko p(t) 2Q ' (0) ~ Ks i 4a ' (0)/3 = Ks t e , 
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where the exponent e depends on the two types of triple collision, namely by 

(i) T* (0) = po : e = ^a = ~ (%/l3 - 1) « 0. 434 26 (314) 
3 6 

(m) T*(0) = e, : e = ~6 = -(-VuM- 1) » 0.98079 
3 6 5 

Now, as is the case of p(t), let us assume differentiation of s(t) also commutes 
with taking asymtotic limit, namely s(t) ~ en s t e ~ l . Then by Ij313(l 

!-(.) ~ i(», e ,V/—.. . { £ <§,= (315) 

8.7 Numerical solutions of triple collision motions 

We shall describe a modified approach to provide numerical C 1 -data for the 
1-parameter family 6(po) of shape curves representing non-collinear triple col- 
lision motions, under the standing assumption of equal masses m, = 1/3 and 
zero total energy. Recall that these curves T*{s) arise from solutions (a,T*) of 
the system l|221[) of ordinary differential equations, where T* = Tg starts from 
the north pole po on the 2-sphere with initial (longitude) direction 8 , and the 
(radial inclination) angle a > has the initial value a = 0. The basic idea is 
to obtain numerical data close to the initial point, by the analytical method, 
and use them as the initial data for the remaining integration by means of a 
Runge-Kutta method. With some more efforts we believe it is possible to settle 
Coniecture l50l bv carefully combining numerical analysis and theory along these 
lines. The following numerical analysis serves at least to illustrate the geometry 
of those shape curves. 

8.7.1 Outline of a numerical approach 

As usual, (9,ip) are the spherical coordinates of the 2-sphere M* ~ S 2 , with 
ip = at the initial point po- By symmetry (as explained earlier) we need only 
consider curves whose initial direction 9q lies in the interval < 9q < 7r/3, and 
moreover, we may as well use tp to parametrize each curve since <p increases with 
the arc-length s. Thus, let 

r^o = l(<P,O0 o (<p)),tp > 0; lim 0e o (<p) = O j 

denote the shape curve with initial angle 9o, and let ag (ip) denote the corre- 
sponding angle a as a function of tp. We shall compute the values of 9e , j^9e 
and ag for < <p < tt/2 and < 9$ < 7r/3. 

Elimination of the arc-length parameter s in the system l|221f) is achieved 
by using the third equation to rewrite the first two equations with tp as the 
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independent variable. In standard (explicit) form the new system reads 



/ 1 A , ■ 2 7^2 cot a, dU* dU* nt . 

1 rdU* , <9t7* 



L esc 2 a esc 2 i/3 — (4J7* cot ip + — — csc 2 a)#' (316) 
2(7* ou uip 

BJJ* 8TJ* 
+ (^7T csc2 a ) ' 2 - ( U * sin 2l P + csc2 « sin2 VW' 3 } 

09 Oif 

where a', 9', 9" means differentiation with respect to (p. 
We consider the power series expansions at tp = : 

9g Q (<p) =9 Q + p(f + f 1 p + f 2 p 2 + ...) (317) 
a e ((fi) = <^(.9o + gnp + g2(fi 2 + •••) 

whose coefficients can be calculated recursively as functions of 9o by the method 
of undetermined coefficients. This is similar to the calculation of the expansions 
in (|2TU|) using the system ODE*. Let 

<W 4>), ^flff (*o (sis) 

be the polynomials in <^ of degree n+1 (resp. n for the third polynomial), where 
the first two are obtained by substituting the calculated expressions for fi,gi 
(as functions of 8q), < i < n, into (|3 1 Tf> and truncating higher order terms. 
The last polynomial is the derivative of Og^((p). Then for sufficiently small 
<p the "true" functions 9g (<p), ag (<p) and -^Sg g (ip) — 9' eQ {ip) will be closely 
approximated by the polynomials (|318|l . and the approximations can be made 
arbitrarily accurate by taking n sufficiently large, i.e. by computing enough 
coefficients fi,gi. 

Given a value of 9q with < #o < 7r/3, we fix a small ip and compute 
the polynomials l|318|l . to be regarded as approximations of 9g (tp), ag (ip) and 
9'g on the interval < ip < ip . In particular, the values of the three 
polynomials at <p — (p will serve as (approximate) initial values for the functions 
a, 8 and 9' , whose further development on the interval (p < p < ir/2 is governed 
by the system (I316|l . This allows us to obtain numerical solutions for 9 and a 
on this interval, using any of the standard iterative methods. Pieced together, 
these data furnish us with the C 1 -data of the triple collision shape curves Tg 
within the interval < ip < tt/2. 

8.7.2 C x -data for a selection of triple collision motions 

As in Section 8.5.3 we perform symbolic computations to calculate successively 
the coefficients fi and gi of the expansions (|317|l . As before, these are trigono- 
metric polynomials of 30o, namely polynomials of sin3#o and cos3#0) and we 
have calculated the exact expressions for i < 9. Beyond that they tend to be 
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rather untractable in their exact form. The exact expressions for the first few 
coefficients are listed below for the sake of reference : 



5(10 + 713) . 25(113 + 20713) . 
/o= 232 Sm3g "' /l = 53824 ^ 



_ 15(544702206 + 58374421V13) . 
h 201243199488 Sm 3 ^ 



5875625(1390 + 313V13) . „„ „„ 
+ 201243199488 ^ ocoseft ^ 



_ 2(1909168577687- 51730231240V13) . 
h ~ 1318947929444352 Sm6 * 



1422740625(17969+ 4520V13) . „„ 

sm Wn COS 6t>n 

1318947929444352 



-15(l + 3y/l3) 

Qi = COS 30n, 

y 1856 ' 



543509 + 352091V13 -8925(49 + 31V13) „_ fnnn . 

qo = 1 cos60 n (320) 

y 29280256 29280256 v ' 



10(253093537+ 56946315713) 
53 134162132992 C ° S 3 ^ 



-3525375(893+ 359713) 

- cos 36*0 cos 6#o 



134162132992 

The coefficient go equals the constant ao in (|266J) and hence is omitted here. 
To illustrate the (decreasing) magnitude of the coefficients of the trigonometric 
polynomials fi , <?i we list a few approximate expressions 

/i w (8.6)10- 2 sin66» , f 3 ~ (2.6)10~ 3 sin 6(9 + (3.7)10~ 2 sin60 o cos6(9 , 
gi re -(9.6)10" 2 cos36» , g 3 ~ (3.1)10~ 2 cos3<9 - (5.7)10" 2 cos36> cos66'o 

In this way one obtains the approximating polynomials Ij318(l for n = 9, say. 
Hence, to obtain approximate numerical data for the solutions T* 6o of the system 
l|316|) . we have chosen <p = 0.05 and 9q — /c7r/300, k = 0, 1, 100, and we have 
computed the numerical solutions using the Runge-Kutta method. These 101 
solution curves outline the general behavior of the shape curve 1^ of triple 
collision motions parametrized by the initial longitude angle 9q- We refer to 
Table 1, Table 2 and Table 3 which list the calculated values of a$ (ip), 0e o {<f) 
and 2^9g (ip) for 6q = kn/30, k = 0, 1, .., 10, and for 6 different values of <p. All 
angles are measured in radians. 

Table 1 : Inclination angle ag (ip), for 9q — fc7r/30 
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k \<p 


7T 

4 


3tt 

8 


15tt 

32 


63tt 
128 


51 2 


2048 





.1947 


.2292 


.2072 


.1926 


.1883 


.1871 


1 


.2132 


.3470 


.5938 


.7243 


.7705 


.7835 


2 


.2464 


.4554 


.7917 


.9694 


1.0379 


1.0586 


3 


.2742 


.5170 


.8895 


1.0914 


1.1747 


1.2014 


4 


.2943 


.5542 


.9465 


1.1643 


1.2597 


1.2922 


5 


.3083 


.5780 


.9827 


1.2121 


1.3179 


1.3566 


6 


.3181 


.5938 


1.0067 


1.2247 


1.3599 


1.4054 


7 


.3248 


.6043 


1.0226 


1.267? 


1.3904 


1.4433 


8 


.3292 


.6110 


1.0328 


1.2816 


1.4116 


1.4723 


9 


.3317 


.6147 


1.0385 


1.2898 


1.4243 


1.4917 


10 


.3325 


.6159 


1.0403 


1.2925 


1.4285 


1.4989 
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Table 2 : Longitude angle 6g a ((p), for 9q — kir/30 



k \<p 


7T 
4 


3-rr 
8 


IbTT 

32 


t)37T 

128 


2bb7r 
512 


10237T 
2048 























1 


.23558 


.40978 


.66916 


.77595 


.80954 


.81866 


2 


.42158 


.61693 


.83930 


.92403 


.95180 


.95968 


3 


.55924 


.73215 


.90859 


.97447 


.99697 


1.0037 


4 


.66448 


.80830 


.94747 


.99907 


1.0173 


1.0231 


5 


.74914 


.86481 


.97349 


1.0137 


1.0283 


1.0333 


6 


.82059 


.91025 


.99302 


1.0236 


1.0350 


1.0391 


7 


.88344 


.94901 


1.0089 


1.0310 


1.0394 


1.0426 


8 


.94075 


.98370 


1.0226 


1.0370 


1.0425 


1.0448 


9 


.99475 


1.0160 


1.0352 


1.0423 


1.0450 


1.0462 


10 


1.0472 


1.0472 


1.0472 


1.0472 


1.0472 


1.0472 



Table 3 : Values of -£z9e (<p), for = /br/30 



k \<f 


IT 

4 


3tt 

8 


15tt 

32 


63t7 
128 


255tt 
512 


2048 























1 


.29933 


.63143 


1.2517 


1.7182 


1.9466 


2.0210 


2 


.40322 


.60768 


.99427 


1.3887 


1.6602 


1.7708 


3 


.39183 


.50135 


.76959 


1.0998 


1.3930 


1.5427 


4 


.34216 


.40251 


.60029 


.87390 


1.1710 


1.3623 


5 


.28289 


.31737 


.46599 


.68801 


.97142 


1.2023 


6 


.22295 


.24306 


.35361 


.52754 


.78004 


1.0406 


7 


.16476 


.17645 


.25526 


.38365 


.58941 


.85626 
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.10857 


.11502 


.16580 


.25043 


.39595 


.62716 
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.05389 


.05675 


.08166 


.12369 


.19908 


.33781 
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8.8 An outlook on the general case 

Finally, let us briefly consider the more general case of non-equal masses and/or 
non- vanishing total energy h. Namely, by Theorem 021 the system ()221l) must 
be replaced by 

' f = ^ + i^(l + 2cota^ln({/*)) 

K; S in 2 a= !^^ha(t/*) (321) 
l = (^ + (sinV)(f) 2 

where u(s) — U*(T*(s)) and T*(s) is the arc-length parametrized shape curve 
12181) on the unit sphere <S' 2 (1). Now the size function p in the moduli space M 
appears explicitly, so the system involves the two "auxiliary" functions a(s), p(s) 
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which are still related by 1)2071) and 1)2081) . Their initial value at triple collision 
is a(0) = p(0) — 0, and the shape curve T* emanates from the physical center 
Po = r*(0), namely the minimumspoint of U* on the northern hemisphere. Due 
to the space-time scaling symmetries of the Newtonian equation (cf. Chapter 
7), for h ^ there are essentially only two cases, h > and h < 0, and our 
case h = may be viewed as the limiting case between negative and positive 
energies. However, for /i^Owe may scale and assume h = ±1. 

We point out the open problem of finding the appropriate version of Theorem 
G (or Gi) in the two cases h = ±1 or when the masses are unequal. For this 
purpose, it is natural to try first the following two special cases. 

• h = ±1 and equal masses (i.e., mj = 1/3). Many results in Section 8.5 
still apply and there are symmetries as before, e.g. it suffices to consider 
shape curves whose angular direction at the north pole is in the range 
< Go < 7r/3. However, the initial value problem is "essentially" singular, 
in the sense that the solutions of (|321|l are singular at s — 0. For example, 
if we assume a series expansion 

a = a s q (l + ais + ..), p = p Q s p (l + ns + ..) 

then we would have 9=1 and p would have the leading exponent p = 
l/2ao- Moreover, from the first equation of (|321|l it follows that do has 
the value from (|266[1 . but this equation also tells us that T*(s) is singular 
at s — 0. 

• h = and the masses are not equal. The system (|321ll is the same as l|221|) . 
Note that the mass distribution {mi] affects the system ()321l) solely via 
the potential function U* , but the trigonometric series development of 17*, 
similar to that in Section 8.2, remains to be done for non-equal masses. We 
also seek a convenient coordinate system on the sphere near the point po . 
In Section 6.6.1 we actually worked out a series expansion of U* centered 
at p , involving coefficients which are symmetric functions of the masses, 
but here we rather need its spherical polar coordinate version. 

On the other hand, in Chapter 6 there are expressions for U* in terms of 
spherical polar coordinates ((p, 9) centered at the north pole N '. Therefore, one 
approach is to find the coordinates (<p, 9) of po in this coordinate system and 
then use spherical trigonometry to determine the transformation from (ip, 9) to 
spherical polar coordinates (<p,6) centered at po- For this purpose, we recall 



cosy = cf. (tT77a (322) 

m 

and the value of 9 (which depends on the choice of zero meridian, 9 = 0) can 
be determined from the longitude differences uii = ±(i 9 — with < u>i < it, 
i = 1,2,3, where 9i denotes the longitude angle of the binary collision points 
hi. To this end, consider the spherical triangle with vertices A/", po and and 
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note that cjj is the angle at the vertex A/". The arc opposite to AT, connecting 
Po and hi, has length 



2a t = 2arccosW- — - — = 2arccos Jl - l -— cf. 1171 . fTTOll 

V 1 — Tit V (1 - mi)m 

and the two other sides have length ir/2 and (p. Hence, by the spherical cosine 
law we deduce the formula 

. . 1 _ Iriij 
C0S(2er. t ) 1 (l-m 4 )m , QOQ ^ 

coscj, = — = r^==- (323) 

sm (p . / 1 — 3 ™ 



Knowing the expansion of U* as a trigonometric series in terms of the angles 
ip, 9, we believe the initial value problem at s = for the system (I221f) can be 
investigated in the same way as we did for equal masses. But this time we 
expect regularity issues to depend crucially on the mass distribution. 

NOTE : The figures belonging to this memoir are available at the website 
http://www.math.ntnu.no/~eldars/allfigures.pdf 
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